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ABSTRACT 


We derive the closed form expression for the probability of bit error in direct detection, 
dense wavelength division multiplexed (WDM) fiber optic systems employing 00K 
as the modulation technique and single-cavity Fabry-Perot (FP) filters in the receiver 
as demultiplexers. The expression is derived in the time domain using the impulse 
response of the single-cavity Fabry-Perot filter and the complex baseband equivalent 
of the received dense WDM optical signal. The two are convolved to produce the FP 
filtered output signal s(t). We then integrate 7l\s(t)\ 2 over one bit period, where TZ 
is the responsivity (A/W) of the photodetector following the FP filter in the receiver 
structure. This integral is the deterministic X of the decision variable Y where 
Y = X + N. N is the postdetection thermal noise (amplifier generated), a zero mean 
Gaussian random variable with variance N 0 T where N 0 is the noise current spectral 
density (A 2 /Hz). Both X and N are combined into an expression for probability of 
bit error. A limited case of the complete model is assumed, and probability of bit 
error graphs are generated. 
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I. INTRODUCTION 


Direct detection optical frequency division multiple access (FDMA) networks 
are increasingly becoming an attractive alternative to coherent optical FDMA net¬ 
works [1]. One of the primary reasons is that noncoherent systems do not require 
expensive synchronization circuitry for proper operation. Also, present optical filter 
technology allows designers to closely pack the channels in frequency, resulting in 
dense wavelength division multiplexed (WDM) systems that can provide aggregate 
bit rates of many terabits per second (IT b/sec = 10 12 bits/sec) [2]. The possible 
uses of dense WDM systems are many (local area networks, undersea surveillance, 
etc.), but one can easily see the enormous economic benefit and importance of being 
able to transmit aggregate bit rates of terabits per second on a single fiber without 
the economic burden of expensive synchronization circuitry imposed by a coherent 
system. In this thesis, we derive the complete closed form expression in the time 
domain for the probability of bit error of dense WDM systems employing OOK as 
the modulation technique and single-cavity Fabry-Perot (FP) filters as channel de¬ 
multiplexers. In our derivation, we make no simplifying mathematical assumptions 
and use the impulse response of the single-cavity Fabry-Perot filter, which is an in¬ 
finite sum of delayed impulses whose intensities decrease geometrically. The results 
of our work are presented in Chapters II, III, and IV. A detailed derivation of the 
deterministic signal component of the decision variable appearing at the output of 
the integrator of the Channel 0 (Channel of Interest) receiver appears in Appendix A. 
See Fig. 1 for a sketch of this receiver. Appendix B shows how the complete model 
is reduced for a limited case and how probability of bit error calculations are made 
for this limited case. Appendix C shows the strategy and programs used to obtain 
the probability of bit error for four values of free spectral range-bit period product 
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and a given range of signal-to-noise ratios. It is interesting to note at this point 
that, although the model presented in this thesis is strictly derived, complete, and 
done without approximation, its major weakness is that it is computationally very 
intensive. In fact, it takes several months of computer time on several SPARC-10 
workstations working simultaneously to generate a single graph of the probability of 
bit error for this system. At the end of the thesis we conclude that, although our 
model is mathematically correct, a discrete time approach to the problem is probably 
a more efficient investigative tool in analyzing the performance of dense WDM fiber 
optic networks. 



Photodetector Output Current = 7J|s(<)| 2 


x = n£ \s(t)\ 2 dt n = J\{t)dt 

Figure 1: Channel 0 00K receiver structure. 
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II. ANALYSIS 


An OOK, dense WDM system utilizing single-cavity Fabry-Perot (FP) optical 
filters as demultiplexers consists of M + 1 transmitters (fixed wavelength lasers) con¬ 
nected over a fiber link to M + 1 receivers. Each receiver contains a frequency selective 
Fabry-Perot filter to demultiplex one of the M + 1 channels. In our derivation of the 
closed form expression for the probability of bit error for this system, M is the even 
number of adjacent channels symmetrically placed in frequency around the Channel 
of Interest (Channel 0) whose carrier is transmitted on an arbitrary wavelength A a . 
The Channel 0 receiver is shown in Fig. 1. The dense WDM optical signal is received 
by the Fabry-Perot filter tuned to A a , which allows the Channel 0 data signal to pass 
and rejects signals in the adjacent channels. A photodetector then converts the fil¬ 
tered optical signal to a current. The photodetector has an arbitrary responsivity TZ 
(A/W). The current is then amplified by a low-noise amplifier which adds a postde¬ 
tection thermal noise n(t ) with two-sided current spectral density N 0 (A 2 /Hz). Both 
the signal and noise currents are now passed to an integrator (low pass filter). The 
output of the integrator is the decision variable Y, which is compared to a threshold 
Vt to determine whether an “0” or “1” was sent. 

A. THE DENSE WDM SIGNAL AT THE INPUT OF THE FABRY- 
PEROT FILTER: r(t) 

Again we note that Channel 0 is our channel of interest. We can write the 
expression for the Channel 0 data signal in the complex baseband as 

b 0 (t)= E bo, iPT (t-iT) (1) 

t=—Lo 

where T is the data bit period (s), and 6 0 , t is the bit in Channel 0 during the time 
period [zT, (z + 1)T]. Note that b 0 ,i e {0,1). L 0 is a positive integer which represents 
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the number of bits in Channel 0 that are trailing the detected bit or bit of interest 
bo,o- Pr{t ) is the rectangular pulse function defined as 


Pr(t) = 


1 , 

0 , 


0 < t < T 
otherwise 


( 2 ) 


Channel k is any arbitrary adjacent channel. There are M adjacent channels placed 
symmetrically in frequency around Channel 0. M is an even integer. The indices for 
k are as follows k = —M/2, • • •, —1,1, • • •, +M/2. We may now write the complex 
baseband expression for the k th channel data signal 

h(t)= £ b k ^ kt p T (t - £T) (3) 

e=-L 

where co k is the radian frequency spacing between Channel 0 and Channel k and 
uj k — ~u_ k . We have already defined the pulse function pr(t) in Eq. (2) above. L is 
a positive integer which represents the number of bits in Channel k which trail the 
bit in Channel k that is the 0 th bit, b k ,o. We now note that b k / is the bit in Channel 
k during the time period [IT,(l + 1)T] and that b ki £ where <j> k is Channel 

&’s phase offset from Channel 0. (f> k is assumed to be a uniformly distributed random 
variable between 0 and 2x (<p k ~ t/[0,27r]). We can now use Eqs. (1) and (3) to write 
the complex baseband equivalent dense WDM received optical signal which appears 
at the input of the Channel 0 Fabry-Perot filter (see Fig. 1) 

+M/2 

r(t) — VP b 0 (t) + £ VP b k (t) (4) 

k=-M/ 2 


and 

/ 0 +M /2 0 

r(t) = # E KiPrtt -iT)+ Y. E hje^PT(t - IT) 

V, '=- io k=-M/2 e =~ L 

k^O 

where P is the received optical power. 


( 5 ) 
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B. THE OPTICAL OUTPUT OF THE FABRY-PEROT FILTER: 

s(t) = r 0 (t) 

To arrive at the optical output of the Fabry-Perot filter s(t), we must convolve 
the input signal r(t ) with the impulse response h(t ) of the filter 

s(t) = r 0 (t) = r(t) * h(t ) = h{t) * r(t) (6) 


First, however, we need the expression for h(t). The Fabry-Perot filter is a causal, 
linear, time-invariant (LTI) system. It can be shown [1] that 


h{t) = (1 -p)f^p 9 8 ft - | 


( 7 ) 


5=0 

where 6(t) is the Dirac Delta Function defined in two parts 

m = 

and 


0, 

0, 


t < 0 
t > 0 


(8a) 


/° + 6{t)dt = 1 (8b) 

Jo- 

Looking at Eq. (7), we note that p is the power reflectivity of the single-cavity Fabry- 
Perot filter and f3 is the filter’s free spectral range (Hz). Performing the convolution 
operation yields the following (see Appendix A for details) 


s(t) = r 0 (t) = y/P (1 — p) 


K 


P 3f) o,oPT (t ~ 


-9=0 


0, 


— loo / 

+ E E p 9 b o,iPr (t 

i=-Lo o=0 V 


0 


SB(t) 


iT 


+M/2 


5/S/(0 
0 CX) 


+ £ EE P 9 h,ee jwk{t ~ {9m PT U-^-lT' 

k=-M,2 V P } 

0 

>- v - 


( 9 ) 
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where ) is the desired signal (signal of interest), sisi(t ) is the intersymbol inter¬ 
ference signal, and SACi{t) is the adjacent channel interference signal. Then 

s(t) = r 0 (t ) = K (sg(t) + s/si(t) + SAci(t)) (10) 

We are now interested in writing an expression for s(t) during the detection 
interval 0 < t < T. We do this by sketching s#(£) from Eq. (9) and then writing 
another expression which accounts for how the rectangular pulses behave during the 
detection interval. We then sketch the pulse functions for the first three values of i 
for sjsi(t) in Eq. (9). We note how they appear during the detection interval, and 
recognizing a pattern, we write another expression for s [SI (t ) for 0 < t < T. A similar 
process is followed to arrive at SACi(t) for 0 < t < T. We now present $#(£), sisi(t), 
and SAci(t) for 0 < t < T (see Appendix A for details) 


7-1 

s B(t) = E P 9 bp,oPT 

3=0 


rji 9 


N V 


for 0 < t < T (11) 


where q = T/(l//3) = 0T is an integer. Recall that T is the data bit period (s), and 
P is the single-cavity Fabry-Perot filter’s free spectral range (Hz). 


f '"A"' ,l ( PT(t) 

SlSl[t) = E E P b 0,iPT - , n , -x 

i=-L 0 |_3=—(l+t)9+l \9 + (l+ *)? 


-(*'-1)9-1 

+ E P 9 t>0,iPT 


l±2S_ 

0 


9±iq_ 

P 


, for 0 < t < T (12) 
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/ 


+Mf 2 

SACl(t) = E 

k=-M/2 
k± 0 


/ 


lg=0 


V 


t-i 


r - 


P, 


+ E 


l=-L 


(-*»)-1 

E P 9 bk,ie j “ k{t - (9/l3)) PT 

g=-(l+e)q+l 


u /?; 

mo 




-(/-1) 9 -1 

+ 2 x 

g=-£q 


( 


Pt 


V 


£+i£ 

, p > 


t - 


' g + tq ' 

, p , 


\ 


) 


, for 0 < < < T (13) 


) 


Recall Eqs. (9) and (10) where 

s (0 = ^( 5 b(0 + S isi(t) + s ACl(t)) 


and 

# = ^( 1-/0 

We can see that s(t) in the interval 0 < f < T can now be calculated by substituting 
Eqs. (11), (12), and (13) into their appropriate positions in Eq. (10). 

C. COMPUTATION OF X = 11 [ T \s(t)\ 2 dt, THE SIGNAL COMPONENT 

Jo 

OF THE DECISION VARIABLE 

Looking at Fig. 1, we see that the filtered optical output of the Fabry-Perot 
filter s(t) is passed to a photodetector. The output current of the photodetector 
is lZ\s(t)\ 2 where 1Z is the responsivity of the photodetector (A/W). Thus, we can 
see that the deterministic signal component of the decision variable at the integrator 
output is 

x = n [ T | s{t)\ 2 dt (14) 

Jo 
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We have already derived s B (t), s IS i(t ), and s AC i(t ) for 0 < t < T. Directly 
substituting these three signals into Eq. (10) yields s(t ) for 0 < t < T. Now our true 
task is to compute |s(f)| 2 and / 0 T \s(t)\ 2 dt as 1Z is an arbitrary constant. 

s(t) = VP (1 — p)[s B (t) + sisi(t) + s,icj(0] (15) 

v -v-' s -v-- 

K s f (t) 

Then 

\s(t)\ 2 = K 2 s’(t)s'(ty (16) 

and dropping the ( t ) notation for convenience in the three terms of s'(t) and noting 
that complex conjugation is a linear operator 


s'(t)* — S* B + s* ISI + 


’ACI 


So 


S ( t)s ( t ) — (sb + Sisi + S A Cl){s B + S*ISI + S *ACI ) 


Multiplication yields 


s (t)s'(t)* = s b s* b + s B s* IS[ + s b s* aci 

+SISIS* B + S IS lS* ISI + SISIS* ACI 
+SAClS B + SAClS* ISI + SACIS*ACI 


Rearranging and simplifying yields 


(17) 


(18) 


(19) 
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\s{t)\ 2 = K 2 s'(t)s'(ty = I< 2 x 


sbs*b + srsis*isi + s acis\ci 


|sfl| 2 k/S/P \sACl \ 2 

+ SgS* IS j + S[SI S *B 


2i?e[sB"S/ 5/ ] = 2 Re[sisis* B ] 


+ s bs\ci + s acis*b 


2 Re[s* B s AC i] = 2 Re[s B s* ACI 


+ s ISI s *ACI + S ACI S *[SI 


2Re[s isi s aci. — 2 Re[sjsi s aci\ (20) 


We will integrate each of these terms over 0 to T to compute Jq \s(t)\ 2 dt. 

Substituting the appropriate definitions into the six terms of Eq. (20) yields 
21 terms (“clusters of summations”) which are integrated from 0 to T. Also note 
that we had to create several “gating” functions to turn off integrals when the pulse 
function products in some of the 21 terms fail to overlap for certain summation 
indices. Although intuitively obvious, see Appendix A [below Eqs. (104) and (112)] 
for definitions of the max and min functions which appear in some of the forthcoming 
equations. 

Recall that K = \/P (1 — p ) and K 2 = P(1 — p) 2 . Then we have 

[ T \s(t)\ 2 dt = K 2 [ T s'(t)s'(tydt = I< 2 [ T \s\t)\ 2 dt (21) 

which has the following terms [see Eq. (20)] 


K 2 [ T \s B \ 2 dt = P( 1 - p)X, a E E/ + ”(r-( 

o=0 m=0 V V 


max(< 7 , m)' 


( 22 ) 
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(23) 


where 


in which 


K 2 f \ s isi\ 2 dt = Msi + Bisi + Cisi 

Jo 

-1 (-»g)-l -1 (-rq)-l 

A-ISI = J°(l - pf Y Y Y Y Vo 

2=—Lo g=-(l+x)g+l r=-Lo m=—(l+r)g+l 

Y>„ = / + ”4„A,, ( minb + (l +, ^ m + (l +r)<? ) ) 


-1 (-w)-l -1 -(r-l)g-l 

B,s, = 2P(\-pY E EE E ^ m bo,ibo. r G l (g,i,m,r) 

2 — — Lq g~ — (1+ 1 )g-f-1 r— — Lo m— — rq 

in which 

( -^[(^ + (1 +*)?) - (m + rq)], for m + rq < g + (1 + i)q 


Gi{g,i,m,r) = 


otherwise 


-1 -(i-l)q-l -1 -(r-l)g-l 

c,s, = p(i-p)\e E E E / +m v.v t 

i~ — Lo g——iq r= — Lo m~-~rq 


max(g + iq,m rq) 

P 


where 


K 2 f \sACi\ 2 dt — Aaci + Baci + Caci + Daci + Eaci + Faci ( 24 ) 

J 0 


+M/2 q— 1 +M/2 q— 1 

a aci = p(i-pY E E E E^x^ 

k=-M/2 9 ~° n--M/2 m=0 


in which 


^,0^n,0 •/ \ x ‘' ‘ 

’ j(u k -cu n ) 

\ e j[(uk-un)T-Uk( 9 /P)+Un(m/ 0 )] _ „j[(w k -ui„)(l/f3)max(g,m)-w k (g/0)+uj n (m/f3)] 


b k ,obt fi e^/P){m fl)] 7 1 _ I max(< 7 , m) , 


, for & ^ 


for k = n 
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+ M /2 -1 (-tq)-l +M/ 2 -x (-rq)-l 

b A c, = p(i-/>Y E E E E E E r*?. 


k——M /2 

k^O 


e—Lg—(l + e)q+l n= _ M /2 
n^O 


v — — L m = — ( l - f - r ) g - j-l 


in which 


J[-^fc(3//3)+w n (m//3)] 
;(^ -^n) 


^ e i[(^-‘*'n)(l//3)min(g+(l+^)g,m+(14-r)9)] _ ^ ^ 


for k 7^ n 


h kA b lj [e 3[(uJkl(3)( - m 9)] ^ min(fif + (1 + m + (1 + r)<?) ^ , 


for k = n 


+M /2 -1 -(^- 1) 9 -1 +M /2 -1 -( r - l ) g-l 

c, c , = p(i -/.) 2 E E E E E E r*». 


k=—M/2 ra=—M/2 r =“'- 


in which 


, - u „ )T-w k (g/ 0 )+u n (m/P)] 


h,eK,r — ( -T | e 

’ bfWjfc - u> n ) L 


_ J _ e i [(^ fc - Wn )( l // 3 ) max ( ff +< 9 , OT + rg )- a / fc ( 3 // 3 )+ u / n ( m // 3 )]| ^ f Qr £ n 


h,eK,r e ^ u>k ^^ m 7 1 — — max(<jf + A?, m + r<?) , for A: = 

r* 


11 




+M/2 q— 1 +M/2 -i — (r — 1)7 — 1 

D ACI = P(i-p ) 2 E E E E E 


k=—M/2 9=0 n=-M /2 r =- L " = _ "' 


in which 


p 9+m bkob* _ \ e 3[(vk-Vn)T-u> k (g/P)+w n (m/P)] 

’ n,T [j(u k - u> n ) L 


J[(vk-u n )(l/P)Tnax(g t rn+Tq)-uj k (g/f3)+<jj n (m/l3)\ 


_L n 9+mv* L _ \ p -j[(uk-Urx)T-uj k (g/P)+u> n (m/P)] 

~r u k,O u n,r r •/ c 


_^-i[( w fc-^n)(l//?) max(g,m-frg)-w*.(s//3)+w n (m//?)] 


for k ^ n 


b kt0 bl r e j ^ Wk ^^ m - 9 ^ + e _J ^ a, ' t I /3 ^ m-s ^] /> 5+m x • • • 


T ~ -jj ma x(g, m + rq ) 


for A; = n 
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+M/2 q -1 +M/2 -X (-rg)-l 

Eaci = p(i-pY E E E E E *>> 

k=—M/2 5=0 n=-M/2 r=-Lm=-(l+r),+l 


in which 


-- e 

- w„ ) L 


,i[(wfc-w n )((m+(l+r)9)//?)-a/ fc (3/^)+w n (m//3)] 




+ l /+ ” l ^’° 6 n ’ r [ -jK--n )] I" 


> -j[(vk-Un)({m+(l+r)q)/P)-uJk(g/P)+Vn(m/P)] 


—e j[(^n//?)(m 5 )] j G 2 (^, m,r), 


for & 7^ n 




-[(m + (1 + r)q) - g]j G 2 {g, m, r), 


for k = n 


where 


G2(9,m,r) = 


1, for g < m + (1 + r)g 

0, otherwise 
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■j-M/2 — l (—lq)—l -\-Mf 2 —1 —(r — l)g—1 

F iC , = p(\- P f E E E E E E ve 

,. = -t//2 <=-!»—(!+<)«+! ' = -t ” = 


A:^0 


in which 


¥>6 = < 




71^0 


3 i[(wfc-w„)((3+(H-^)9)//?)-w*(fl//3)+w„(7n//3)] 


_^i[(“»:-^n)((m+rg)//3)-w k ( 3 //3)+a/ rl (m//3)] 


+^ +m ^/ n 


1 


-j{iOk - U> n )] 


0 -j[(uh-Un)({g+{l+t)q)/{3)-w k (g/l3)+uj n (m//3)] 


— e -j[(uk-Un){(m+rq)/P)-v k {g/0)+u n (m/0)] \ ^ m? 


for A: ^ 77, 




p 9+m x ... 


-p[{9 + (1 + £)q) -{m + r?)]^ G 3 (g, t, m, r), 


for k = 


n 


where 


G 3 (g,£,m,r ) = 


f 1, m + rq < g + (1 + £)q ) 


{ 0, otherwise 

K 2 f 2 Re[sBs}si]dt — K 2 [ 2 [sBSisi\dt = 

J 0 Jo 


mi-P) 2 E E ( sT 

5=0 i=—Lo m=—(l+t)?+l 


9-1 -1 —(*—1)9—1 

+ 2 C(i-rt ! E E E / + ”v«v, 

5r=0t = —Lq m — — iq 


T — — ma x(^ ? m + zg) 


(25) 


where 




-t[(m + (1 + i)q) - flr], g < m + (1 + i)q 


0, 


otherwise 
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K 2 f 2Re[s* B s A ci]dt = K 2 f 2Re[sBSACi]dt = Ab-aci + Bb-aci + Cb-aci 
Jo Jo 


where 


?-i +A //2 7-1 

Ab-aci = P{ 1 - pf 12 12 12 P"? 

5=0 k=—M/2 m=0 
k± 0 


in which 


V?7 


_ 0 b k u ^ Lj a, fc( 3 ’ - ( m /^)) _ e J w k(( max (5.™)-m)//3) 




+/+ ™ Vo6 - 0 __ 


0 -jw k (T-(m/P)) _ -jw k ((max(g,m)-m)/(3) 


7-1 +M/2 -1 (-<?)-! 

Bb-aci = ^(1 “ p) 2 H S 

3=0 k = — M/2 i =~ L m=-(l+e)q+l 
0 


in which 


¥>8 = P 9+m bo,obk,e— 


JUk 


0 iwjb((H-^)g//3) _ 7u; fc ((5-m)//3) 




+/> 3+m &o,o&M 


[-M] 


-jwfc((l+^)g/0) _ p ~3Uk((g-rn)lP) 


G 5 (g,t,m) 


where 


Gs(gJ,m) = 


l, g < "i + (l +% ] 

0, otherwise 


7-1 +M/2 -1 -(<- 1 ) 7-1 

Cb-aci = ^*(1 — p) 2 X] SZ 

5=0 k=-M/2 e =~ L m =~ lq 
k^O 


( 26 ) 
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in which 


'P’o — p 9+m bo,obk,e~ 

JUk 


+P 3+m bo, 0 b] 


k,t 


[~juk] 


,ju k (T-{m/p)) _ jw k ((max(g,m+eq)-m)/P) 


,-jv k (T-(m/p)) _ -jui k ((max(g,m+tq)-m)/P) 


K 2 f 2Re[s* ISI SAci]dt — K 2 f 2Re[sjsiSAci]dt 
Jo Jo 

— A IS I-ACI + Bjsi-ACI + C[si-AC I 
+Disi-ACI + Eisi-ACI + Fisi-ACI 


where 


in which 


<£>io 


-1 (—*«) —1 +M/2 q — 1 

AlSI-ACI = -P(l - pf Y Y Y Y V10 

i=-L 0 g=-(l + i)q + l k = _ M j 2 m=0 
k±0 


= P 3+m b 0t ib kfi — [ e ^t(i//5)( 5 +(i+.»-m//3] _ jl Ge ^ m) 

J^k L J 

+p3+™bo,ibl i0j Y— _ xj Ge(g,i,m) 


where 


G 6 (g,i,m) 


1, m 


0, otherwise 


< 9 + (1 + i)q 


-1 (-t'g)-l +M/ 2 -1 (-eq)-l 

Eisi-aci = p^-pYY Y Y Y Y Vn 

i=-L 0 j=-(l+t)g+l k=-M/2 * =-i m =-(l+*)?+l 

kq£0 


in which 


¥*11 — P 9+m bo,ibk l i 


3&k 


Ju k [(\!p)mm(g+(l+i)q,m+{\+t)q)-m/P] _ jw k [-m/P] 


+p 9+m b 0 ibt . _ L-jw fc [(l/^)min(3+(l+t') ? ,m+(l+^) 9 )-m//3] _ & - 3 w k [-m./P] 


( 27 ) 
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in which 


¥>12 — 

where 


in which 

¥*13 


in which 

¥>14 = 


-1 +M/2 -l —(£ — l)9—l 

ClSI-ACI = P( 1 - p) 2 E E E E E 'i’ 12 

i=-Lo 3 =-(l+i)g+l k——M/2 ^ m= ~^ 

k±Q 


p 9 + m b Q M~r 


ju k 


e i<^fc[(l//3)(£f+(l+*)9)—m//3] _ e joj h [tq/P] 


+P 9+m boM,i 


1 


L-JWfcj 


3 -ju; fc [(l//?)(^+(l4-09)- m /^] __ p -jvk[*<l/0\ 


G r (g,i,£,m) 


G 7 (g,iJ,m ) 


I" 1, m + < g + (1 + i)q 

[ 0, otherwise 


-1 —(*—1)9—1 +M/2 q -1 

Disi-aci — ^(i-/>) 2 E E E E ¥>13 

k=-M/2 
0 


= /9 s+TO 6ot^o~ — [e^^-W/?)] — P)raax.(g+iq,m)-mi p] 


JUk 


-LnP +m h A* _L_ \ p -jUk{T-(m/p)] _ jw*[(l//?) max(s+*9,m)-m//J] 


-1 -(*-1)9-1 +M/2 -1 (-£9)-! 

Eisi-aci — p( 1 - />) 2 E E E E E v>i 4 

t=-L 0 g=-iq k——M/2 m=-(l+^)g-j-l 

A:^0 


p9+ m b() t ib k ( — 






[-M] 


0 jw>:Kl//?)(l+£)9] _ [(l//3)(p+i9 )—ttt//?] 


-twfc[(l//9)(l+£)9] _ „-i"fc[(l//3)(9+*9)-™//3] 


G a (g,i,Z,m) 
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where 

1, g + iq < rn + (1 + t)q 

G s (g,iJ,m) = < 

0, otherwise 

-1 —(*'—1)9—1 +M/2 -l —(£—l)g—1 

Fisi-aci = P(l - p) 2 Y 2 2 £ il V15 

t=_Lo 3= ~ ! ' 3 k=-M/2 t= ~ L m= ~ eq 

k^O 

in which 

(vj 15 — p 9+m b 0 ibk t—— — e -J a '*[( 1 /^) max (fl+*9> m +^9)- m // J ] 

’* ’ JUk [ J 

+ / 0 9 +m 6 0 t - 6 ^_i_ e -jW[T-(m//?)] _ e ~jw k [(l/P)meix(g+iq,m+eq)—m//3] 

’* ’ [“JWfc] L J ^ 

We have now laid down all terms of 

pT pT y 

/ |s(t)| 2 dt = A" 2 / s'(t)s'{t)*dt = K 2 |s'(()| 2 dt 

•/O Jo Jo 

This will be used to compute the deterministic signal detection statistic X where 

X =n [ T \s(t)\ 2 dt 

Jo 

This will be used to compute probabilities of bit error for the dense WDM system. 

D. THE DECISION VARIABLE Y 

The decision variable Y appears at the output of the integrator (see Fig. 1). 
It consists of a signal component X, which was presented in Section C, and a noise 
component N 

Y = X + N (29) 

where 

T 

X =71 [ |s(i)| 2 dt (30) 

Jo 
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and 


N = [ T n(t)dt (31) 

Jo 

We note that n(t) is a zero mean Gaussian random process with two-sided current 
spectral density N 0 (A 2 /Hz), and N is a zero mean Gaussian random variable with 
variance N 0 T. Y is compared to a threshold Vj to determine whether a “0” or “1” 
bit is presented at the output. 


E. PROBABILITY OF BIT ERROR FOR THE DENSE WDM 
SYSTEM 


For a detection threshold Vj and an ACI/ISI bit pattern ip p = {&*,*, & 0 ,;} for 
k = -M/2, • • •, +M/2, k ^ 0, l = -A, • • •, 0, and i = —L 0 , • * •, —1; the conditional 
probability of bit error for the dense WDM system utilizing the decision variable Y 
described by Eqs. (28)-(30) is given by [2, 4] 




i ( Vt-MM ' 

2 V l vW , 


(32) 


where 


Q(x) = -L. jf" t-^dy (33) 

Xq is the value of X given by Eq. (14) or (30) with the value of / 0 r \s(t)\ 2 dt obtained 
by summing Eqs. (22)-(27) when 6 0i0 = 0. Xi is obtained in the same manner with 
bofi = 1. The average probability of bit error Pb is given by taking the expected value 
of ■Pe(V’p) given in Eq. (32) over all possible bit patterns xp p [2]. Let us define an 
ACI/ISI bit pattern set V’ = {V’p}; P = 1>‘” ,NPAT. NPAT is the total possible 
number of bit patterns in the set ?/> = {V’p}- Then 


a = {aw,)} (34) 

{4>p} 

If we count up the number of independent bits in ?/> p , there are M(L + 1) + L 0 bits. 
Using the assumption that a 0 or 1 is equally probable, we see that there are two 
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possible ways to fill each bit position. Then, by the generalized multiplication rule, if 
we let NPAT be the total number of possible bit patterns in the set ip = {ip p } 

NPAT = (2),(2) 2 ■ (2)„,i +1)+It = (35) 

If we assume that all bit patterns in ip = {ip p } are equiprobable (an excellent assump¬ 
tion), then by the law of total probability 

Pb = 2 M ( L + 1 )+ L ° jPe ^ 1 '* + 2 M ( i+1 )+ Lo ' Pe ^ 2 ^ ^ 2 m ( l + 1 )+^o P ^npat) 

and 

i NPAT 

P > = 2 M(L+i)+Lq £ P.(M 

p=l 
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III. NUMERICAL RESULTS 


We attempt to generate four probability of bit error graphs. We limit our 
expression for X given by Eq. (14) or (29) and the sum of Eqs. (22)-(27) by setting 

Lq = 1 L — 0 (f>k — 0 

and 

2nkl 

“k = —jT~ 

where k is an integer and I is the normalized channel spacing integer (/ > 0). When 
we do this, our 21 terms (“clusters of summations”) reduce to ten terms (“clusters 
of summations”). In the interests of brevity, we refer the reader to Appendix B for 
the details on how we came up with the ten programmable terms of X for this case. 
Referring to Appendix B, the probability of bit error equations we use to generate 
these graphs mirror the equations presented in Chapter II, Section E. Referring to 
Appendix C, one finds the programming strategy for the ten terms and the final 
graphs, the actual programs for each of the ten terms, and the final four graphing 
programs which produce one graph each. Each of the four graphs plot probability of 
bit error versus signal-to-noise ratio Z — 1Z P\JTjNQ in dB. Each graph has five traces. 
One trace shows either Single Channel (SC) operation without Fabry-Perot (FP) 
filtering or SC operation with FP filtering and without ISI and ACI for comparison 
with the other four dense WDM traces. The other four traces show probability of 
bit error versus Z (dB) for four selected values of the normalized channel spacing 
integer /, or equivalently the number of adjacent channels M. See Appendix C for 
the relationship between / and M. Each of the four graphs is plotted for a different 
value /3T, the free-spectral range-bit period product. We also note that the four 
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graphs are generated with 


p = Single Cavity Fabry — Perot Filter Power Reflectivity = 0.99 


We now present the values of (3T the four values of M, and the corresponding 
values of I used for each of the four graphs. Each separate value of M, or equivalently 
/, produces one of the four dense WDM traces for use in comparison with the fifth 
trace showing single channel operation. 

For Fig. 2: (3T = 500 and 


M 


For Fig. 3: /3T = 1000 and 


M = 


For Fig. 4: 0T = 1500 and 


M = 


For Fig. 5: /3T = 2000 and 


M = 


124 ,98^ ^ ^24 

L 1=4 1=5 1=8 1= 20 J 


198, 164 110 ^48, 

L 1=5 1=6 1=9 1=20 


212 164 124 74, 

L 1=7 1=9 1= 12 1 =20 


248 220 164 ^98, 

L 1=8 1=9 1=12 1=20 


Looking at Figs. 2-5, we see that only Fig. 2 (/?T = 500) is a complete graph. 
For @T = 500 we were able to compute all ten terms of X (see Appendix B) for all 
four given values of M , or equivalently /, and completely compute and simulate all ISI 
and AC I effects of our complete model given the constraints given at the beginning of 
this chapter (L 0 = 1, L = 0, <j> k = 0, u k = 2ivkl/T). For Figs. 3-5 (f3T = 1000, 1500, 
and 2000, respectively), we were unable to arrive at a solution for Eq. (24) (see Term 
#5 in Appendix B), which is the very important K 2 / 0 T \s A ci\ 2 dt term. It is very com¬ 
putationally intensive as it involves computing a quadruple summation multiple times 
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Figure 2: 


13 14 15 16 17 

Z(dB) 

of bit error vs. Z (db) for f3T = 500. 
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Figure 3: Probability of bit error vs. Z (db) for (IT = 1000. Only I = 20 trace 
accurate (within 1/16 dB). 
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Figure 4: Probability of bit error vs Z (db) for (3T = 1500. Only I = 20 trace 
accurate (within 1/16 dB). 










for the various adjacent channel bit “loadings” (see Appendix B for the details), 
as well as having to repeat the same very lengthy calculations for each of the four 
given values of M, or equivalently I. However, we will still be able to draw some 
conclusions for the I — 20 traces of these graphs as the channel frequency spacing 
between adjacent channels given by 

= S = ? = < 37 > 

where Rb is the channel data rate or channel bit rate. Note that A/ is proportional 
to I. At these larger frequency spacings, ACI effects become quite negligible. Equa¬ 
tion (37) explains why we call 1 the normalized channel spacing integer as adjacent 
channels are separated from each other at an integer multiple I of the bit rate Rb- 
Let us assume a system design of a filter free-spectral-range f3 = 3800 x 10 9 Hz = 
3.8 GHz. Then we can say that for f3T = 500 our bit rate Rb = 7.6 Gb/sec, for 
f3T = 1000 our bit rate is Rb = 3.8 Gb/sec, for /3T = 1500 our bit rate is Rb ~ 2.5 
Gb/sec, and for fiT = 2000 our bit rate is R b — 1.9 Gb/sec. It is noteworthy that, in 
the time domain model, calculating probability of bit error as presented in Chapter 
II and Appendix A is easiest for the highest channel data rate Rb, or equivalently the 
smallest bit period T. 

We define power penalty as the increase in signal-to-noise ratio required for the 
dense WDM system to achieve a 10 -15 bit error probability over the Single Channel 
(SC) system which achieves this goal operating at Z = 12 dB. We can now make 
some power penalty statements. 

Looking at Fig. 2 (/ 3T = 500) we see that the power penalty increases as I 
decreases, or equivalently the number of adjacent channels M increases. Looking at 
Eq. (37), we see that as I decreases, the interchannel frequency spacing A/ decreases. 
Thus, it is logical that ACI effects would increase. In short, the channels become much 
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more tightly “packed”. Looking at the results for Term #5 in Appendix C when 
(3T = 500 we see that the magnitude of the contribution for I = 20 is only about two 
to four percent of the contribution for / = 4. Returning to Fig. 2 (j3T — 500), we 
see that the power penalty for I = 20 and M = 24 is approximately 2 1/8 dB; for 
1 — 8 and M = 60, the penalty is approximately 2 3/8 dB; for I = 5 and M — 98, the 
penalty is approximately 2 7/8 dB; and finally for / = 4 and M = 124 the penalty is 
approximately 3 1/2 dB. It is in this case of / = 4 and M = 124 that we see the true 
beauty of dense WDM. If we assume /? = 3.8 G Hz, then the channel bit rate R b is 
7.6 Gb/sec, and the aggregate system bit rate R b is 124 x 7.6 Gb/s = 942 Gb/sec, 
which is quite close to 1 Tb/sec (1 Tb/sec = 10 12 b/sec). Thus, for slightly more than 
doubling the signal-to-noise ratio Z = 7Z P^tJNq, we have increased the aggregate 
bit rate of one single fiber close to a terabit per second. Also, the dense WDM 
system utilizes a noncoherent receiver and, thus, is much cheaper to implement than a 
coherent optical FDMA network with expensive synchronization circuitry. Looking at 
Z , however, we see that increasing signal-to-noise ratio is not that easy to do. Upper 
limits have already been approached in making low-noise current amplifiers; thus, the 
true technical problem to be solved in the dense WDM system utilizing single-cavity 
Fabry-Perot filters is to place an optical amplifier before the FP filter (see Fig. 1) 
to boost the received optical power P. With this problem solved, economical dense 
WDM systems utilizing single-cavity FP filters with aggregate bit rates with many 
terabits per second are possible. 

Now, in Figs. 3-5 we have already seen that at I = 20, the effects of ACI 
become quite negligible as the channels are spaced farther apart in frequency. Thus, 
we can conclude that the I = 20 probability of bit error traces are accurate for 
(3T = 1000, 1500, and 2000. The rest of the traces for / < 20 are not. A manuscript in 
preparation by Tri T. Ha entitled, “A Discrete Time Approach to Sensitivity Analysis 
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of Direct Detection Optical FDMA Networks with 00K Modulation,” considers the 
same problem as this thesis, but uses a discrete time approach. The probability of 
bit error graphs in this future paper are generated with the same parameters used for 
this thesis. For the discrete time approach, the N = (3T = 1000 and 1500 traces for 
/ = 20 match the I = 20 traces for this thesis to within less than 1/16 of a dB, and 
the N = (3T = 2000 trace for I = 20 matches the one in this thesis to within 1/5 of 
a dB. 

Finally, using the case /3T = 500 we feel that our time domain model is very 
accurate as Fig. 2 agrees with the corresponding discrete time graph for N = 500 to 
within 1/4 of a dB for I = 4 and M = 124, agrees to within 1/8 of a dB for I = 5 
and M = 98, and matches almost exactly for I = 8 and M = 60, and I = 20 and 
M = 24. These differences are well within the realistic bounds of numerical error as 
our time domain approach requires many billions of calculations. 
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IV. CONCLUSION 


We have presented a complete model and equation for the probability of bit 
error of an 00K, dense WDM system employing single-cavity Fabry-Perot filters 
as channel demultiplexers. Our expression completely models Inter-Symbol Interfer¬ 
ence (ISI), Adjacent Channel Interference (ACI), and phase offset between Channel 
0 and Channel k for all M of the adjacent channels. No simplifying mathematical 
assumptions have been made to arrive at the final answer contained in Eq. (14), Eq. 
(29), and the sum of Eqs. (22)-(27). Others [1, 3] have worked on this problem; 
however, the solution is arrived at in the frequency domain and often involves sim¬ 
plifying mathematical assumptions [1, Eqs. (15), (36), (39) and 3, Eq. (5)]. The 
closed form expression we have presented for probability of bit error for the dense 
WDM system, although mathematically rigorous and an extremely complete model, 
has one significant drawback—it is extremely computationally intensive. The single 
complete graph for probability of bit error when /3T = 500 took three months to 
compute using multiple SPARC 10 workstations working 24 hours per day. Even for 
our limited case of L 0 = 1, L = 0, u k — 2TrkI/T, and (f) k = 0, we never arrived at 
solutions for the other three values of f3T = 1000, 1500, and 2000. The solution is to 
use a discrete time approach. Again, we are confident of the accuracy of our model 
as the same problem was attacked independently via discrete time approach, and our 
graph for /3T = 500 matched the discrete case of N = f3T = 500 to within the bounds 
of numerical error. Finally, we again mention that to practically implement dense 
WDM systems with aggregate bit rates of many terabits per second, large increases 
(doubling, tripling,...) in the SNR of the dense WDM system over that required for a 
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single channel system are necessary. To accomplish this an optical amplifier will have 
to be placed in front of the Fabry-Perot filter to greatly increase the filter’s received 
optical power. 
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APPENDIX A 

DERIVATION OF ft |s(<)| 2 <ft FOR USE IN X = 11 |s(t)| 2 A 

Before we begin the derivation of / 0 T \s(t)\ 2 dt, we must first define the variables, 
equations, and terms involved: 

Y — The decision variable appearing at the output of the integrator (see Fig. 1) 
according to 

Y = X + N 

X — The deterministic signal portion of the decision variable 

x = n [ T \s(t)\ 2 dt 

Jo 

7£|s(t)| 2 — Output current of the photodetector 

TZ — Responsivity of the photodetector (A/W) 

n(t) — Postdetection thermal noise with two-sided current spectral density N 0 (A 2 /Hz) 

N — Random signal component of the decision variable Y. N is a zero-mean, Gaus¬ 
sian random variable with variance N 0 T 

N — [ n(t)dt 

Jo 

p — Power reflectivity of the single-cavity Fabry-Perot filter 

j3 — Free spectral range of the single-cavity Fabry-Perot filter (Hz) 

T — Data bit period (s) 

P — Received optical power (W) 
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r 


w(() “ { J; !therwise } (38) 

Desired Channel/Channel of Interest — Channel 0 

Adjacent Channel — Channel k : k = —M/2, * • •, — 1, 1, • • •, +M/2 where M is an 
even integer 

6 0 ,o — The detected bit of interest in the interval 0 < t < T 

b 0} { — Bit in Channel 0 during the z th time interval [zT, (i + 1 )T] where 

bo,i C {0,1} 

bk,t — Bit in Channel k during the I th time interval [£T, (£ + 1 )T] where 

b k ,e G {0,e^*} 

and 

j = y/^i 

<f>k — Phase offset between Channel 0 and Channel k. (pk is assumed to be a uniformly 
distributed random variable between [0,27r] 

4>k ~ U[ 0, 2tt] 

u>k — Radian frequency spacing between Channel 0 and Channel k. Channels are 
symmetric around Channel 0, i.e., 

U k = ~LU_ k 

h(t) — Impulse response of the Channel 0 single cavity Fabry-Perot filter 

m =(!-/>) (t - (39) 
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*(*)- 


The Dirac Delta Function defined in two parts 

' 0, t < 0 


m = 


0 , 


t > 0 


(40a) 


and 


f 6(t)dt — 1 (40b) 

J o- 

bo(t) — Complex baseband equivalent data signal in Channel 0 

o 

(41) 


L 0 — 


h(t) 


L 


r(t) 


b 0 {t) = b o,iPT(t ~ iT) 

i=—Lo 


An integer greater than zero representing the number of bits in Channel 0 that 
are trailing the detected bit &o,o- 

Complex baseband equivalent data signal in Channel k 

b k (t)= £ bk'te’ Wkt p T (t — IT) (42) 

l=-L 

An integer greater than zero representing the number of bits in Channel k which 
trail the 0 th bit in Channel k , bk,o- 


Received complex baseband signal at input of the Channel 0 Fabry-Perot filter 

(43) 


+M/2 

r(t) = VPbo(t)+ v(P6fc(t) 


fc=-Af/2 
k^O 


We now begin the derivation of fjf \s(t)\ 2 dt which allows the computation of 

x =n [ T \s(t)\ 2 dt 

Jo 

where s(t) is the output of the Fabry-Perot filter with an impulse response of h(t) 


r(t) 


m 


s(t) 
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As stated before 


M<) = a -/>) (i-|) 

-f M/2 _ 

r(t) = y/Pb Q (t)+ E y/Ph(t) 

k— — M/2 
k^O 


Substituting Eqs. (41) and (42) yields 

( 0 +M /2 o N 

E KMt-iT)+ E E h^prit-iT) 

l= ~ L ° k=-M /2 f=-L ' 

A: 9^0 

We are interested in the detected bit 6 0 ,o during the detection interval 0 < 
we need only evaluate s(t) = r 0 (i) in the interval 0 < t < T 


/ oo r oo 

h{t — T)r(r)dT — / r{t — T)h(r)di 
-OO J — OO 


s(0 = r 0 (t) = r(£) * /i(0 

= h(t ) * r(i) 

/ +M/2 

= MO* (a/PMO + v/P E MO 

k=-M/2 
k^O 

As convolution is a linear operator and distributes over addition 

/ + M/2 

s(t) = r 0 (t) = l\/p(h(t) * b 0 (t )) + \fp E MO * MO 

k=-M/2 
k± 0 

/ +M/ 2 

MO = MO = MO * MO + E MO* MO 

k=—M/2 
k^O 


I (44) 

t < T, so 

(45) 

(46) 

(47) 

(48) 
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Now looking at the term within the parentheses 

+ M/ 2 0 

h(t) * b 0 (t) + E MO * MO = MO * E bo,iPr(t - iT) 

k=-M/ 2 ,=_L ° 

k^O 

+M/2 0 

+ E MO* E bk^prit-lT) 

k=-M/2 l= ~ L 


Applying the distributive property of convolution 

+M/ 2 o 

h(t)*b 0 (t)+ E M0*M0 = E h(t) * (b 0ii p T {t - iT)) 

k=-M /2 *' =_Lo 

+M/2 0 / \ 


+M/2 0 / \ 

+ E E '■(<) * K<‘ n, PT(t - it) 


k=-M/2 

0 


Expressing convolution in terms of its integral definition 


+M/2 o -oo 

M 0 *M 0 + E M 0 *M 0 = Y, / M< “ r ) [ViPrlT - *T)] dr 

i=-L 0 • / -°° 


k=z-M/2 

fc#0 

+M/2 o 


u -oo 

+ E E / M* - r) b ktl e JWkr Pt {t - IT)] dr 

„ r — OO ^ J 


k=-M/2 


Substituting the expression for h(t — r) [Eq. (39)] into term Ax, above 

Ax = E / U “/OE^f* _t “ fl b o,iPr(T -iT)dr (52 

i=—Lo fl =0 V PJ 

Factoring (1 — p) and interchanging the order of integration and summation yields 

Ax = (1 — p) E E PH** [ 6 (* ~ T ~ ') M r - ^)dr (53 

i=—Ln 3=0 ■ / -°° V PJ 
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Using the identity: x(t)6(t — t 0 )dt = x(t 0 ), the integral evaluates to 


px(t — g j (d — iT). We then have the following for A\ 


£ roo 

M = £ / - T ) [ b o,iPT(r - iT)] dr 

r J —OO 


i=—Lo 


= (*-/>)£ J2p 9b ^PT\t 


i=-L 0 g=0 \ / 

Recalling the expression for 

+ M / 2 0 

B x = £ £ / h(t-r)\b k , e e^p T (T-iT)}dT 

/>__ r J-OO 1 J 

k=-M/2 u 


Substituting the expression for h(t — r) given by Eq. (39) into Eq. (55) above, we 


+M/2 0 


~ f - v poo w 

EE/ (i-o)E. 

fi r •/ —OO _ri 


k=-Mj 2 


p 9 d | t _ r _ | 6 M e^ T p r (r - UT) dr (56) 


Factoring 1 — p, interchanging the order of integration and summation, and factoring 
bk,t and p 9 from the integral 

+M/2 0 oo /-oo / n\ 

Bi = (l-p) £ £ Y<p 9 h,i / <5 t-r-- e J ^ T p T (T-^r)dr (57) 

k=-M/2 ‘=~ L °=° ^ V N 

fc# 0 


Applying the property x(t)S(t — t 0 )dt — x(t 0 ) 


e ju,fcr p T (r - iT)dr = e JUJk(t ~ (al(i)) p T f f - | - £t) 


Substituting Eq. (58) into Eq. (57) yields 

+ M /2 0 roo r , 

= £ £ / h(t — r) {b k ,ee 3UJkT p T {T - £T) dr 

/?_ r d — oo L J 

k=-M/2 L 
k^O 

-fM/2 0 oo / 


-fM/2 0 oo / 

(i - /=>) £ £ £ t 

l=—L 3=0 V 


k=-M/ 2 
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Substituting Eq. (54) and Eq. (59) into Eq. (51) and then substituting this result 
into Eq. (48) yields 


s(t) 


r 0 (t) = VP 


0 co 

(i-ri £ 

i=— Lq g —0 


Ai 



+M/2 0 oo (a \ 

+ {i-p) £ £ (t-l-ir) 


k=—M/2 
k± 0 


B 1 

We can now factor out (1 — p) from our expression above. Then 

0 oo / 


s(f) = r 0 (t) = VP(1 -/>) 


I] Xy & o,iPr -* r 

U=-L 0 fl=0 \ P ) 


+M/2 0 oo / n \ 

+ £ £ E/Swe^f-WOMp,. «-£-«■ 

A:^0 


(60) 


(61) 


Looking at Eq. (61). We separate the terms involving the detected bit of interest in 
Channel 0, 6o,o, the trailing bits in Channel 0, 6o,i i ^ 0, and the terms involving bits 
in the other channels, 6^, k / 0 


s(t) = r 0 (t) = y/P (1 — p) 


K 


p 31) o,opt yt - 

s B (t) 


+ £ £/V.w(i-f-ir) 

l = ~Ln 0=0 \ P / 


-L 0 g 

vr-“ 

s/S/W 

+M/2 0 co 


+M/2 0 oo (a ^ 

+ £ £ £ U - | - IT 


k=-M/2 
0 


SAC/W 


(62) 
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where s B (t) is the desired signal, s IS i(t) is the Intersymbol Interference (ISI) signal, 
and s AC i(t ) is the Adjacent Channel Interference (ACI) signal. Using this compact 
notation, we have 

s(t) = r 0 (t) = I< ( s B (l) + Sisi(t ) + s A ci(t )) (63) 

Sketching s B (t) for h 0 ,o = 1 



We can also sketch the intersymbol interference due to the trailing pulses in Channel 

0 


♦ 



ISI contributions to the pulse of interest 


Now let us assume in one pulse interval T that 


T 

j = q 
0 


( 64 ) 
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where 1/(3 divides evenly into the period T making q an integer quantity. We are only 
interested in the detection interval 0 < t < T, so we only need to evaluate s(t) = r 0 (f) 
in the interval 0 < t < T. Looking at the desired signal, which is the system response 
to the bit of interest 6 0 ,o 



We can see that, during the interval 0 < t < T, the response is a sum of scaled, 
shifted rectangular pulses. Therefore 

I 

g-1 

SB{t) = £ P 9f, 0,0PT 

3=0 

We now direct out efforts towards developing an expression for s/s;(f) during the 
interval 0 < t < T. We begin by examining the ISI term for i = —1 in Eq. (62), 
denoted by — 1 ISI. 

- 1 ISI = £ p 9 bo,-iPT (t - | + IT) , for - T < t < oo (67) 

3=0 V P ) 
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By examing the sketch of — IISI 



we see that 1 1 SI will also be a sum of shifted, scaled rectangular pulses for 0 < t < 
T. The -1/5/ expression will have two parts, the first of which is 

- 1 /5/(l) = P 9 K-\Pt (^) , for 0 < t < T (68) 


But with the q lh shift we no longer sum scaled pulses beginning at t = 0 and end¬ 
ing at 1//?, 2//?, • ■ but now we sum pulses which progressively scale and begin at 
]■/Pi 2//?> • • •, and end at t = T. This second part of the expression looks much like 
the expression for s B (t) during the interval 0 < t < T [Eq. (66)]. Thus, we may write 
the second part of — 1 IS I 




+ P 9 bo,-iPT 


g-q 


HTJj ■ 

/ 


for 0 < t < T (70) 


Examine the I SI term for i — — 2 in Eq. (62), —21 SI 


2ISI = J2 K -2 P 9 Pt U - i + 2T \, for 


We sketch —2/57 


-2T <t<oo (71) 


-2 T -2 T + - -T 


-2T+- 


-2T + 


(?- 1 ) 


■rl 


-2T + 


(2g-l) 


-r + 


-2T + 


(3? — 1) 


By analyzing this sketch, we can obtain 

- 2 IS 1(1) = 52 1 P’I’v.-iPT (£z%) - 


for 0 < t < T (72) 


-2 757(2) = 3 £Vft 0 ,-2Pr 


*-2<r\U l P 


for 0 < t < T (73) 


Note: These rectangular pulses shift into the interval 0 < t < T on the 2q th shift. 



Similar to — 1 ISI: 




we may write directly 


-3 IS I = J2 p 9 b 0 -3p T 

g=2q+l 


( mi ) 

\9-2q) 


4g—1 

+ £ p g K-ZPT 

9 - 3g 




for 0 < t < T (76) 


Then from Eqs. (70), (74), (76), the expressions for the -IISI, -2ISI, and -3ISI 


pulses, we are ready to use the pattern developed to write an expression for the ISI 



between 0 < t < T. Recalling 


sisi(t ) = E E P 9 bo,iPT [t---iT\ , 

i=-L 0 g=O \ P ) 

we may write directly from the pattern developed 


for - L 0 T <t<oo (77) 


wo = t T 

i=-L 0 S =-(l+i)g+l \9 + (1 + l)q 


-(i-l)q-l 

+ E p 9 Kpt 


g + iq 


g + iq 


for 0 < t < T (78) 


We now have expressions for sb( 0 and srsi{t), during the interval 0 < t < T. We 
have only the ACI term s^cii^) to calculate 

+A//2 0 oo / \ 

*Aci{t) = E EE P 9 h,it ilVk(t - (9lP)) p T < - § - £T) for - LT < t < oo 

*=-Af/2 5"- Lg= ° _ ^ _E 

fc^O * 


* Note: This portion of SAci(t ) is almost exactly like the ISI term with which we 
just worked [see Eq. (62)]. The only difference is the complex exponential factor. 
Looking at the 0 th term of the ACI in the k th channel 


OACI = p^h^-^pT ^ - |j , 


for 0 < t < oo (80) 


we sketch OACI 



We see, that in a manner similar to the sg(t) signal, the pulses add as before, but 
now the amplitude of the pulses is modulated by pSg3 u k(t-(ff/P)) versus just p 9 before. 
Then 


o ac i = Y J p 9 K^ Uk(t ~ (9l(i)) PT 

3=0 


for 0 < t < T (81) 


The —1 AC I will be similar to the —1 ISI term [see Eq. (70)]. However, now we not 
only have p 9 but also the complex exponential gated by the pulses in the summation. 

-1 AC I = J2^P 3 h,-ie^ k{t - {9/p)) p T 


+ £ p 9 b kt . l e j ^ t ~^p T 


T _ 

9 ~ 9 


g-q 

(3 


, for 0 < t < T 


The derivation of —2 AC I and —ZAC I will proceed in a similar manner as —21 SI 
and —31 SI [Eqs. (74) and (76)]. We see that two different type pulses gate the 
product of p 9 and a complex exponential term. These gated terms are then summed 
over the period 0 to T. We then use our expressions for 0 AC I, —1 AC I, ... to write 
an expression for SACi(t) during the interval 0 < t < T in a manner similar to the 
sisi(t) expression. We may now write the final expression for the ACI 


SAClif) = 


+M /2 

[ 9-1 

E 


k=~M/ 2 

U =0 

k^O 

\ 

-1 

(-M-i 

+ E 

E 

i——L 

3 =-(l+%+l 


g + (1 +I)q 


+ E x 

3=~*9 
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PT 


V 


T - 


, fi , 


t - 


' g + tq ' 


) 


for 0 < t < T 


(83) 


/ 


We have now derived s B (t), s ISI (t ), and s AC i(t) during the interval 0 < t < T. We 
now develop the expressions needed to compute |s(f)| 2 and Jq \s(t)\ 2 dt. Substituting 
the expressions for ss{t), sjsi(t), and s A ci(t ) f° r the interval 0 < t < T into Eq. (63) 
yields s(t) in the interval 0 < t < T 

s(t) — \[P (1 — p)[ 'SB(t) + sisi j t) + s AC i(t) ], for 0 < t < T (84) 

' K /S 7{t) ™” 

where s B (t), sjsi{t ), and s A ci(t ) are described by Eqs. (66), (78), and (83) respec¬ 
tively. Then 

\s(t)\ 2 = I< 2 s'{t)s'(ty (85) 

Dropping the (t) notation for convenience in the three terms of s'(t) and noting that 
complex conjugation is a linear operator 

s'(t)* — s B + s* ISI + s* ACI (86) 


So 

s'{t)s'(t)* — (SB + SJSI + S A Cl)( s *B + s *ISI + S A CI ) 
Multiplication yields 

s'(t)s'(t)* = s b s* b + s B s* ISI + s b s* aci 

+ s isi s *b + s isi s *isi + s isi s A ci 


+ S A CI S *B + s ACI s *ISI + S A CI S * A CI 


Rearranging and simplifying yields 

\s(t )\ 2 = K 2 s'{t)s'(ty = I< 2 x 


SbS*B + $ISI S *IS I + S A ClS* A CI 


bsl 2 h/s/l 2 


l s >4Cj| 2 


(87) 


( 88 ) 
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+ S B S ISI + SjsiS* B 

s -V-' 

2 Re[s B s* ISI ] = 2 Re[s IS is* B ] 

+ s B s *ACI + s ACI s *B 
- ---' 

2Re[s* B SAci] = 2 Re[sBS* ACI ] 

+ S ISI S ACI S ACI S ISI 

S -V-' 

2Re[s* ISI SACi] = 2 Re[sjsis* ACI ] (89) 


We will integrate each of these terms over 0 to T to compute / 0 T \s(t)\ 2 dt. To avoid 
confusion, we must realize that constant K 2 must be carried through to the final 
calculation of all quantities. Recalling Eq. (66) 

/ \ 


7-1 

SB p 9 b 0fi PT 

3=0 


T 


T-l 

\\ 


*_£ 


) 


We know 6 0 ,i E (0,1), p 9 is real, and pr(-) is also real. Then 


(90) 


Sb = S B 


(91) 


and 




l ( 

q-l 

Y P 3 KoPt 
g =o 


v 


\\ 2 


rp J 

W ~P. 


9\{ ft, 


J 


Factoring 6 0) o as it is a constant to the summation, we obtain 


(92) 


/ 


M a = % fi 


{ 


< 7-1 

Yp 9 Pt 

g=0 




\ 


' T _9_ 

W 


9 

t - - 

v ft. 


\ 


(93) 
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Performing the squaring operation yields 


/ 


I „ 2 l 2 

l s B - Oq,0 


/ 


g-1 q — l 

£ £ 

5=0 m =0 


V 


vv Pj 


t- 9 - 

K 0, 


Pt 


! m m’ 

T - 

vv /?, 


m 

* ~ J. 


( 94 ) 


Then 


A' 2 |sb | 2 = P(l-p)X.o 


q -1 9-1 

212 'EE/**w 

5=0 m=0 


f T-l 

H /?, 


v 0> 


Pt 


T-™ 
\\ Pj 


f m 1 


/j 


for 0 < t < T (95) 


Now we want to compute Jq K 2 \ss\ 2 dt [see Eq. (89)]. 

[ T K 2 \s B \ 2 dt = P(l-f) 2 tEE/**x- 

*^0 5=0 m =0 


f 


Pt 


' t _ £ 

VV /?, 


0 

t - 

K P< 


PT 


_ m ' 
T _ 

VV /*, 


f m 1 


dt 


(96) 


Looking at Eq. (96), we have interchanged the order of integration and summation. 
As the pulse functions are the only functions dependent on time in the expression, 
we look at the product of the two unit amplitude pulse functions in r] 
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| | => Product is 0 

Y/J => Product is 1 



We have arbitrarily assumed that m > g, but our logic would also work if we assumed 
g > m. The product of two pulses is 1 for ra//? < t < T and 0 otherwise. Thus, we 
see that the integral of the product of the two pulse functions in g will be a square 
area of height 1 x (T- largest of g/(3 or m//3). Thus 

£ K'utfdt = p( i - P )Xo ), for o < ( < T 

(97) 

where max(flqm) is defined by 

max(xi, x 2 ): Choose largest of x\ or x 2 , which are both positive. 

If xj = x 2 , then max(x l5 x 2 ) = x t = x 2 . 

We now need to compute / 0 T K 2 \s [SI \ 2 dt. We see from Eq. (78) that s IS i is real and 
not complex since f> 0 ,< € {0,1}. Thus 


$isi — S/s/ an d M/s/| 2 = s 2 IS i 


(98) 


M/s/1 2 = 


-i (—»?)-i 

E E 

i=-L 0 3=-(l + i)? + l 


f^bO'iPT 


a 


( \ 

\9 + (1 + i)q) 
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To perform the squaring we set up another version of a + b with the g-index changed 
to m and the Tindex changed to r to account for cross-product terms. We call this 
new representation c + d. Thus 


|s/ 5 /| 2 — (o T b)(c T d) — 

A B C V 

Performing the multiplications yields 

-l (-w)-i -i (-»•?)-1 

A = ac = Y! Y1 XI P 9+m b 0ti b 0 ' r 

i=—Lo g^=z—(l+i)q+l r=— Lq m=—(l+r)?+l 


X 


PT 


pT(t) 


g + { 1 +i)q) PT \™ + (1 + r)q / 


0T{t) 


( 100 ) 


B = ad = 


-l (-«<?)-! -i -(*—1)9-1 

E E EE p“+ m K<Kr 

i=—Lo g=— (l-f-z)q'-hl t— — Lq m——rq 


( mi \ 

PT U+(i+oJ PT 


l T l t _ f m + r< i\\\ 


0 


v 


T - 


m + rq 


) 


C = be 


-1 -(‘-1)9-1 -1 (-r?)-l 

E E E E x 

i—-L 0 g—~iq t=—Lq m = —(l+r)g+l 


/ T f t _ (9 + iq\\\ 


Pt 


P 


\ 


g_ + ui_ 

< 0 > 


Pt 


pm 


m + (1 + r)q / 


) 


V=bd = 


-i —(*-1)9-1 -1 —(’'—1)9—1 

E E E E / + "V,V x 

i=—L q 9 ——‘9 t——Lq m=—rq 


Pt 


Tit 


' g±jg \\\ 

, p )j 


g±jg 
, P t 


Pt 


T It- 


' m + rq ^ ^ ^ 

~P~ 


m + rq 

n~, 


( 101 ) 


( 102 ) 


(103) 
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Integrating over the interval of interest: I< 2 \s IS i\ 2 dt = ■ • • 

I< 2 J q \sisi\ 2 dt = K 2 £ Adt + £ Bdt £Cdt + £ Vdt (104) 

The equations for A, B, C, and V are Eqs. (100)-(103), respectively. We now analyze 
term A [Eq. (100)] by sketching the multiplication of the pulses 

PT f- ff g) ) 

/ \m + (l +r)qj 

□ =* Product is 0 
—_ i Y/A =>• Product is 1 


0 + (l + O? m + (l+r )<7 

P P 


am 


\9 + (1 + i)q) 

Again, we have arbitrarily assumed one pulse lasts longer than the other. Since the 
height of each scaled pulse is 1, we see that the area under the product is 1 times the 


minimum of 


q + (l+i)q m + (1 + r)qf 

—<— - or ---— 

0 t> 


We can apply the same logic as used before to arrive at the value of the term 


I < 2 / Adt 


Then 


rl -1 (-w)-l -1 (—rq)— 1 

l Adt = P(1 - „f £ E £ £ *>0 

t=-Lo g=-(l+i)q+l r-~L 0 m=~(l-fr) 9 +l 


(105) 
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where 


— P 9+m bo,ibo,r 


min (g + (1 + i)q , ra + (1 + r)q) 


where 


min (g + (1 + i)q, m + (1 + r)q) is defined by 

min(x!,X 2 ) : Choose smallest of x x or x 2 , which are both positive. 

If Xi = x 2 , then min (x l5 x 2 ) = X\ = x 2 . 

We now compute K 2 / 0 T Vdt. Looking at the expression for V [Eq. (103)], we sketch 
the product of the pulses, arbitrarily assuming m + rq > g + iq 


| | =>■ Product is 0 

Y/A => Product is 1 



g + iq m + rq 

0 0 


We see the product of pulses only exists from the maximum of 

g + iq m + rq 

___ or ___ 

Thus, the integral of the pfoduct of the pulses will be 

max (g + iq, m + rq) 


Then 


-l -(*—1)7-1 -1 —(*—1)7—1 


K*j vdt = £ E E E p’+ m KAr 

® i=-Lo g=—iq r=-Lo m—-rq 


max (<7 + iq, m + rq) 


( 106 ) 
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giving 

-1 (-rg)-l -1 —(>—1)9—1 

b = E E E E 

r—~Lo m=-(l-{-r)q-\-l 2 —— Lo 9=—iq 

VT ( ) PT ( r (‘~ (zzill 

(m + (l+r)qj P [y /i+iijj 


(108) 


If we look at the expression above Eq. (108), we see that 0 is identical to C [Eq. 
(102)] in every way except the order in which the summations appear. The order of 
the summations can be rearranged because the terms inside of the summations are 
completely separable in relation to the sets {m,r}, Thus, we can rearrange 

the order of the summations and conclude 


B = C (109) 

and that 

rj~i rp rj-i 

I< 2 [ Bdt + I< 2 [ Cdt = 2I< 2 [ Bdt = 2K 2 [ Cdt (110) 

Jo Jo Jo Jo 

We choose to integrate twice the value of B [Eq. (101)]. Interchanging the order of 
integration and summation yields 

T / -1 ( *9) 1 -1 -(r-l)9-l 

2K 2 J Bdt = 2/W £ E E E P 9+m Kib o,rX--- 

® \i=—Lo g=— (l+t') 9 +l r=—Lo m=—rq 
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(Ill) 



gm \ 

g + (l + i)q) 


Pt 


( 


\ 



m + rq\\ 
1 ~ 


m + rq 

~T~ J 


\ 


) 


dt 


Sketching the product of the two pulses, we see the product will exist for m + rq < 
g + (1 + i)q and will be zero otherwise. 




Product is 0 
Product is 1 


We create a “gating” function to account for this 



_jm_\ 

9 + ( 1 + i)q) 


Pt 



Gi(g,i,m,r) = 


Then finally 

2 K 2 [ T Bdt = 
Jo 


| ^[{a + (i+ *)?)- + 

I ' o, 


for m -f rq < g -f (1 + i)q 
otherwise 


( 112 ) 


I< 2 [ T Bdt + K 2 f T Cdt 
Jo Jo 

-1 (-«?)-1 -1 — ( r — 1) 7 -1 

2P(l-p)\Z E E E r) 

t=-Lo g~— (l+O^+l r —~ Lo m=—rq 

( 113 ) 
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Finally, combining the results of Eqs. (105), (106), and (113) we have 

rj-t rji rji » /ji 

I< 2 [ s ISI s* ISI dt = I< 2 f \s ISI \ 2 dt = K 2 [ Adt+ [ Bdtf C dt + [ Vdt 

Jo Jo Jo Jo Jo Jo 

(114) 

where 

T -1 (-tg)-l -1 (-rq)-l 

K>[ Adt = P(l - P f ■£ E E E Vo 

^ t=— Lq g~ — (l-|-z)g+l r~~Lo m=— (1-fr)g-f-l 

in which 

- rw., pto + a + 0t.m + (l+r),) j 

I < 2 f T Bdt + I< 2 f T Cdt = 2K 2 [ T Bdt 

Jo Jo Jo 

-1 (-w)-l -1 -(r-l)g-l 

= 2P(l-p) 2 E E E E v' + ”ViVrG.(s,',m,r) 

i=-L 0 g=-(l+i)q+l r=~L 0 m~-rq 


and 


T _1 -(t-l)g-l -1 -(r-l)g-l 

W Pdi = />(l-p) 2 X; E E E / + "«bAr X 

0 »— — r /I-— in r--- r.y. 7T7 — — TCI 


T- 


i——Lo g = — iq r = — Lq m~ — Tq 

'max ((7 + iq, m + rq) 
k 0 . 


Now we must compute l^^c/l 2 and A 2 / 0 T |s/ic 7 | 2 dA We see that we have a complex 
bk,e as bk,i £ {0,e J '^*}, as well as a complex exponential [see Eq. (83)]. 

We begin by noting that the complex conjugate operator has the following properties 

(uv)* = u*v* 

(u + v)* = U* + V* 


Taking the complex conjugate of Eq. (83) 


5 ACI 


+M/2 

E 

k——M/2 
kji 0 


/ 


V 


9-1 


■3=0 


T - - 


t - 


0, 
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Letting 


and 


-l 

+ E 

t=-L 


£ p 9 blee- ju,k{t - {9/0)) p T 

■ S =-(l+/)9+l 


( m) \ 

\9 + (1 + £)q) 


+ E P 9 K^~ iUk{t ' i9ll3)) PT 
g=-£q 


L (H 5 )) i 



c 


c 


f 


d 


d' 


PT ( m ) 

w U+( i+w 

T fm + (1 + r)q ) 


Pt 


Pt 


4-mn 

H’-tW 




T- 


m + rq 


Then, with a change of indices, we can express s* ACI as 


MCJ 


+Af/2 g — 1 

£ E a' 

,,. m=0 

n=—M/2 


n^O 


A* 


(115) 

(116) 

(117) 

(118) 

(119) 

( 120 ) 

( 121 ) 
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-f M /2 -1 (-rq)-l 

+ E £ £ c' 

n=—M/2 r=-L m =-(l+r) g +l 


+M/2 _i -(r-l)g-l 

+ E E E p m 6; )r e- j ^ (t - (m//3 ^d' 

r=-L m——Tq 

n=—M/2 
n^O 


( 122 ) 


-j-M /2 <7 — 1 

E E^^,oe^ (< - (5//J)) a 

k=-M/2 3=0 
k^O 


+M/2 _i Mg)-1 

+ E E E 

k=-M/ 2 e =~ L 9=-(l+t)q+l 


+M/2 _i — (^—l)g—1 

+ E E E 


k=-M/2 

k^O 


-L g=—£q 


(123) 


Then 


l^c/l 2 = *ac,s- 4C , = (A + B + C) (A- + B- + C-) 


(124) 


Now, distribution of multiplication over addition yields 


AA* + AB* + AC* + BA* + BB* + BC’ 


+CA* + CB* + CC* 


(125) 
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Rearranging yields 


= AA* + BB* + CC* + 2 - (AB* + BA* 


+2 [- (AC* + CA*)j + 2 (BC* + CB" 
.2 . .2 


Using the identity Re[Z] = (Z + Z *)/2 we have 


AA* + BB* + CC* + 2Re[AC* 


(126) 


+2 Re [AB*] + 2Re[BC* 


(127) 


Multiplying the terms above yields 

+ A//2 q — 1 +M/2 q — 1 

AA* = V Y) V V] p3+rn bk ob * noe j[(^-^)t-^U9/P)+^(m/P)] a& > 

k=-M/2 9=0 n——M/2 m= ° 


(128) 


+M/2 -i I-'*))-1 +M/2 -i (-rq)-l 

BB-= E E E E E E 7 

fa-M/2 <=-t»=-(l+0«+ 1 n =—M/2 '=-*■ m=-(l+r),+l 


(129) 


where 


- _ n 9+mt L* _j[(t*»jt-Wn)t-Wfc(a//3)+w„(m//3)] / 

7 — P u k,e Ll n,r c ^ 


+M/2 -l -(<-l) 9 -l +M/2 -1 —(r—l)?—1 

E E E E E E c 

fa-M/2 ‘=~ L ‘=-‘< .— M/2 r =“ L 


(130) 


where 


£ _ p9 +rn }) ke })* e j[( u k-^n)t-Uk(9/0)+Un(m/0)]^l 


To compute 2Re [AC*], we know Re[-] is a linear operator so it may be moved across 


the summations. After multiplication of terms A and C* we get 

+M/2 q- 1 +M/2 -1 -(r-l) 9 -l 

2Re [AC*] = 2 E E E E E « 

fa-M/2 >=° n=—M/2 1 ” = 


(131) 
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where 


k, — Re /+ m 6 fci0 6; r e i[ ( a ''=- w ") < - u;fc ( 5/ ^ + ^( m / /3 )lad'' 

We know bk , 0 , 6* r , and e i[(^fc-“n) £ - w fc(3/^)+ Wn ( m //3)] are complex. We will calculate k in 
Eq. (131) using the identity 

Rc[Z] = i (Z + Z‘) 

Since (to)* = u*v*, by extension it can be shown that 

(uvw)* = u*v*u >* 

Also note that if S — a x + a 2 + ge 3 + a 4 -\ — ■ where cq, a 2 , a 3 , a 4 , • • • are complex, then 
i?e[5] = i?e[ai] + Re[a 2 ] + Re[a 3 ] + /?e[a 4 ] + • • • 

We will also use the fact that if z is complex and a is constant, then 

Re[az] = aRe[z\ 


Using the facts above, we can see 

K — J_ [a. A* e j{(vk-Un)t-w k (g/P)+uj n (Tn/l3)] 

n,T 

+ i t, 0 ^,re" i[( “*" u ' n)1 ""' [(sW+u " (m//S)1 ] p 3+m ad' 

Then 

K = Q b* e M u/ k-u n )t-u k (g/f3)+uj Tt (m/p)] 

+^,o i «,r e " j!( ° Jfc ' Wn)i " Wt(3//5)+ “ n(m//3)1 ] p 5+m ad / 
with a and d' defined by Eqs. (116) and (121), respectively. 


(132) 


( 133 ) 



Now, we will turn to computing 2ile[AB*]. This will have five summations just like 
the 2 Re [AC*] term above. 

+M/ 2 q-l +M/2 -l (-rg)-l 

2Re [AB*] = E E E E E < ( 134 ) 

fc=-M/2 3=0 n=—M/2 r =- Lm =-(l+r) q +l 
k^O n^O 

By inspection of Eq. (123) 

q - Re [p 3+m 6 fci o6; r e j[ ( Wk_w " ) ‘ _a ' A: ( fl//3)+a ' n (m ^ ) ] aC / 

After letting the 1 /2 generated when applying Re[-} operator to the argument cancel 
the 2 outside of the summations, we get 

C — [ o /c,0 o „,r e 

_j_£« ^ re -j[(uk-Un)t-uk(g/P)+un(m/f3)] p9+ m ac ' (135) 

With a and c' defined by Eqs. (116) and (119), respectively. 

We now compute the term 2i?e[BC*]. This term will have six summations. Moving 
the Re[-] operator inside of the six summations to the argument yields the following 
by inspection 

i = Re [b ki eb* n r e j[(wk "^ )t “^ (9//3)+Wn(m//3)1 p s+m cd'] (136) 

After letting the 1/2 generated when applying the Re[-] operator in l [Eq. (136)] 
cancel the 2 outside of the summations, we can then write down the expression for 
2i?e[BC*] 

+M/2 -1 (-tq)-l +M/2 -1 -(r-l)q-l 

2*e[BC*] = E E E E E E ‘ (‘37) 

fa-M/2 (1+0.+1 „ = _ M/2 r—i 

Ar^O n^o 

i = ^b k/ bl tr e j[iuk ~ Un)t ~ Wki9/p)+Un{m/m + b* k>l b n>T x • • • 

g-j[(w fc -Wn)t-Wfc (<?//?) (rn//3)] 1 pg+m c( j/ 



with c and d' defined by Eqs. (118) and (121). 

Now we must find 

/ K 2 \s A ci\ 2 dt = K 2 [ AA*dt + K 2 [ T BB'dt + I( 2 f CC'dt 
Jo Jo Jo Jo 

+K 2 jf 2 Re [AC*] dt + I < 2 £ 2 Re [AB*] dt 

+I< 2 £ 2Re[BC*]dt (138) 

where we recall that 

K 2 = P(1 - p ) 2 


In the term A A* [Eq. (128)], we have the pulse multiplication aa' where 



Arbitrarily assuming m > g, we sketch the product a x a'. 



I 1 => Product is 0 
V/} =$■ Product is 1 


The product is 1 and exists from the maximum of {m/[3,g/p) to T. Either g or m 
can be the largest. Hence, we will be integrating 


£ = e Muk-u n )t-Uk(g/P)+Un(m/P)] 


(139) 


with the lower limit of integration being (1//?) max(< 7 ,m). The upper limit of inte¬ 
gration is T. We consider two cases for Eq. (128). 
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u>k = kAu 


Case 1: k ^ n => ujk ^ u n 
When k ± n, u k ^ u n because 


(140) 


where again u>k is the radian frequency spacing between Channel 0 and Channel k, 
and A u is the uniform radian frequency spacing between adjacent channels. 


K 2 r AA *dt = 
Jo 


+M/2 q -1 +M/2 q-l 

P(i-pf £ £ £ £ P’ +,n h,oK,c x 


k=-M/2 9 ~° n=-M/2 m ° 


i 


(l//3)max(5,m) 


k^O 

£dt 


n^O 


(141) 


Letting 


X = j [(wfc - u n )t - u k (g//3) + u n (m/P)\ 
dx = j[{u k -u n )]dt 


then 


L 


(l/P) max(g,m) 


£dt = 


p j[(vk-Vn)t-Uk(g/P)+Vn{m//3)] 


j(u>k — LO n ) 


1 r e i[(w*-ui„)T-Wfc(g//?)+w„(m//3)] 

j(oJk -w„) 1 

_ e i[( w fc-w n )(l//?) max(5,m)-w fc (5//3)+cj„(m//?)] 


(l//3)max(s,m) 


(142) 


Case 2: & = n =>■ uifc = u> n 

Again similar logic to that used above [Eq. (140)] allows us to conclude Uk = w n . For 
this case, the complex exponential £ reduces to 


£ _ e i[0-t*>fc(s//3)+u<„(m//3)] 


(143) 
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Now, since Uk = u>„, we can further reduce the complex exponential to 


£ = e- 7 ^ 


j[{uklP)(™—g)\ 


( 144 ) 


Therefore, for u>k = 


/ (dt = e ^IP)(m-a)\ / dt = e ;[K//3 )(m-g)\ T _ m ^(g, m 

' (l/P) meLx(g,m) J(l/p) ma.x(g,m) p 


Thus, we may write the expression for K 2 / 0 T AA *dt 


+M/ 2 9 -l + M/2 5—1 


A' 2 / AA-rfi = P (1 - p ) 2 E E E Ef'*’** 1 ' ( 146 ) 

JO - r\ -n 


k=—M/2 9 ° n=—M/2 m_ ° 


where 




l* 


_^i[(‘* , fc- u, n)(l//3)max(g,m)-a-/ k (5//3)+a- n (m/)3)] 


for A: / 


6 fc ,o^, 0 ^)] T — ^ max(j, m) , for A; 


Now we turn to the BB* term [Eq. (129)]. We have the pulse multiplication cc 1 where 


c = Pt 


mt) 

g + (1 +£)q 


c =Pt 


m + (1 + r)q 


Arbitrarily assuming m + (1 + r)q > g + (1 + l)q , we sketch the product c x c'. 
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) 

□ 

m 


Product 

Product 


is 0 
is 1 


The product exists and is 1 from 0 to minimum of 

(9 + (!+% m + (l+r)q\ 

l 0 ' fi ) 

Hence, we will be integrating the same complex exponential ( [Eq. (139)] as for the 
AA* term [Eq. (128)], but the upper limit of integration will be 

min (g + (1 + L)q, m + (1 + r)q) 

while the lower limit is 0. We again consider two cases. 


Case 1: k ^ n =$■ u>k ^ qj n 

After interchanging the order of integration and summation, we will be integrating 
[see Eq. (129)] 


£ _ e jl(vk-Vn)t-wi'(g/j3)+w n (m/0)] 


(147) 


Then 


/•(l//9)min(g+(l+r) ? ,m+(l+r)7) 

/ tdt 

Jo 
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|(l/0) min(g + (l+^)<?,m + (l+r)g) 


j(uj k - UJ n ) 


j(uj h -u n ) 


o j[(uk-Vn)t-u>k(g/P)+Wn{m/P)] 


J[(uJ k -u> n ){l/P) mn(g+(\+t)q,m+(\+r)q)-w k (g / (3)+w n (ml (3)] 


_J[-Uk(g/P)+Un{m//3)] 


( 148 ) 


Then, factoring the result yields 


r(l/P) min( 3 +(l+£)g,m+(l+r) 9 ) 


£dt = 


j(u k - u n ) 


J[-“k{g/P)+u„(m/0)} e i[(w fc -w n )(l//3)min(s+(l+^)ij,m+(l+r) 9 )] _ (149) 


Case 2: k = n => to = cu n 

This is similar to the AA* term for w* = [see Eq. (144)] 


£ — e i[(^fc//?)(m-a)] 


(150) 


/■(l/(J)min(g+(l+/) 9 ,m+(l+r)<j) /-(l/PJnunlg+ll+flg.nH 

/ £ dt = e J[K//9(™-s)l / 

Jo Jo 

= e^K/W"*-*)] i min (5 + (1 + *) 9 , m + (1 + r) 9 ) - 0 


f(l//3)min(s+(l+%,m+(l+r)g) 


(151) 


Thus, we may write the expression for K 2 / 0 T BB *dt 


j' +M/ 2 —1 (—^g)-l +M/2 —1 (—rg) — 1 

A' 2 / BB'<A = p(i—p) 2 E E E E E E 


fc=-Af /2 
A: 5^0 


_-L 5 T— — (lH-^)g + l r ~-L m=-(l-j-r)g+l 


( 152 ) 
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where 


h,(K,r 


Jl-Uk(g/P)+Vn(m/0)] 

j(u k - UJ n ) 


i[(uj k -u ) n )(!//?) min(s-f-(l+%,m-f (1+ 


r)?)1 -l) , 


for k ^ n 


T2 — 1 

= h,tK,r ^eM WklP ^ m ~ 9) } min (<7 + (1 + t)q, m + (1 + r)q) ^ , 

for k = n 

Now we integrate the CC* term [Eq. (130)]. We have the pulse multiplication dd ; 


where 


W,_ A+£?)'\'\ l T (t- 


d' = p T 


L £ 

m + rq 


d _ B d '-„ li_L J-Jl 

Looking at the pulse multiplication graphically where we arbitrarily assume m + rq > 


g + tq for the purpose of sketching the situation, we have 



9±?q 

P ' 


m + rq 

~T~ 


□ =£■ Product is 0 
|/V| =» Product is 1 


Looking at the equation for CC* [Eq. (130)], we see that after interchanging the order 
of integration and summation, we will be integrating the same complex exponential 
as for the AA* and BB* terms 




i = e- 


(153) 
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where the lower limit of integration will be 


1 

- max(# + tq, m + rq) 


while the upper limit will be T. Then for 
Case 1 : k ^ n => uik / c o n 

rT 1 

idt = — -t 1 

(1/p) max(g-\-£q 1 m+rq) Jy^k ~ ^n) 


L 


l 


j(u k - UJ n ) 


j[(u> k -u; n )t-w k (g/P)+tj n (m/P)] 

\(l/P) ma.x(g-\-£q,m+rq) 
J[(uk-Vn)T-uj k (g/P)+v n (m/P)] 


_ e j[(vk~Un)(l/P) m3Lx(g+£q,m+rq)-w k (g/P)+u n (m/P)] 


(154) 


Case 2: k — n => u>k = lo u 

Again, similar to the AA* term, for k = n, the term 


_ pj[(wk-w„p-Wfc(g//3)+wn(m//3)] 


£ = e 


(155) 


reduces to 


— J[(“k/0)(m-g)} 


( = e 


(156) 


and then 


[ T (dt = e j[(wfc/ ^ )(m_5)] [ T dt 

J(l/0)ma.x(q+tq,m+rq) J( 1 


— e i[(“fc//3)( m -3)] 


'(1//3) max(g+* 5 ,m+rg) 

T - -t max(# + tq, m + rq) 


(157) 


Thus we may write the expression for K 2 / 0 r CC *dt 


+M/2 _i -(*-l) 9 -l +M/2 -1 -(r- 1 ) 5-1 


a : 5 


/ cc** = p(i - E E E E E E / + ” 

jur/n 5=“^ _»>r/n m ——rq 


k=—M/2 
k^O 


n=-M/2 
n^O 


X </?3 


(158) 
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<^3 - < 


b k,eb* n>r 


J[(u k - Wn )T- (Jfc (s//3)+u>n (m//3)] 


i(^fc - w„) 

^ e 3[(vk-Un)(l/P)m<ix(g+tq,m+rq)-uj k (g/l3)+LJn(Tn/P)] 


for k ^ n > 




1 


T - - max (5 + m + rg) 


for A: = n 


Now we turn to the 2/?e[AC*] term [Eq. (131)]. We have the pulse multiplication of 
ad' where 

I 

m + rq\\ \ 



frit 

9 ’ 


/ 

Pt 

V 

fij 


d' = pt 

(T- 





\ 

f>) 

) 

' \ 


W 




('- ( t 3 )) 


Looking at this graphically by arbitrarily assuming m + rq > g, we sketch the product 
axd' 



□ =* 


Product is 0 
Product is 1 


Then we can see that the product ad' will exist from (1//3) ma x(g,m + rq) to T, which 
are the lower and upper limits of integration, respectively. 

Case 1: k ^ n u>k / u n 

This expression has the two complex exponentials 


c 


j[(w*-w„)t-w fc (s//3)+u;„(m//3)] 


(1 cm 




~ p-j[(vk-wn)t-<jj k (g//3)+u n (rn/l3)] 


g = e 


(160) 


to integrate. Letting 


u = ~j [(w* - w n )t - u k (g/fi) + u n (m/p)] 


du = -j(u k - w„) 


the integration yields 


/(1/P)max(g,m+Tq) 


[~j(Wk -<*>„)] 


-j(u k - UJ n )] 


D -j[(uk-Un)t-Uk(g/P)+U n (m/P)] 


(\/P) maLx(g,m+rq) 


~j[{vk-u n )T-Lu k (g/P)+tu n {m/0)] 


o-3[{uk-Vn)(l/P)ma.x(g,m+rq)~-w k (g/0)+w n (m/0)] 


(161) 


r(l/0)ma.x(g t m+rq) 


£dt = --- e 


j[(vk-Vn)T-w k (g/0)+w n (m/0)] 


j(u k — CJ„) 1 

^ P 3\(u k -un){ 1/0) max(s,m+r< 7 )—u/ fc {g/0)+w n (m/P)] 


(162) 


Case 2: k — n => u k — u n 

Looking at the expression for 2 Re [AC*] [Eq. (131)], we see that for u k = u n the two 
complex exponentials reduce to e^ Wk /^ m ~ 9 ^ and which are constants 

to the integration. Thus, we will integrate the product of the pulse functions ad' from 
0 to T. This is the only time function for this case. Then, looking at the sketch below 
Eq. (158) 

rT rT J 

/ ad'dt — / dt = T — — max(a, m + rq) (163) 

JO J (1//3) ma.x(g,m+rq) /? 


(163) 


Thus, we are ready to write the expression for 


I < 2 [ T 2Re [AC*] - 
Jo 
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Then, 


+ M/2 g-1 + M/2 -1 — (r — l)g— 1 


A' 2 / 2Ae [AC*] di = P(1 — p) 2 E E E E E ^ ( 164 ) 

- '° , 3=0 __ r=-L m=-rq 


where 


9*4 = < 


P 9 + m K oK,r 


]{u k — LO n ) 


k=-M /2 n=-M/2 

k^O n^O 




0 j[(^A:-^n)(l//?) max(g,m+rg)-u; fc (g//?)+a» n (m//?)] 


, g + mL* £ [_1_ [ e -i[(wfc-Wn)r-w fc (g//?)+a; ri (m//?)] 

i P /c,0 n,r [[~;(u»fc-u;„)] 


_p-3[(wk-Vn)(l/P) ma.x(g } m + rq)-uj k (g/p)+to n (rn/(3)] 


, for k ^ n 


'b k flb* k)T e j[M(3)im - 9)] + &fc,o6fc,re" J ' [( ^ //?)(m " s)1 ] / +m x • • • 


T- 1 - max( 3 -, m + rq) 


iox k — n 


We now turn to the 2Re [AB*] term [Eq. (134)]. We have the pulse multiplication 
ac' where 

( \ 

~ ' 9 


& = p T 


rp £_ 

W ft 


9\\ Pj 


C = Pt 


pm 


,m + (1 +r)q / 

To form a sketch, we arbitrarily assume m+ (1 + r)q > g , yielding 


/ 
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| | =£■ Product is 0 

Y/A =£• Product is 1 



We see that the product will exist for m + (1 + r)q > g and will be 0 otherwise. Thus, 
unless m + (1 + r)q > g, the integral of 2 Re [AB*] is 0 . We will later create a “gating” 
function to account for this. 


Case 1: k ^ n =$> lo ^ u n 

For this case where m + (1 + r)g > g we will be‘integrating the same complex 
exponentials as before 


£ __ e 3[(u k -w n )t-w k (glP)+w n (mlP)) 


(165) 


and 


Q — e -H( UJ k-^n)t-u k (g/0)+w n (m/l3)] 


(166) 


We have computed these integrals before and will just substitute the new upper and 
lower limits in the final answer for this term 

n t • •. (m + (l + r)q) 

P 

Lower Limit : %■ 

P 

If m + (1 + r)q < g, the integrals of the two complex exponentials are zero as the 
product of the pulse functions ac' — 0 . Again, we will create a “gating” function to 
allow for this. 
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Case 2: k — n ==>• Uk = tu n 

Looking at the 2.fte [AB*] term [Eq. (134)]. We see that the two complex exponentials 
reduce to 

e i[{ukl0){m-g)) an( j e -j[(<*W0)(m-s)] 


These are constants to the integration with respect to t. Thus as before, we will only 
be integrating the product of the pulse functions ac' from 0 to T. The integral has 
value if m + (1 + r)q > g and is 0 otherwise. Then, if m + (1 + r)q > g , the limits of 
integration are the same as for the case Uk ^ u n and 

rT r(m-+-(l + r)q)/P \ 

( a c'dt= I dt = -[{m + (1 + r)q) - g] (167) 

Jo Jg/P P 

Now we will create the “gating” function to null the integrals if the product ac' is 0 
from 0 to T. 

To indicate the zeroing/nulling of all integrals if g > m + (1 + r)g, we will create 
another gating function called G 2 (< 7 ,m,r). 

( 1, for g < m + (1 + r)q 1 


G 2 {g, m,r) = < 


(168) 


[ 0, otherwise 


We are now ready to write the expression for K 2 

T +M/2 q-\ 

2Re[kW]dt = P(1 -p) 2 J2 E 

k=-M /2 5 ~° 
k? 0 



/ 0 T 2 Re [AB*] dt. Then 


+M/ 2 -l (-rg)-l 

E E E 


n=—M/2 ' = 
n^O 


^5 


(169) 
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where 


( 0 9+m bunb* _ - _ e i[(uk-Un)((™.+{l+r)q)/(3)-w k (g/P)+Wn(TnlP)\ 

V K,u n,r j( Wk -ui n ) [ 


0 j[(wn//3)(m-g)] 


4 . (n 9+m b* b _1_ \p-j[(uk-U7i)((m+(I+r)q)/0)-uJk(g/P)+u„(m/0)] 

+ \ P °W 0 n ' r [-j(u k -w n) } [ e 


-e j[(un/P){m g)] j G 2 {g,m,r ), for k ^ 

b kfi bl ir e j ^ Wk ^^ m - 9 ^ + bl 0 b k , r e- j[{uJk//s){m - 9)] ] p 3+m x • • • 


;[( m + (1 +r)q) - g] G 2 (g,m,r), 


for k = n 


where 


G 2 (g, m, r) 


1 , for g < m + (1 + r)q 


0 , otherwise 


Now, finally we must write out the last (sixth) term of the ACI. We now need to find 


I< 2 [ T 2 Re [BC* 
Jo 


Again, we look at the expression for 2Z2e [BC*] [Eq. (137)] to see the pulse multipli¬ 
cation cd' where 


c = p T 


jTjt) 

9 + (1 + tjq 


d' = p T 


T \t - 


m + rq 
~ 

m + rq 

~T ~, 


We sketch at this pulse multiplication by arbitrarily assuming g + (1 + l)q > m -f- rq. 
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0 0 


Product is 0 
Product is 1 


Case 1: k ^ n =$■ u>k ^ iu n 

Here we will be integrating the complex exponentials 


£ — e jl(vk-u> n )t-Wk(g/P)+u n (m/P)] 


(170) 


and 


g — e -jK“k-Un)t-u k (g/P)+Wn(m/l3)] 


(171) 


The integrals will have the following limits 


Upper : 
Lower : 


9 + (1 +l)<? 

0 


m + rq 

~T~ 


The integrals will have value for m + rq < g + (1 + l)q and will be zero otherwise. 
We have computed these integrals before and will just substitute the new limits in 
the final answer for this term. We will create another “gating” function to account 
for the integrals “turning ofF” when m + rq > g + ( 1 + t)q 

( 1 , m + rq < g + (1 + £)q 1 

G 3 {gJ,rn,r) = < > (172) 

( 0, otherwise J 

Case 2: k = n ==> u>k — 
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Looking at the 2 Re [BC*] term [Eq. (137)], we see again that the complex exponentials 
reduce to 

e j[(u k l(3){m-g)] an ^ e -j[(u»/fc/^)(m-g)] 

which are constants to the integration with respect to t. Thus, in a similar manner 
as before, we will only be integrating the product of the pulse functions cd' from 0 
to T. The integral has value if m + rq < g + (1 + £)q and is zero otherwise. We have 
already created the “gating” function [Eq. (172)] to account for this turning on and 
off of the integral. Then, if m + rq < <7 + (1 + £)q the limits of integration are the 
same as for the case ui^ ^ u n and 

[T r(g+(\+t)q)/0 1 

/ cd dt= dt = -[(g + (1 +£)q) - (m + rq)] (173) 

•^0 J (m+rq)/@ fJ 

We are now ready to write the expression for 

K 2 f T 2Re [BC*] dt. 


Then 




rp +M/2 

-i 

+M/2 

_1 -(r-l)g-l 

K 2 / 2 Re [BC*] dt = P( 1 - p) 2 J2 

E E 

E 

E E 

J o 

k=-M/ 2 

t=-L g=-(l+t)q+l 

n=-M/ 2 

r~ — L m — — rq 

k^O 


nji 0 



(174) 
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where 


n 3+ m /), ./>* I \ f ,j[(wk—Un)((g+O-+Z)<l)/l3)— l * J i s (g/P)+Vn(m/0)] 

P °k,l°n,r j{ Wk - u,„) [ e 


_ e j[(w k -Un)((m+rq)/p)-uj k (g/0)+u>n(m./())] 


Ij+ml* U _i_ \ p -j[(u k -u>n)((g+(l+e)q)/l3)-Uk(g/P)+Un(m/p)] 

+P [ 

_ e -j[(uk-un)((™+rq)/P)-uk(g/P)+un(m/ 0 )] ^G 3 (g,£, m,r), for 


b k>i b* k re i[(a,fc//J)(m_s)1 + bl e b k}T e- j[ ^ k/0){m - 9)] ] p 9+m x • • • 


for k = n 


-[(^f + (1 + £)g) - (m + rg)] G 3 (g,A m, r) 


where G 3 (g,£,m,r) is defined by Eq. (172). We have now completed integrating all 
six terms of the ACL 

Recalling the expression for \s(t)\ 2 [Eq. (89)], we see that we still have to compute 


the following terms 


2 Rc[s B Sj B j] — 2Re[s/5/s^], 

2Re[s* B s AC i] = 2Re[s B s* ACI ], 


2Re[s* ISI s A ci] = 2Re[sisis ACI ], 

multiply each by the constant K 2 = P(l — p) 2 and integrate each from 0 to T. We will 
use either form of the above three expressions depending on the ease of computation. 
Recall from Eq. (66) that 




for 0 < t < T (175) 


Let us look at s B s* ISI term. Looking at the expression for s B (t) during the interval 
0 < t < T [Eqs. (66) or (175)], we see that s B {t) is real as &o,i 6 {0,1}. Looking at 



the expression for Sjsi(t ) during the interval 0 < t < T [Eq. (78)], we see also that 
s isi(t) is real as again b 0l £ {0,1}. Then sjsi = s }si and s b s *isi = s b s isi, and 


2Re[sBS* ISI ] = 2 sbsisi 


(176) 


Performing the multiplication yields 
2Re[s B s* ISI ] = 2 [ sb - s / 5 /] 


9-1 -l (—*«)—i 

= 2EE E P s+m bo,obo,iPr 

g—Q i=-Lo 77i = -(l + i)g+l 


,_s_\ „ (_jm_ 

/}) T (m + (l+i)? 


g -1 -1 —(*— 1)9 — 1 

+2E E E P s+m 6 oA,pt 

g=Oi=-Lo m= — iq 


T t- 


T _ g_ 


m + iq 

OlZ 

m + iq 


(177) 


We see that when we multiply by K 2 and integrate, the integral will distribute across 
addition to the two terms e and f. The integrals then move inside of the summations 
to the pulse function products. Now we need to figure the proper limits for the 
integration of the two pulse function products in e and f [Eq. (177)]. We note 
that the products of pulse functions are the only time functions in the expression 
2 Re[sBS* ISI ], so they are the only integrands in the expression. We examine the pulse 
multiplication in term e [Eq. (177)]. Let 


E = pt 


T\t-i 


m + (1 + i)q 


(178) 


Sketching E, we arbitrarily assume m + (1 + i)q > g. 
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| 1 =$> Product is 0 

Y/) =$> Product is 1 


0 0 

The integral / 0 T E dt will have value for g < m + (1 + i)q and will be zero otherwise. 
If g < m + (1 + i)q 

[T /•(m+(l+j')i J )//3 i 

I Edt= Lt dt= p 1(m+ (1 + ,)q) ~ s] (179) 

Since only the pulse function product is being integrated, we will create another 
gating function, which also gives us the value of the above integral when it turns on. 
This is similar to Gi(g,i,m,r ) [Eq. (112)]. 

T f 4l(m + (l+ i)q) - g ], g < m + (1 + i)q j 

Edt = G 4 (g,i,m) = ^ P l (180) 

[ 0 , otherwise J 

Look at the pulse multiplication in term f [Eq. (177)]. Let 




Sketching F, we arbitrarily assume m + iq > g. 


(181) 


79 




We see that 

f T 1 

Jo ==T ~ max (#,™ + *'<?) (182) 

We can now write down the following integral 

r T 

K 2 J 2 Re[sBS* ISI ]dt 

Substituting integrals / 0 r E dt and Jq F dt for the pulse function products in the 
2f?e[s#sJ 5 j] expression [Eq. (177)], and multiplying by K 2 yields the following result. 

Note: 

rT 

J ^ E dt is substituted for the pulse function product in e in Eq. (177). 
fT 

F dt is substituted for the pulse function product in f in Eq. (177). 
Recalling that K 2 = P( 1 — p) 2 we have 

A 2 / 2 Re[sgs* lsi ]dt = K 2 f 2[sgs/s/]<A 
Jo Jq 

= 2P(\ - p) 2 Y^ Y Y p 9+m bo,ob 0ii G 4 (g,i,m ) 

5=0 «'=— Lq m=-(X+t)5+l 

9-1 -1 -(.-1)9-1 r , -1 

+ 2P(1 -pfY Y Y P 9+m bo,ob 0ii T - - max(#,m -f iq) (183) 

g=Q \ — 77i— — iq P 
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G 


q-1 +M/2 -1 

+ E E E E P 9 +m b 0 fl b k ^ k(t ~ (mlP)) 

g=0 k _._ M j 2 t=-Lm=-(l+e)q+l 
k^O 



H 


g-l +Mf 2 _l -(^-1) 9 -1 

+ E E EE / +m Mi/ w * ( ‘- (ro/ 3 )) x- 

3=0 k=-M/2 fc " L m= - tq 
k^Q 
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Pt 





Pt 



■v 

I 


\ 


/ 


(184) 


The pulse products are G, H, and I as shown above. Then we note again that 6 0 ,i is 
real. Also, b ^ £ {0, e ^ k } and are com pl ex . We then have 


2Re[s* B s A ci\ — 2Re[sBSACi] = 

9-1 +M/2 9-1 

E E E + 

3=0 k=-M/2 m=0 

0 


9-1 + M /2 _i (-€9)—1 

+E E E E (/ t "Vo 


5 ~° k~~M/2 m=-(14-%+l 

k^O 


9-1 + M /2 -1 -(^- 1 ) 9-1 

+ E E EE (/ + ”Vo bk t ee^ Wk ^~^ m ^^l + p 3+m & 0 ,o^e- J "* (t - (rn//3)) l) 


3=0 *=- M /2 m= -* 9 

A:^0 


(185) 

We see that computing K 2 / 0 T 27?e[sgS > ic'/] will involve integrating each of the six 


terms inside of the three clusters of summations. The integrable part of each of these 
six terms is a complex exponential of the form 


$ _ e juk(t~(m/l 3 )) 


(186) 


or 

H = (187) 

The pulse products G, H, and I will provide the limits of integration. 
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Looking at G [Eq. (184)] 


G = pt 


ru-4 


Tit- 


T _ 9_ 


(188) 


Sketching G, we arbitrarily assume m > g. 



Product is 0 
Product is 1 


9_ ™ 

P P 


Since the height of the pulse product G is 1 from (l//3)ma x(g,m) to T, we see that 
we will have the following limits of integration for the complex exponentials 


^ _ e jw k (t-(m/f3)) 


(189) 


^ - - e -j“k(t-(m/0)) 


(190) 


which are the integrands in the first cluster of summations in the 2Re[s* B SACi] term 
[Eq. (185)] 

Upper Limit : T 


Lower Limit : — max(g, m ) 


Turning to H [Eq. (184)] 


H — pt 


T\t-^ 


T _ 9_ 


pm 

m + (1 + £)q 


(191) 
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Sketching II, we arbitrarily assume m + (1 -f £)q > g. 




% 


Product is 0 
Product is 1 


We see that for g < m + (1 + t)q, the upper and lower limits of integration for the 
integrands = e i"W-(W/3)) and ^ _ e -j* k (t-(m/p)) contained in the sec0 nd cluster of 
summations in the 2 Re[s* B s AC i] term [Eq. (185)] will be 


Upper Limit : •^(m + (1 + £)q) 

Lower Limit : ^ 


We will now create another “gating” function to null the integrals of d = (<-(”»//3)) 

and g = e i n the second cluster of summations in the 2Re[s* B SACi\ term 

[Eq. (185)]. This nulling of the two integrals will occur if g > m + (1 + £)q. 

So we define 


f 

G 5 (gJ,m) = { 


g < m + (l + £)q | 


[ 0, otherwise J 

Now we compute the integrals of d = e j w kU-(m/p)) and g = 

upper and lower limits of integration 7 X and 


(192) 


for arbitrary 



J a\ 


ju k 


171 


on 


(193) 
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Similarly for /u 


fidt = e-’^-WMdt 

J ct\ J a\ 


_ L r -JVk(t-(m/0)) 

l-M) 


71 

Ofi 


Now we turn to the pulse product term I in Eq (184) 



Sketching I, we arbitrarily assume m + £q > g. 


(194) 


(195) 



We see from the sketch that the upper and lower limits of integration of the integrands 
$ — and fi = in the third cluster of summations in the 

2/?e[sgS4c/] term [Eq. (185)] are 


Upper Limit : T 

Lower Limit : — max(^, m + £q) 

Now that we have the three sets of limits of integration for each of the clusters of 
summations and the indefinite integrals of d and fi, we are prepared to write the 
expression for 

K 2 f 2Re[sgSAci]dt = K 2 f 2 Re[sgs^ci]dt 

Jo Jo 
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Then, recalling that K 2 = P( 1 — /o) 2 



where 

in which 

ip? 


in which 

98 = 

where 


lRe[s* B SACi]dt — K 2 [ 2Re[sBSAci\dt = + Cba ( 196 ) 

Jo 


q -1 +M/2 9-1 

Ajm = B(1 - «) 2 E E E»' 

o =0 . m =0 

y k=-M/ 2 


= P 3 +m & 0 , 0^,0 — 
ju k 


e ju k (T-(m/(3)) __ e i^((max( 5 f,m)-m)//?) 


+p 9 +m 6 0 i 0 6 ^ 0 - 


-M] 


e -jw fc (r-(m//8)) 


e -J“*r((max(g,m)-m)//3) 


9-1 +W/2 -1 M?)—1 

B B ^ = r(i-rt 2 i: E E E ¥>» 

5=0 k=-M /2 m— — (\+Z)q+\ 

k^O 


P 3+m bo fl bkj — 

ju k 


e juk((l+t)<l/ 0 ) _ e ju k ((g-m)l( 3 ) 




+^ + "Vo^j—, 


,-j“'k((l+%// 3 ) _ „-juk(.( 9 -™)// 3 ) 


m ) 

G 5 (g,£,m) 





£( < m + (1 + i)q 
otherwise 


9 -l +M/2 -1 -(*- 1 ) 9-1 

Cba = -P(l - E E E ^ 

5=0 k=—M /2 e =~ L m =~ lq 
kpo 
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in which 


<£>9 = p g+m b 0t0 bk,£- 


jv k 


1 


e jvk(T-(m/0)) _ gjw*((max (g,m+Cq) — m)/P) 

-juj k (T~(m/f3)) _ -ju k ((ma.x(g < m+Cq)-m) / 0) 


The last term we have to compute is I< 2 / 0 r 2ile[sJ SJ s^ C /]c?t. We have already shown 
that s/s/ is real; thus 

2i2e[s} s/ s > ic/] = 2Re[sisis A ci\ (197) 


As before, we will multiply s IS i and s AC i- 

Look at sisi [Eq. (78)], we will let the first term be A and the second be B. Looking 
at s A ci [Eq. (83)], we will let the first, second, and third terms be C, D, and E, 
respectively. Then 


sisiSaci — (A + B)(C + D + E ) — AC + AD + AE + BC + BD + BE (198) 


Then, we compute all six of the above terms individually 
-1 (—i?)-l +M/ 2 g-1 

AC = 52 E EE P 9+m bo,ih,oe jlJk(t - (mm X • • ■ 

i=Lo £=-(l+t')g + l kzz—M/2 m=0 


w (<,+a+oj w 


' ± 31 ' 

b- 5 )/ 


(199) 


-1 (—W)-1 +Af/2 -1 

ad= £ 52 EE E x-- 

i=—Lo 5= —(l+t)g+l k——M/2 7Tl — — 

A :^0 


( WO \ n f WO 

PT I — ) PT 


K g + (l + i)qJ \m + {l+t)q t 


( 200 ) 
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-1 (-»'?)-1 +M/1 _1 -(e-l)q-l 

AE= Y t E E E P 9+m bo,ih^ jwk(t ~ {mm x--- 

z=-L 0 s=-(l+t')g+l k=—M/2 l ~~ L 771 
k± 0 


n f W‘) v, 

PT L+d+oJ PT 




V 


m -f 

3~, 


) 


( 201 ) 


-1 —(*—1)9—1 +M/2 q-1 

BC= Y. E E E ^ + ”V.<'4.oe , '“* < ‘■ (m/ ‘ !,, X • • • 


i=-Lo g=-iq k= _ M/2 ™=° 


fc^O 


Px 


f iNizSI 

H’-t 1 )); 


\ /( m\\ 

T v-j) 


Pt 


m 


p ) 


( 202 ) 


-1 —(*—1)9—1 +M/2 -l (-£?)-l 

bd= y E EE E 

i——Lo g=-iq k=—M/2 m=—1 


Pt 



£z!D 

( 9 + 

K P JJ 


\ 

/ 


Pt 


( m ) \ 

\m + (1 + O?/ 


(203) 


-1 -(*-1)9-1 +M/2 -1 -(£-l) 9 -l 

be= Y E E E E 

i=-L 0 </=—.</ k=—M/2 i= ~ L m= ~ eq 

k^O 


Pt 


f ina l 

'-(*¥■)) 


m + ^\\ 

V V V f> J ) 


H 


(204) 
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When we distribute 2/?e[-] across each of the terms AC + AD d-and use Re[Z ] = 

[Z + Z*] 2 , the 2 will cancel with the 1/2, and we will be integrating the same complex 
exponential terms and g = e -JW W-( m /0)) as we did before for the 

2Re[s* B SAci] term [see Eqs. (186) and (187)]. Then, all we need to do is investigate 

the pulse products for each of the six terms AC + AD d-to figure out the proper 

limits of integration. Also, we may need to create other “gating” functions to null the 
integrals if there are specific index occurences which make the specific pulse product 
0 . 


The AC pulse product is 


Pt 


( mi ) 

\9 + (! +i)q) 



Arbitrarily assuming g + (1 + i)q > m, we sketch the pulse product 



| j Product is 0 
Y/A =* > Product is 1 


We see that the limits of integration will be 


Upper Limit: 
Lower Limit : 


+ (1 + i)q) 

m 

J 
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Now we define the appropriate “gating” function 


G 6 {g,i,m) 


< g + {l + i)q 


0, otherwise 


( 205 ) 


This nulls the two applicable integrals when the pulses do not overlap. 

The AD pulse product is 

+ + \ m +(1 + ^)9/ 

Arbitrarily assuming m + (1 + £)q > g + (1 + ») 9l W e sketch the pulse product. 



□ =r> Product is 0 
V/\ => Product is 1 


9 + (1 + i) q rn + (l+£)q 

~~ P V 


We see that the limits of integration will be 


Upper Limit : - min(^ + (1 + i)q, m + (1 + i)q) 

Lower Limit : 0 


The AE pulse product is 


B ( mi \ 

T U + (! +0?J ?T 


' m + iq 

PEL 

m + tq s 

~T~, 
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Arbitrarily assuming g + (1 + i)q > m + tq, we sketch the pulse product. 



| [ =£• Product is 0 

V/X =* > Product is 1 


m + tq ^ + (l + i )9 


We see that the limits of integration will be 


Upper Limit : 
Lower Limit : 


j{9 + (1 +i)q) 

jj( m + £ q) 


Now we define an appropriate gating function 


G 7 (g,i,£, m) 


1, m + £q < g + (l + i)q 

0, otherwise 


This nulls the two applicable integrals when the pulses do not overlap. 


The BC pulse product is 



Arbitrarily assuming g + iq > m, we sketch the pulse product. 


( 206 ) 
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1 


[ I ==» Product is 0 
Y/A => Product is 1 



m g + iq 

We see that the limits of integration will be 


Upper Limit : T 


Lower Limit : max(^ + iq, m) 

H 


The BD pulse product is 

/ 

PT 


T [t — 


'9 + iq\\\ 


(3 


V 


T- 


9±jq 
, S3 j 


PT 


(3T{t) 


,m + ( \+i)q y 


) 


Arbitrarily assuming g + iq < m + (1 + t)q, we sketch the pulse product. 



Product is 0 
Product is 1 





93 



We will now apply the methods outlined below Eq. (204) and perform the 2Re[-} 
operation on each of the six terms of s* isi saci — s isi s aci ■ We will then integrate the 
resulting d = e^^W-W/ 3 )) and ^ terms in each of the six clusters of 

summations and multiply by the proper gating function, if necessary, to compute 

K 2 [ 2 Re[s* ISI SAci]dt = K 2 [ 2 Re[sisiSAci]dt 
Jo Jo 

Recalling that K 2 = P{ 1 — p ) 2 

T j 1 

K 2 [ 2Re[s* ISI SAci]dt = A 2 / 2Re[sisi s ACi\dt = ACr+ADr+AEr+BCr+BDr+BEr 

Jo Jo 

where 

-1 (- t ' g)-l +M/ 2 q -1 

ac r = p(\- p ) 2 E E E Ev.o (208) 

t=-L 0 5 =-(l+t)q+l k~~M /2 m= ° 

Ar^O 

in which 

Vio = P 9+m bo,ibko -— 
ju k 

-\-p 3+m b 0 ib* k 0? — -—- e -^k[(i/P)(9+^+i)o)-m/P] _ x G 6 {g,i,m) 

’ ’ [-JUk] [ ^ 

where 

'1, m < g + (1 + i)q 

Go(g, i, m) = < 

0, otherwise 

* Note: = e j0 = 1 

-1 (-«<?)-1 +M/2 -i (-£g)-l 

AD R = P(l - pf E E EE E vn (209) 

i=—Lo g=-(l+«)<?+l k= _M /2 &=-£ m=-(l+*)g+l 

in which 

(p n — p 9+m boib k e - e i u k[^/0)knin(g+( 1 + i )'l<' m +( 1 + e )<l)-m/P] _ e jw k [-m/0\ 

’* ’ jWk L J 


0 j ^ k [( l /^)( s +( l + i )?)- m // 3 ] 


G & {g,i,m ) 
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, o+m l L* ^ r^-jw fc [(l//3)min(5+(l+i)q,m+(l+^)g)-m//3] _ e —juj k [-m/0\ 

+p boA ’ e H^]\- e J 


-1 (-tg)-l +Mf 2 -l -(*-l) g -l 

ae r = p(i-p ) 2 Y. £ £ £ £ 

izz — Lo g--(l+i)q+l k -_^fj2 771 ~ 


( 210 ) 


in which 


<^12 = P 3+m bo,ibk,e -— e- 


,jwic[(l/(3)(3+(l+t)q)-TO//3] _ p jw k [iq/p] 


+p 9 +m b 0 ibie —-— [ e -' J ' u "=t(i/^)(fl+(i+‘)9)-™/^ _ e -WV/3]l G 7 (g,i,£,m ) 
’ [-JWfc] 1 


where 


G 7 (g,iJ,m ) = 


1, m + iq < g + (1 + *)? 
0, otherwise 


in which 


_1 -(t-l)g-l +M /2 9-1 

BC„ = P(1 - ri 2 £ £ £ £v !3 

t =-Lo g=-ig k——M/2 m=0 


( 211 ) 


„ 9 +mh l 1 Ljwjfe[T-(m//?)] _ p ju k [(l/{3)max(g+iq>7n)-m/l3] 

P VQ,i°k,Q- e c 

]U k L 




e -jo>fc[T-(m//3)] g—jW[(l//?) ma,x(g+iq,m)—m/P] 


_1 —(1—1)9—1 +A//2 -1 (-*?)-l 

bd r = p(i-p ) 2 EE E E E ( 212 ) 

i=-L 0 g=-*g k=-M /2 <=-£ ’ m=_ (l+*)'J+ 1 


in which 


SP14 o,* L J 

+/ +m 6o ,-6* r L-^[(l//5)(l+^)9l _ e -jc*[(l/^( fl +i,)-m//3]j Gg ^ j ^ £) m ) 

’ ’ [-JWfcJ L J J 
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where 


G 8 (g,i,l,m ) = 


1, g + iq < m + (1 + i)q 
0, otherwise 


-1 -(i-l)g-l +M/2 _l -((-l)q-l 

BEn = P(l-p) 2 E E E E E ¥>18 


i=-Lo g=-iq fc= _ M/2 t=-L m=-lq 


in which 


( 213 ) 


V15 : t> 3+m b 0 ,{b k: e -— e 


Ju k [T-(m/P)] __ e jw k [(l/P)ma.x(g+iq,m+£q)-m/l3] 


+P ,+m bo.,K 


l-JUkl 


e -jvk[T-(m/0)] _ e -ju k [(l/p)m<ix(g+iq,TTi+£q)-rn//3] 


We have now computed all terms of 


f \s(t)\ 2 dt = K 1 [ s\t)s f (t)*dt [see Eq. (89)] (214) 

Jo Jo 


Now we recall from Eq. (89) that 


f \s(t)\ 2 dt = K 2 [ s f (t)s (t)*dt — K 2 [ \s'(t)\ 2 dt 
Jo Jo Jo 

= K 2 [ \ss\ 2 dt + K 2 [ \sisi\ 2 dt + K 2 [ \sACi\ 2 dt 
Jo Jo Jo 

+ I < 2 [ 2 Re[s B s* ISI ]dt + I < 2 [ 2 Re[s* B s AC i}dt + I < 2 f [s* ISI s AC i]dt (215) 
Jo Jo Jo 

We will now extract the answers for the above six integrals from the previous rather 


lengthy derivation so that we will have the final answer for 


[ T \s(t)\ 2 dt = K 2 [ T \s’(t)f 
Jo Jo 


( 216 ) 


in one single place. Now, before we consolidate the terms we once again restate 

max(xi,x 2 ): Choose largest of x\ or x 2 , which are both positive. 

If Xi = x 2 , then max(x!,x 2 ) = xi — x 2 . 


min(x 1 ,x 2 ): Choose smallest of Xi or x 2 , v ;ich are both positive. 

If Xj = x 2 , then min(xi,x 2 ) = xi = x 2 . 



Also recall that K = y/P (1 - p) and I< 2 = P( 1 - p ) 2 . Then we have 


|s(i)| 2 <ii = K 2 r s'(t)s'(tydt = K 2 jf T I s(t)\ 2 dt 


( 217 ) 


which has the following six integrals given by Eq. (215) 


<7-1 9-1 


k 2 [ T \s B \ 2 dt = p(i - pfbl 0 E E p s+ " (r - 

^0 5=0m=0 V 


max(< 7 , m) 

0 


(218) 


K 2 f \sisi\ 2 dt 

Jo 


we recall the result of Eq. (117) 

A' 2 J T s,s,s’,si dt = K ‘ [ \‘isi\ 2 dt = K 2 £ Ml + £ Bdt £ Cdt + £ Vdt 

(219) 

T -1 (-ig)-l -1 (-T-g)-l 

k 2 Adt = p(i-p) 2 E E E E vo 

® i——Lo g=— (l+t)9+l t=—Lq m=— (l+r)g+l 

in which 

o+m j , (mm(g + (1 + i)q,m + (1 + r)q)\ 

¥>0 = Oo.fOo.r I -^- I 

K 2 F Bdt + K 2 [ T Cdt = 2 1< 2 [ T Bdt = 

Jo Jo Jo 

-1 (-*?)-1 -1 —(r—1)9-1 

2P(1 -/>) 2 E 2 / +m 6o,,-6 0) rG ! i(^,i,m,r) 

i=-L 0 p=-(l+t)g+l r=-L 0 m~-rq 


where 


Gi(^,i,m,r) 


-[(# + (1 + 0?) - (m + rq)}, for m + rq < g + (1 + i)q 
— \ J 

0. otherwise 


-1 -(*'- 1 ) 9-1 -1 —(7—1)9— 1 


/f 2 / d* = p(i - p ) 2 £ E E E p s+ ” i 

•'° t=—Lo g=—19 r——Lo m=—rq 


ma,x(g + iq,m + rg) 




We know from Eq. (138) that 


I< 2 F \s A ci\ 2 dt = I< 2 [ T AA*dt + K 2 [ T BB *dt + I< 2 C CC'dt 
Jo Jo Jo Jo 

+I{ 2 [ T 2Re[AC*]dt + K 2 [ T 2Re[AB*}dt 

Jo Jo 

T 

+K 2 I 2Re[BC*]dt 

Jo 

So we have the following six ACI components 

j-£ + M/2 q—\ + M/2 q— 1 

a - 2 / AA'dt = p(i — pf Y. E E E/’ J+m *v. 

J0 , 9=0 m=0 


k——M/2 y n= — M/2 


k± 0 


n^O 


where 


<P i 


^fc,O^n,0 • 


1 


j(u>k - UJ n ) 


0 j[{uk-u n )T-u> k (g/(2)+w„(m/(3)] 


= < 


_ ( =j[(wfc-Wn)(l//3)max(3,m)-a-* ; (s//?)+a; n (m//?)] 


for k ^ n 


bkflb* k 0 e^ u k/P)( m [T — ^ max(j, m)] , for k = n 

l Kol 2 


( 220 ) 


( 221 ) 


K■ 


+M/2 

-1 M?)-l +M/2 

-i (— r ?)~ i 

E 

E EE 

E E 

:=—M/2 

te-t »=-(!+<),+1 „. M/J 

t — ~L m=-(l+r) 5 -(-l 


n^O 

(222) 


where 


h,fk 


n,r 


" e j[-uk{g/P)+Vn(m/f3)\ 
j(u> k -UJ n ) 


^ e i[( u ’*-‘ i '")( 1 /( 9 ) min (5+(l+^)9.m+(l+ r )9)] 



V^2 


for & ^ n 


b k /b 


fc,r 


oi[(wfc //5)(?rz—^)] 




min(# + (1 + ^)g, m + (1 + r)q) 


for k = n J 
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x <P3 


+M/2 -1 —(<—1)9—1 +M/2 -1 —(r— 1 ) 9 —1 


a 2 / ccdt = p(i — p) 2 y. E E E E E e 


kzz-M/2 

k^O 


= -M/ 2 
n^O 


where 


,,7 [(<*>* -w n )T-wic ( 9 / 0 )+u>n (m//?)] 


h,eK,r —/ - 7 | e 

’ [j(u k - 0J n ) L 


_ J _ e j[(w fc -a-„)(l//3)max(9+<9,m+r9)-w*(ff//3)+w n (m//3)] ^ £ or fc -j. n 


bk tb*k e^ w klP){ m -9)] T — — max(<7 + £q, m + rq ) , for A; = 


(223) 


+M/2 9-1 +M/2 -1 —( r — 1 ) 9—1 


A 2 / 2Ae[AC*] = F(1 — p) 2 E E E E E 

JO n — n r —— f, rn = — 1 


k——M/2 3=0 n—M /2 r= " L m= ~ Tq 


(224) 


where 


l* 


j[(uj k -w n )T-u; k (g/(3)+u n (m //?)] 


_ e j[(vk-Vn)(lf0)max(g t m+Tq)-uJk(g/0)+Vn(rn/0)] 




-i[(wfc“Wn)T-W fc (5//?)+^n(m//?)] 


_ e “'i[(^-^n)(l//?)max(5,m+r(7)-Wfc(^//?)+u;n(m//?)] ^ £ or fc ^ n 

^,0^it,r ej[(U ' fc//3)(m-5)] + &fc,o^,re _jt(U,fc//3)(m_3)] ] /> 3+m X • • • 


T - —ma,x(g,m + rq) , 


for k = n 
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+M/2 q -1 +M/2 -l (-rq)-l 


f 2i?e[AB*] = P(1 - p) 2 £ £ £ £ £ 

^0 _n _r I * , 


k=-M/2 9=0 n=-M/2 r=-£m=-(l+r)9+l 


( 225 ) 


where 


p 9 + m b k 0 b* _ \ e 3[(uk-Vn)((m+(l+r)q)/P)-Uk(g/P)+u ri (m/P)] 

’ n,T j(u>k -w n ) L 






-i[(w k -a-n)((m+(l+r)g)//?)-u;fc (g / ())+w n (m / (3)\ 


_ e -j[(«n//J)(m-p)]l ) G 2 ( 5 , m , r ), 


for k ^ n 


&M^,r ei[K/ ^ )(m ~ 9)] + 6fc 0 6 iir e- j ^/ /3 )( m - 9 ^] /3 9+m x • 


where 


-[(m + (1 + r)g) - #] I G 2 (^, m, r), 


<* 2 ( 0 , m,r) 


for g < m + (1 + r)q 


0, otherwise 


for k = n 


f j +M/2 -1 (-^g)-l +M/2 _1 -(r-l)cj-l 

K 2 J 2i?e[BC*] = P(1 - pf E E E E E E ¥>« 

0 k=-M/ 2 *=- L 9=-( l+0«+l n= — M/2 r= ~ L m = ~ rq 


( 226 ) 
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where 


i e 


JK^k-Wn)((g+(l+t)q)/P)—Wk(g/P)+w n (m/(3)] 


P 3 [(— Un ) ( ( rn+r q ) /( 3 )—Wk (g /P) +^n (m //? )] 


, n-\-m L* L ^ \ r ~ if(wfc—■^n)((p + (l+^)g)//5)—Wfc(^//?)+tJ n (m//?)] 

_g - it("*-Wn)((»n+rq)/^)-u"*(fl// 3 )+u<n(m/^)] ^ £, 777 , 7 ), for k ^ tl 


^,^/t,r ej[( ^ //3)(m ’ S)1 + fyfcA.r 6 ”^^” - ^ 1 ] ^ +m X 


-[(c/ + (1 + £)g) - (m + r?)] J G 3 (< 7 , t, m, r), 


where G 3 (< 7 ,£, m,r) is defined in Eq. ( 172 ). 
Now, for the last three terms of 


for k — 


which are 


[ T \s(t)\ 2 dt = /i 2 [ T s'(t)s'(t)*dt = A " 2 [ T \s'(t)\ 2 a 
Jo Jo Jo 

t »r 

I < 2 [ 2 Re[s B s* ISI ]dt] I < 2 2 Re[s* B s AC i]dt 

Jo Jo 


K 2 [ T 2 Re[s* ISI s A ci]dt 
Jo 

T «j 1 

K 2 f 2 Re[s B s* ISI }dt = K 2 / 2 [s B s IS i]dt = 

Jo Jo 


2P{l-pfY Y 12 P 9 +m bo,obo,iG 4 {g,i,m ) 

g= 0 i=—Lo m=—(l+i)g+l 


q-i -1 —(*-1)9—1 r i 

+ 2 P (1 - p ) 2 Y 12 12 P 9+m bo,obo,i T - - max (< 7 , m + iq) 

g—0 i—-Lq m — — iq L r 


( 227 ) 
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where 


G 4 (g,i,m) = 


l 

P 


[(m + (1 + i)q) - g], 


0, 


g < m + (1 + i)q 

> 

otherwise 


T T 

K 2 / 2 Re[s* B s ac i\dt = K 2 f 2Re[sBSACi]dt — Aba + Bba + Cba 
Jo Jo 


q -1 +M/2 q -1 

Aba = ^(1 - p) 2 E E V 7 
9=0 k=-M/ 2 m=0 


in which 


ip 7 


p 9+m b 00 bko - \e 9UJ A T -( m /0)) _ f>j w k(( max (g >m)—m)/p) 

' ' ju k *■ J 

4 -n9 +m h h* - \ p -j“k(T-(m/(3)) _ -jw k ((ma.x(g,m)-m)/0) 

+f VoV„ hM| [' 


where 

<p8 


where 


9-1 +M/2 -1 (-iq)-l 

Bba = p{i-p) 2 E E E E V* 

g=0 k __ M i 2 e=-Lm=-(l+e)q+l 
kqt 0 


P 3+m bo,ob k ,, 


JUk 


J^ k ((l+t)q/P) _ pju k ((g-m)IP) 


G 5 {g,£,m ) 


+p 9+m b 0 ob* ke7 — \e~^A(WW0) _ e -ju k (( g -m)/p) 1 
’ ■— ju k \ L J 


Gs{g,i,m) = { 


1, 

0, 


C B a = P(1- P ) 2 E 

5=0 


g < m + (1 + £)q 

> 

otherwise ^ 

+M/2 -1 -(*-1)9-1 

E E E ^ 

k=-M/2 l= ~ L m= ~ tq 
A# 0 
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_ p 9+m boobk(. ~ \ C 3UJklT - (m/l3)) c 3 Uk({ma.x(g,m+eq)-m)/ 0 ) 


+p 9+m b 0 ob* k£ - e -i“k(T-(m//3)) _ e -ju k ((max(g,m+eq)-m)/0) (228) 

0,0 k ' l [-jUk\ ^ 


T rT 

K 2 f 2 Re[s* ISI SAci\dt = I < 2 / 2i?e[s/s/s > ic7]^ 

7o «/o 

= AC# + ADrAEr + + BDr + BEr 


-1 (— tqf) — 1 +M/2 q- 1 

/ic B = p(i-p ) 2 E £ £ E» 

i=-L 0 sr=-(l+«')g+l k=-Mj 2 m ~° 


(229) 


where 


V’io = p 


' 9+m boib ko — e jw * [(1//3)(ff+(1+i)?)_,Jl//31 - 1 G e (g,i,m) 

’ ’ [ *. 

+p 3+m b 0 ibl or ^ e -i«-*[(i//»)(»+(i+*»-m/« _ 1 G 6 (5, h m) 


where 




1, m<5r + (l + i)? 

0, otherwise 


Note: e ±i“k(rnlP-ml!3) _ e i0 _ 1 


-1 +A//2 -1 (-*?)-! 

ai»h = j p(i -/>) 2 EE E E E ¥»n 

i—-Lo g=-(l+i)q+l k=-M/2 i- ~ L m= ~( 1 +0?+ 1 


(230) 


in which 


<ai = /> 9+m MM— k 

L 


\f>jv k [{l/P)min(g+(l+i)q,m+(l+e)q)-rri/@] _ ju k [-m/0) 


t _ Q+mL L* ^ f iwfc[(l/^) min(c;+(l4-09i m +(l+^)9) — m //^] m /^] 

+'’ l e ] 
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( 231 ) 


-1 (-tg)-l +M/2 -i -(*-l) 9 -l 

ae r = p(i-p? ■£ E E E E v.2 

i= — Lo g= — (l+i)q+\ k — — Mf2 Tn-=.—lq 


in which 


^12 — p 9+7n bo,ibk/- £' 


Ju k [(l/(3)(g+(l+i)q)-m/0\ _ ju k [tq/0\ 


where 


% ,-6T, r 1 1 — ~ [ e -^t(l//5)(3+(l+>)9)-™//3] _ e -J^[^//3]| G 7 ( 5) m ) 

’ ’ l-JWjfcj 1 J . 


G 7 (g,i,t,m) 


m + £q < g + (1 + 


0, otherwise 


in which 


_1 —(« —l)?—X +M/2 q-1 

BC r = P(1-pY E E E Ev>!3 

i= ~ Lo 3 =- tq k=—M/2 m=0 

k^O 


(232) 


5 +mL L _ Liw*[r-(m//?)] ja;^[(l//3)max(5+?g,m)-m//3] 

P °0,i°k,0 ■ e e 

J^A: 


I .S+mJ ./,* 1 L-jw*[T-(m/0)] _ ^-jwfc[(l//?)max(5+tg,m)-m//3] 

6o ' ! ^«h Wfc ] L e 


-1 -(i-l)fl-l +M/2 -1 (-2*)-! 

BB R = p(i-rt 2 ■£ E E E E ( 233 > 

t=-L 0 3=-tqf fc—— Af/2 m=-( 1 +^)g+ 1 


in which 


<^14 = 


i^k[(i/iS)(i+f)9] _ < oj‘*'k[(i/0)(s+«'g)-"i//?] 


f .-—- [ e -i^[(i/^)(i+^] _ e -j« fc [(i/«(ff+w)-m/p]l G 8 {g,i,l,m) 

' ’ l-jwjfcj L J 
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where 


Gs{g,i,l,m) 


g + iq < m + (1 + £)q 
otherwise 



where 1Z is the responsitivity of the photodetector (A/W). This will be used to com¬ 
pute probabilities of bit error for the dense WDM system. 
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APPENDIX B 


LIMITED CASE OF THE COMPLETE MODEL 

(10 PROGRAMMABLE TERMS) AND RESULTING 
PROBABILITY OF BIT ERROR EQUATIONS 


In Appendix A we derived the complete expression for \s(t)\ 2 dt. Referring to the 
OOK receiver structure in Fig. 1, we note that the deterministic signal component 
of the random variable appearing at the output of the integrator is 

X = K [ T \s(t)\ 2 dt (235) 

Jo 

where 'll is the responsivity of the photodetector (A/W) and s(t) is the complex 
baseband output of the Fabry-Perot filter. Note that the photodetector detects s(t) 
and produces an output current Tl\s(t)\ 2 . We will call the integrator output the 
decision variable Y where 

Y = X + N (236) 

We note that N is a Gaussian random variable with mean zero and variance N 0 T 

N ~ N(0, N 0 T ) (237a) 


and 

N= C n(t)dt (237b) 

Jo 

where nit) is the postdetection additive white Gaussian noise and N 0 (A 2 /Hz) is the 
two-sided current spectral density of n{t). 

Now we will set up a fairly limited but realistic case for which we will program 
the appropriate terms derived in Appendix A and generate the probability of bit 
error. Letting 


L 0 = 1 L = 0 


<t>k — 0 




2 ivkl 
T 
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u>k = 27 xklfT is a special case of the radian frequency spacing between Channel 0 and 
Channel k. Recall again that I is the normalized channel spacing integer (7 > 0), and 
T is the data bit period (s). L = 0 and 4>k — 0 means that we will only model and deal 
with adjacent channel interference (ACI) from the bits in the adjacent channels which 
occur during the detection window 0 < t < T. Recall that the adjacent channels are 
symmetric in frequency around Channel 0, the channel of interest. L 0 = 1 means that 
we will model the effects of a single ISI bit trailing the detected bit of interest 7 0 ,o in 
Channel 0. <j>k = 0 means that the phase offset between Channel 0 and Channel k is 
0, or that all adjacent channels are bit synchronous with with Channel 0. 

Then, since <p k = 0, L = 0, and b kt( E {0, e j<t>k }, we let 

b—M/ 2 ,o = b-M/ 2 + 1,0 — ' ■ ■ — fr-i,o = b and b E {0,1} (238) 

b~ is the left adjacent channel bit pattern and all M/2 of the left channel bits will be 
1 or 0 simultaneously during detection window 0 < t <T. And again, since 4> k = 0, 
L = 0, and b k / € {0, e^*}, we let 

bM/ 2,0 = b M/ 2 _i t o = • • • = b lt o = b + and b + E {0,1} (239) 

b + is the right adjacent channel bit pattern and all M/2 of the right channel bits will 
be 1 or 0 simultaneously during the detection window 0 < t < T. Also, we know 
L 0 = 1, and bo,i E {0,1}, so using Eq. (1) or Eq. (41) we have 

bo -i E {0,1} and 6 0 , 0 G {0,1} (240) 

We define an ACI/ISI bit pattern 

= {fr,6 + , 6 o,-i} (241) 

Note that b ~, 6 + , and 6 0 _i are 0 or 1 with probability 1/2 yielding eight possible 
values of ip p . 
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Now let use define an ACI/ISI bit pattern set 0 


fp = {^i, 02, * ■ ■, V’s} (242) 

We can denote each individual element in the set x/j as xp p where p = 1,2, 3,..., 8 as 
there are eight possible values of xp p , and 


0 = {</> P } 


(243) 


where p = 1, • • •, 8. Now we define 


Xo(V’p) = X(xp p , 6 0 ,o = 0) or X evaluated at the current value of V’p with &o,o = 0 

(244) 

Xi(xp p ) = X(xp p ,b 0l o = 1) or X evaluated at the current value of xp p with 6 0 ,o = 1 

(245) 


It can be shown that the conditional probability of error for the dense WDM 


system given the ACI/ISI bit pattern xp p is [2, 4] 

1 1 \ 1 (V T -X^„y 

P (errors) = -Q + -Q I ^ ~ 


(246) 


where 


There are eight possible members of the set xj) 


x\, = 


Vr is the detection threshold where 



= 0 

0 

0 

V>5 = 1 

0 

0 1 

i>2 

= 0 

0 

1 

Xp6 = 1 

0 

1 

</>3 

= 0 

1 

0 

V>7 = 1 

1 

0 


= 0 

1 

1 

V’s = 1 

1 

1. 


(247) 


(248) 


V T = 


X 0 max +Xi min 


(249) 
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and 

X 0 max = max (X 0 (ip p )) (250) 

Hp} 

i.e., to find Xq max, compute all eight values of Xo [see Eq. (244)] and then choose 
the maximum, and 

Xi min = min (Xi(^ p )) (251) 

Hp) 

i.e., to find X\ min, compute all eight values of X\ [see Eq. (245)] and then choose 
the minimum. Then, by the law of total probability, the probability of bit error is 

Pb — P (error|z/>i) P(V’i) + P (error|^ 2 ) P{i> 2 ) H-b P (error|^ 8 ) P(ip 8 ) (252) 

Assuming all bit patterns are equiprobable 

P(M = P(M = ■■■ = P(i>s) = | (253) 


Then we can easily see that 


Again, we note that 


Pb = £ H P (error| xjj p ) 


x = n | s(t)\ 2 dt 

JO 


(254) 


(255) 


where / 0 T \s(t)\ 2 dt was computed in Appendix A and consisted of 21 terms (“clusters 
of summations”) [see Eqs. (218)—(234)]. Given our parameters Lq = 1, L = 0, <j>k = 0, 
and u>k = 2 nkl/T, 11 of the 21 terms drop out and we are left to compute in order 
ten to generate the probability of bit error. If we name these ten terms (“clusters of 
summations”) SUMx, SUM 2 , • • •, and SUM 10 , then for Lq = 1, L = 0, <f>k = 0, and 


u>k = 2irkI/T 


[ T \s(t)\ 2 dt = SUMi + SUM 2 + • • • + SUM 10 
Jo 


(256) 


and using Eq. (255) 


X = 1l{ SUMj + sum 2 + • ■ ■ + SUM 10 ) 


(257) 
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Now PT may be factored out of each sum. So we define 


SUM, = FZXQUOT,-) j = 1, - - •, 10 (258) 

So 

X = KPT( QUOTJ + ftPT(QUOT 2 ) + • • • + 71PT(QUOT 10 ) (259) 

and 

X = n PT[QUOT x + QU0T 2 + • • • + QUOT 10 ] (260) 

Also, letting QUOT x + QUOT 2 + • • • + QUOT 10 = QUOTSUM, we also see 

X = 71 PT[ QUOTSUM] (261) 

We can also say that, in general, QUOT x , QUOT 2 , • • •, QUOT 10 , and QUOTSUM are 
dependent functions of ip p and 6 0 ,o- Thus, we may write in general 

QUOT 1 (i/’p, vo), QUOT 2 (V> p , b 0 ,o), • ■ •, QUOT xo (0 p , Vo) 


and 

QUOTSUM (V> P ,M 

Now we can see that 


XMp) = ftPr[QUOTSUM(tf p ,Vo = l)] (262) 

X 0 {ip p ) = 7lPT[QUOTSUM(^ p ,6 o ,o = 0)] (263) 


Recalling that 


Ur 


A 0 max +Xi min 

2 


we let 


ip p max = Value of ipp which causes maximum value of 

of Xq (or maximum value of A" with &o,o = 0). (264) 
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and 


Vv min = Value of ipp which causes minimum value of 

of Xi (or minimum value ofX with 6 0 ,o = 1). (265) 


Then using Eqs. (249), (262)-(265) 


Vj = 


■ftPr[QUOTSUM(^ max ,6o,o = 0) + QU0TSUM(i/> p ^Ao = 1)] 

2 


(266) 


Now using Eqs. (246), (262), (263), and (266) 
P(e vror/ipp) = 


5 Q {[tS) (|Q uotsum (^-V« =')] 

QUOTSUM (V’p ,„».,Vo = 0) + QUOTSUM (</>„ = 1) 


+ 


Hi 


iipt\ 

7^f) 


QUOTSUM A max , b 0 , 0 = 0) + QUOTSUM (^ mm, 6 0 , 0 = 1) 


—[QUOTSUM(W Ao = 



(267) 


Thus, we see that P(error|0 p ) is directly related to the signal to noise ratio which we 
will call Z. 

„ upt m „ HF 

VNo 

Thus, when we compute a probability of the bit error graph, we will choose a suitable 
range of values for Z , a value of free spectral range-bit period product f3T, and several 
values of the number of adjacent channels M. For each of these values of M, we will 
compute a probability of bit error trace. To compute a point on a trace, choose a value 
of Z , compute all eight values of P (error|?/> p ), sum all eight values of P(error/?/> p ), and 
divide by 8 [see Eqs. (248), (254), and (267)]. The point is then plotted. 
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Recall that 


then 



q = (3T 

We will use this and the fact that [see Eq. (259)] 


(269a) 


(269b) 


X = n PT( QUOT,) + n PT(QUOT 2 ) + • • • + n PT{ QUOT 10 ) (270) 

V_ v _✓ v- - -/ '-„-' 

Term #1 Term #2 Term #10 

to compute the probability of bit error graphs. We also use Eqs. (256), (257), (258), 
(260), (261), (262), (263), (264), (265), (266), and (267) in this endeavor. 

We now proceed to write out the ten terms of X , which we will use to compute 
the probability of bit error graphs. Before we begin we recall that max (xa,x 2 ) 
means choose the largest of X\ or x 2 . If x 2 = x 2 , then max (xi,x 2 ) = x\ = x 2 . 
Also recall that min(xi,x 2 ) means choose the smallest of X\ or x 2 . If Xx = x 2 , then 
min (xi,x 2 ) = x x = x 2 . The first term to be programmed will be Eq. (218), Appendix 
A. Factoring out a T we get 

SUM, = PT( 1 - p)X,of XV +m (l - ( maX ^ )) (271) 

Using q = (3T and multiplying by 1Z we get 

Term #1 = nPT(i - „)X fi £ (l - (=^)) (272) 

With L 0 = 1 and q = f3T, the first term (“cluster of summations”) in Eq. (219), 
Appendix A, becomes 

PT- 1 pT—i /i \ 

P(l-P) 2 Y1 (V-i) 2 (^min(9,m)) (273) 

g =1 m=l \P J 

Factoring out a T and multiplying by 7 Z we get 

(3T — 1 (3T— 1 

Term #2 = TlPT(l - p) 2 (b 0 -i) 2 ^ / +m ^min(^,m) (274) 

g =1 m =1 P 1 
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With Lq = 1 and q = j.3T , the second term (“cluster of summations”) in Eq. 
(219), Appendix A, becomes 

0T -1 2/3T—1 

2P(l-rt 2 X: E / + ”(V-.) 2 G,( 9 ," = -l,m,r = -l) (275) 

17=1 m={3T 


where 


Gi(g,i = —l,m,r = -!) = < ^ 


b - (m - (3T)], 


for m — (3T < g 


otherwise 


Factoring a T and multiplying by 1Z, we arrive at the expression for Term #3 

(3T-1 2/3T-1 r ■, 

Term #3 = RPTx 2 (l E E / + ” f X =-l,m,r =- 1 ) 

g =1 m-fiT 

(276) 

With T 0 = 1 and <7 = (3T, the third term (“cluster of summations”) in Eq. 
(219), Appendix A, yields Term #4, after factoring a T and multiplying by 7 Z. 

2/3T-12/3T-1 r 1 

Term #4 = 77 PT (1 — />) 2 ( 6 o,-i ) 2 E E / + ” 1-™™ x(g - PT,m-0T) 

g=(3T m=0T L P 1 

(277) 


The next applicable term is Eq. ( 221 ) of Appendix A. Our given parameters 
include <j>k = <l>n — 0. We can then say by recalling that b k ,t G (0, e^ k ) 


^ 71,0 — b n ,o and b k Q — 


(278) 


We will also use the facts 


27 rkl 

u k = ——- and u> n = 


2irnl 

T 


(279) 


where / is a positive integer called the normalized channel spacing integer. Now, by 
slightly reworking the u k 7 ^ cu n expression by substituting the results given in Eqs. 

(278) and (279) into the appropriate places, we get 


^1 - ^,0^,07777—771-Tvf e 


j[2irl(k - n)] 


j[2irI(k—n)—(2irkI(g)/f3T)+(2wnI(m)/0T)] 


J[(2TvI/pT)(k-n) ma.x(g,m)-(27rkI(g)/l3T) + (2TvnI(m)/pT)] 


for k 7 ^ n (280) 
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Note: The 2irl(k — n) term in the first phasor of Eq. (280) contributes integer 
multiples of 2ir to the phasor, and so it is neglected in the next rewrite of the term. 
Then, by factoring the common factor e f( 2 ’ r/ / PT)ftnm)-(kg)]^ we arr j ve 

Tb k ,ob n<0 e^ I/(ST ^ nm) - (k9)] 


<Pi = 


1 _ gi(2ir///3T)(fc—n) max(j,m) ? for k ^ U (281) 


j[2irl(k - n)] 

Finally, after we factor out a T from both the ujk ^ u n and u>k = ui n arguments to 
the summations, rearrange the summations for more logical programming, and use 
the fact that q = f3T, we arrive at the final form for Term #5. 

+M/2 +M/2 PT-10T-1 

Term #5 = TZPT( 1 - pf E E (M(M E E P 3+m x ¥>i ( 282 ) 

k=-M /2 n=-M/2 9=0 m= ° 

k^Q n^O 


where 




= < 


ej(2ir I / (3T)[(nm)—(kg)] 

j[2irl(k - ra)] 


,j[(2wkI/0T)(m-g)\ 


1 - 


^ _ gj(27r// /3T)(k—n) max(s,m) 

max(< 7 , m) 


f3T 


for k ^ n 


for k = n 


The next applicable term is the first term (“cluster of summations”) in Eq. 
(227), Appendix A. Utilizing L 0 = 1, factoring a T, multiplying by K, and noting 
that q = (3T yields 

0T—10T—1 fi ( n ; _ _i \ 

Term #6 = 7lPTx 2(1 - p) 2 (6 0 ,o)(V-i) E E P 9+m T - ~ ( 283 ) 

g- =0 m~\ 


where 


Gj{g,i = ~l,m) 
T 


= < 


1 

Jr 


[m-g], 


for g < m 


0, otherwise 

Looking at the second term (“cluster of summations”) in Eq. (227), Appendix 

A, and after utilizing L 0 = 1, factoring a T, multiplying by 71, and substituting 
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q = /3T, we have 


/3X—1 20X-1 


Term#7 = KPTx2(l-rf J (V„)(V-i) E E / +m 1 - ™ max(y, m - /3T) 

5=0 m=/3T L P 1 


(284) 


The next applicable term is the first term (“cluster of summations”) in Eq. 
(228), Appendix A. Since <f>k = 0 


^k,0 ~ ^k,0 


(285) 


We can then work with the term inside of the summations 


<j> 7 = p 9+m b 00 b k0 _ \ e ^k(T-(m/0)) _ e jw k ((max(g,m)-in)//3) 

’ ’ jUk 1 


+P 9+m KoKpy ——y [e _jw * (T ' (m//?)) - e -><c((max(5,m)-m)//3)J (286) 


Then, since b* k 0 = b k 0 and after some rearranging, we have 


P 3+m b 0 , 0 b kfi 

<P7 = - : -- X 

u k 


Juk(T-(m/P)) _ e -juk{T-{m/P)) _ e jw k ((max(g,m)-m)/P) e -jw k ((max(g,m)—m)/P) 


Then factoring a —1 from A yields 


(287) 


0 juj k ((max(g,m)—m)//3) _ - ju k ({mzx(g,m)-m)/(3) 


(288) 


Euler’s Relationship defines 


e je - e~ je 


(289) 


Applying this relationship yields 


2p 9+m b 00 b ki0 \ . ( ( m\\ . ( ( max(< 7 , m) — m 

v, =---Um ^ IT - - 1 I - sm U I- J - 


( 290 ) 




and since u k = 2irkI/T 

ip 7 =2 pS '/^j\- sin 2 * kI ( l ~Jf) _sin ^r( max (^> m ) " m ) ( 291 ) 


After factoring a T, multiplying by ft, utilizing q = /?T, and rearranging the summa¬ 
tions for more logical programming, we arrive at the expression for Term $8. 

+M/ 2 PT—1 0T -1 

Term #8 = ftPT(l — p) 2 Y S H ^7 

k=—M/2 9=0 TO= ° 


where 


p 9+m b 0fi bk,o ■ 0 , r A m 

9,7 = ,ti sm 2M r"^r 


sin (max(«, m) - m) (292) 


The next applicable term (“cluster of summations”) is Eq. (229) of Appendix 
A. Since <f> k = 0, we know &*, 0 = &£ i0 - Using L 0 = 1, we can rework the argument of 
the summations. Making the appropriate parameter substitutions yields 

¥>io = P 9+m bo,-ih,o—U j “ kl{9 - m)/(3] -l]G 6 (9,i = -l,m) 

JUk 1 J 

+/ +m V-iVop^j [e-^UU- m )/tf - l] Gete,* = -l,m) (293) 

Further symplifying yields 


P 9+m (b 0 .- l )(h,o) „ , . , n e i-*[( fl -m)//3] _ )//?) 

<^io=- G&(g,i — — l,m) - 

Wk J 


2 sin [uk[{g-m)/P]] 


Using ujk — 2irkI/T yields 


Tp 9+m (b 0 \)(b k fl) , . \2xkl 

¥>io = —- VfcJ = ~ 1 » m ) sin ~pf~ { < 9 ~ rn > 


(294) 


(295) 


Substituting <^>io inside of the summations, multiplying by ft, factoring T, using 
q = /3T, and rearranging the summations for logical programming yields the following 
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for Term #9. 


+ M/2 0T-10T-1 

Term #9 = 7ZPT(l — p) 2 £ ^ E ^ 

k=-M /2 9=1 m=0 


(296) 


where 


amd 


p 9+m (b 0 i)(bk 0 ) . 1 , . 

<^10 = -nF-= — 1 , m) sin 


nkl 


2irkl 

~w 


(.9 ~ m ) 


1, 

0, 


for m < g 
otherwise 


G 6 (g,i - -1 ,m) = 

Look at the last applicable term (“cluster of summations”), Eq. (232), Appendix 
A. Again we take advantage of the fact that <f) k = 0, which means b k fl = K 0 . Making 
the appropriate parameter substitutions yields 


¥>13 


p’ +m (K- i)(M i 


(m//?)] _ e -juik[T-(m/p)] 


Ju/ k [(l/0) ma.x(g—0T,m)—m/0] _ -jw k [(l/0)max(g-0T,m)-m/0] 


]] 


(297) 


Using Euler’s Relationship we get 

Tp3 +m (b 0 ^){b kfi ) 


¥l3 


7T kl 


“ 

( m\ 


2irkl 1 

sin 


— sin 

p T (max(j pT,m) m ) 


(298) 


Substituting 3 back inside of the summations, factoring a T, multiplying by 7Z, 
utilizing q = /3T, and rearranging the summations for logical programming yields 


+M/2 20T —1 0T — 1 

Term #10 = 1ZPT(1 — p) 2 £ E E P» 

km-U /2 ’= 0T ’"=° 


(299) 


where 


¥>13 — 


7T kl 



( m \ 


2irkl, , 

sin 

i- 

bo 

^ H- 1 

1 

: 

i _ 

— sin 

p T (ma x(g /3T,m) m) 
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For the given parameters L 0 = 1, L = 0, <j>k = 0, and Uk = ‘Inkl/T, we have 
now derived all ten programmable terms of X where 



Term #1 Term #2 Term #10 


(300) 


We will use these ten programmable terms to compute probability of bit error for the 


dense WDM system. 
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APPENDIX C 


PROGRAMMING STRATEGY AND COMPUTER 
PROGRAMS FOR GENERATION OF PROBABILITY 
OF BIT ERROR GRAPHS 


In this appendix, we present the computer programs written to generate the proba¬ 
bility of bit error graphs for the system under consideration. We attempt to generate 
four graphs, one for each value of /3T considered. We reiterate that (3 is the free 
spectral range of the Fabry-Perot filter (Hz) and that T is the data bit period (s). 
PT is called the free spectral range-bit period product. The four values of PT used 
are 

PT = [500 1000 1500 2000] 


In each of the four graphs there will be five traces as we will present probability bit 
error for four values of the normalized channel spacing integer /, or equivalently the 
number of adjacent channels M, along with a probability of bit error trace for single 
channel (SC) operation without Fabry-Perot (FP) filtering or SC operation with FP 
filtering and without ISI or ACI. We will show the relationship between I and the 
number of adjacent channels M later. Now, however, we present the four values of I 
corresponding to each value of PT: 


For PT — 500 


I — [A 5 8 20] 


For PT = 1000 


I = [5 6 9 20] 


For pT = 1500 


/ = [7 9 12 20] 
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For /3T = 2000 


/ = [8 9 12 20] 


Without a lengthy discourse, we present the mathematical relationship between 
the normalized channel spacing integer / and the number of adjacent channels M 

1 < 301 ) 

For the special case /*. = k Af = kl/T we can easily see that A/ = I/T. Then 

3T 

M = -j - 1 (302) 


For the values of I above M will not always be an even integer. Thus, an algorithm 
to consistently arrive at an even integer value of M which is less than or equal to the 
true mathematical value presented in Eq. (302) was devised. To get the number of 
adjacent channels, we perform the following operation 


f3T 


— Q + R 


(303) 


where Q is the integer quotient of the division operation and R is the remainder. To 
arrive at M we 


1. Subtract 1 from Q if Q is an odd integer and R = 0. 

2. Subtract 2 from Q if Q is an even integer and R = 0. 

3. Subtract 1 from Q if Q is an odd integer and R ^ 0. 

4. Subtract 2 from Q if Q is an even integer and R ^ 0. 

Using these rules, we obtain at the four values of M for each value of (3T. 
For (3T = 500 

M = [124 98, ^ JMJ 

1=4 1 =5 /=8 1 =20 
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For f3T = 1000 

For f3T = 1500 

For (3T = 2000 


M = [198 164 U0 ,48,] 

7=5 7=6 7=9 7=20 


M = [212 164 124 

7=7 7=9 7=12 7=20 


M = [248 220 164 ,98, ] 

7=8 7=9 7=12 7=20 

Now we use these four values of I and four values of M for each value of (3T 
to compute the four graphs of probability of bit error with five traces each. We will 
use the equations and methods developed in Appendix B to compute these graphs. 
Recall that 


x = n pr(QUOT x ) + n pt(quot 2 ) + • • • + n pt{ quot 10 ) ( 304 ) 

'- v .- ' '- V ' '-V-" 

Term #1 Term #2 Term #10 

Realizing this, we program the ten terms given in Appendix B in the following way: 

TERM #1 [Eq. (272)]: We will compute the QUOT x portion of this term four times, 
once for each value of fiT = 500, 1000, 1500, and 2000. We will use each 
of these four values in a separate program to compute the probability of 
bit error according to the equations developed in Appendix B. 

TERM #2 [Eq. (274)]: We will compute the QUOT 2 portion of TERM #2 four times, 
once for each value of f3T = 500, 1000, 1500, and 2000. Each of these 
values will be used in a separate program to compute the probability of 
bit error according to the equations developed in Appendix B. 

TERM #3 [Eq. (276)]: We will compute the QUOT 3 portion of TERM #3 four 
times, once each value of (3T = 500, 1000, 1500, and 2000. Each of these 
values will be used in a separate program to compute the probability of 
bit error according to the equations developed in Appendix B. 
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TERM #4 [Eq. (277)]: We will compute the QUOT 4 portion of TERM #4 four times, 
once for each value of (3T = 500, 1000, 1500, and 2000. Each of these 
values will be used in a separate program to compute the probability of 
bit error according to the equations developed in Appendix B. 

TERM #5 [Eq. (282)]: We will compute the value of QUOT s a total of 48 times, 12 
times for each value of /3T = 500, 1000, 1500, and 2000. For each value 
of (3T, we have the four corresponding values of the normalized channel 
spacing integer /, or equivalently, the number of adjacent channels M. 
We will compute the value of QUOT 5 three times for each value of I. 
QUOTs will be computed once for all M/2 of the lower adjacent channels 
being packed with Is and the upper adjacent channels being packed with 
0s. This is the case: ( b~ — l,b + = 0). QUOT 5 will be computed once for 
the lower adjacent channels being packed with 0s and the upper adjacent 
channels being packed with Is. This is the case: ( b~ = 0, b + = 1). Finally, 
we will compute QUOT 5 once more for both the lower and upper adjacent 
channels being packed with l’s. This is the case: {b~ = l,b + = 1). These 
twelve values of QUOT 5 for each value of j3T will be used in each of the 
four separate bit error programs to compute each of the four multiple 
channel probability of bit error traces using the equations and methods 
developed in Appendix B. 

TERM #6 [Eq. (283)]: We will compute the QUOT 6 portion of TERM #6 four times, 
once each for (3T = 500, 1000, 1500, and 2000. Each of the values will be 
used in separate to compute program probability of bit error according 
to the equations developed in Appendix B. 
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TERM #7 [Eq. (284)]: We will compute the QUOT 7 portion of TERM #7 four times, 
once for each value of /3T = 500, 1000, 1500, and 2000. Each of these 
values will be used in a separate program to compute the probability of 
bit error according to the equations developed in Appendix B. 

TERM #8 [Eq. (292)]: We compute the value of QUOT 8 twelve times for each 
value of /3T. These twelve values of QUOT 8 for each value of /?T will be 
used in each of the four separate bit error graphing programs to compute 
each of the four multiple channel probability of bit error traces using the 
equations and methods developed in Appendix B. 

TERM #9 [Eq. (296)]: We compute the value of QUOT 9 twelve times for each 
value of f3T. These twelve values of QUOT 9 for each value of /3T will be 
used in each of the four separate bit error graphing programs to compute 
each of the four multiple channel probability of bit error traces using the 
equations and methods developed in Appendix B. 

TERM #10 [Eq. (299)]: We compute the value of QUOT 10 twelve times for each 
value of /3T. These twelve values of QUOT 10 for each value of (3T will be 
used in each of the four separate bit error graphing programs to compute 
each of the four multiple channel probability of bit error traces using the 
equations and methods developed in Appendix B. 

Note: When we say we are computing QUOT^, we are not being exactly mathemat¬ 
ically correct, as we compute each of these terms without the bit values 6 0 ,o and/or 

6 0 ,-i the factor (1 — p) 2 . These are accounted for in the final programs. For example 

Term #1 = K PT (1 - pf b\ 0 £ £ (l - ) (305) 

'- - -' 

QUOTj 
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However, the computer program for Term (QUOTj) only computes 


s!> (-hr)) 


(306) 


four times for (3T = 500, 1000, 1500, 2000. (1 - p ) 2 and the value of (6 0 , 0 ) 2 are 
accounted for in the final programs which utilize the equations and methods developed 
in Appendix B. The author apologizes for the slight stretch of the truth, but it seemed 
necessary to succinctly explain the general method of computing each of the ten terms, 
and the final graphs. 

For completness we also present the probability of bit error equation for a single 
channel operation 

p b = Q (^) (307) 

where we recall [Eq. (271), Appendix B] that 


Z 


TZPT 

VNoT 


= np 



(308) 


and 

Q(x) = jT e- y2 / 2 dy (309) 

The computer programs for each of the ten terms, the numerical results, and the final 
graph are now presented. 
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xmxmmmmxmmxmmmmxmmxxmxmmmmmmxmxm 


l THESIS COMPUTER WORK 

7. _ 

7. I I 

% | TERM #1 | 

7. I_I 

7 . 

7. COMMENTS: THE IDEA OF THIS PROGRAM IS TO COMPUTE THE FOUR VALUES OF 
TERM #1 AND TO WRITE THESE VALUES TO DIARY FILE FOR LATER 
USE/MANIPULATION IN CALCULATING THE DETECTION STATISTIC 
AND THE PROBABILITIES OF BIT ERROR 


7. 

7. 

7. 

7. 

7. 

7, JOHN A. STUDER 
7. CPT, U.S ARMY 
7. 550-53-7181 

mmxxnmmmxmmxmxmmmmmmmxmmmmmmxmx 


DATE LAST MODIFIED: 11 SEP 94 


rho = 0.99; 

betatau = [500 1000 1500 2000]; 
for i = 1:4 

termlsum(i)= 0; 
for g = l:betatau(i) 
for m = l:betatau(i) 

termlsum(i) = termlsum(i) + (rh.o.' ((g-l) + (m-l))*. . . 

(l-(max([g—1 m-1]’)/betatau(i)))) ; 

end 

end 

end 

diary johnmanl.txt 

diary on 

termlsum 

diary off 

end 
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termlsum = 


1.0e+03 * 


7.0510 

8.5126 

9.0083 

9.2563 

’/, numbers/arrows/words 
*/, below added after matlab 





*/, dumped answers to file using 
'/, the diary command 

500 

1000 

1500 

2000 

- betatau values 


‘/, The numbers displayed above are the computed values of 
'/, Term #1 for the four values of betatau given above. 

‘/, These values will be used later in other programs to compute detection 
'/, statistics and the probabilities of bit error for various signal to noise 
'/, ratios. 


I TERM #1 I 
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% THESIS COMPUTER WORK 

7. 

7. 

7. 

7 . 

7. 

7. 

7, COMMENTS: THE IDEA OF THIS PROGRAM IS TO COMPUTE THE FOUR VALUES 

7. 

7. 

7. 

7. 

7. 

7. 

7. 

7. 

7. 

7. JOHN A. STUDER 
7. CPT, U.S. ARMY 
7. 550-53-7181 


I I 

| TERM #2 | 

I_I 


OF TERM #2 CORRESPONDING TO THE FOUR VALUES OF 
betatau: 500,1000,1500,2000. WE WILL USE BACKGROUND 
PROCESSING TO WRITE THESE FOUR VALUES TO AN OUTPUT 
FILE CALLED term2.out. WE WILL USE THESE FOUR OUTPUT 
VALUES FOR TERM #2 TO LATER, IN ANOTHER PROGRAM 
COMPUTE THE DETECTION STATISTICS AND PROBABILITIES 
OF BIT ERROR FOR THE VARIOUS SIGNAL TO NOISE RATIOS. 


DATE LAST MODIFIED: 11 SEP 94 


rho = 0.99; 

betatau = [500 1000 1500 2000]; 
for i = 1:4 

term2sum(i) = 0; 
for g = 1:betatau(i)-l 

for m = 1: betatau(i)-l 

term2sum(i) = term2sum(i) + ((rho'(g+m))*min(g,m)); 

end 

end 

terra2sumfinal(i) = term2sum(i)/betatau(i); 

end 

term2sumfinal 
exit 
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< M A T L A B (R) > 

(c) Copyright 1984-94 The MathWorks, Inc. 

All Rights Reserved 
Version 4.2 
Mar 29 1994 

Commands to get started: intro, demo, help help 

Commands for more information: help, whatsnew, info, subscribe 


term2sumfinal = '/, The numbers/arrows/words 

7. below were added after 

959.5662 492.4271 328.3413 246.2563 '/, the matlab background job 


*/, dumped the values to the 

7. " ~ “ ’/, file term2.out 

7. 500 1000 1500 2000 '/, All values/text not done by 

7. ~ ’/, Sun Stn is preceded by a "7," 

7, -betatau values- 



29965020 flops. 





THESIS COMPUTER WORK 



mmmmmmmn 


7. 

7 . 

7, 

7. 

7. 

7. 

7. 

7. 

7. 

7. 

7 . 

7 . 

7. 

7. 

7. 

7. 

7. 

7. 

7. 

7. 


| TERM #3 I 


COMMENTS: THE IDEA OF THIS PROGRAM IS TO COMPUTE THE FOUR VALUES OF 
TERM #3 CORRESPONDING TO THE FOUR VALUES OF 
betatau: 500,1000,1500,2000. WE WILL USE BACKGROUND 
PROCESSING TO COMPUTE THE VALUES AND WRITE THESE FOUR 
VALUES TO A FILE CALLED term3.out. WE WILL USE THESE 
FOUR VALUES FOR TERM #3 TO LATER, IN ANOTHER PROGRAM 
COMPUTE THE DETECTION STATISTICS AND PROBABILITIES OF 
BIT ERROR FOR THE VARIOUS SIGNAL TO NOISE RATIOS. 


JOHN A. STUDER 
CPT, U.S. ARMY 
550-53-7181 


DATE LAST MODIFIED: 12 SEP 1994 


format short e 
rho = 0.99; 

betatau = [500 1000 1500 2000]; 
for i = 1:4 

term3sum(i) = 0; 
for g = 1:betatau(i)-l 

for m = betatau(i):((2*betatau(i))-l) 
if m-betatau(i) < g 

term3sum(i) = term3sum(i) + ((rho~(g+m))*(g-(m-betatau(i)))); 

end 


end 

end 

term3sumfinal(i) = (2*term3sum(i))/betatau(i); 
end 

term3sumfinal 
exit 
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< M A T L A B (R) > 

(c) Copyright 1984-93 The MathWorks, Inc. 

All Rights Reserved 
Version 4.1 
Jun 15 1993 

Commands to get started: intro, demo, help help 

Commands for more information: help, whatsnew, info, subscribe 


term3sumfinal = '/, Numbers/arrows/words 

7, below were added after 

1.2211e+01 4.2917e-02 1.8815e-04 9.2720e-07 */, the matlab background 

7. job dumped the values 

*/. ~ " '/. to the file term3.out 

•/. 500 1000 1500 2000 7, All values/text not 

7. 7. done by Sun Stn. is 

7. '/, preceded by a 

'/, -betatau values- 


mmm mm 

7. 7. 

7. TERM #3 7. 

7. 7. 


29990004 flops. 
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xxxxxxxxxxxxxx xxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxx 

THESIS COMPUTER WORK 


X 

% 


X 

X 

X 

X 

•/. 

X 

X 

X 

•/. 

X 

X 

X 

X 


I TERM #4 


COMMENTS: THE IDEA OF THIS PROGRAM IS TO COMPUTE THE FOUR VALUES 
OF TERM #4 CORRESPONDING TO THE FOUR VALUES OF 
betatau: 500,1000,1500,2000. WE WILL USE BACKGROUND 
PROCESSING TO WRITE THESE FOUR VALUES TO AN OUTPUT 
FILE CALLED term4.out. WE WILL USE THESE FOUR OUTPUT 
VALUES FOR TERM #4 TO LATER , IN ANOTHER PROGRAM 
COMPUTE THE DETECTION STATISTICS AND PROBABILITIES 
OF BIT ERROR FOR THE VARIOUS SIGNAL TO NOISE RATIOS. 


X 

X 

X 

X 


JOHN A. STUDER DATE LAST MODIFIED: 13 SEP 94 

CPT, U.S. ARMY 

550-53-7181 


format short e 
rho = 0.99; 

betatau = [500 1000 1500 2000]; 
for i = 1:4 

term4sum(i) =0; 

for g = betatau(i):((2*betatau(i))-l) 
for m = betatau(i):((2*betatau(i))-l) 

term4sum(i) = term4sum(i) + ((rho~(g+m))*... 

(1-(max(g-betatau(i),m-betatau(i))/betatau(i)))); 

end 

end 

term4sumfinal(i) = term4sum(i); 

end 

term4sumfinal 
exit 
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< M A T L A B (R) > 

(c) Copyright 1984-94 The MathWorks, Inc. 

All Rights Reserved 
Version 4.2 
Mar 29 1994 

Commands to get started: intro, demo, help help 

Commands for more information: help, whatsnew, info, subscribe 


term4sumfinal = 

3.0440e-01 

1.5865e-05 

7.2482e-10 

3.2152e-14 

'/, Numbers/arrows/words 
'/, below were added after 
'/, the matlab bkgd job 

X 


- 

- 

X dumped the values to 

X the file term4.out 

X 500 

1000 

1500 

2000 

# / t All values/text not 

X 

X 

X - 

- betatau 

values - 


7, done by Sun Stn is 
'/, preceded by a 


xxxxxxxxxxxxxxxx 

y. */. 

X TERM #4 */. 

’/. */. 

immtmmn 


60010008 flops. 
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THESIS COMPUTER WORK 


TERM #5 
betatau = 500 


COMMENTS: THE IDEA OF THIS PROGRAM IS TO COMPUTE TERM #5 
FOR THE VALUE OF betatau = 500. WE WILL COMPUTE 
FOR THE FOUR VALUES OF I (AN INTEGER). THE VALUES OF I FOR 
betatau = 500 are I = [4 5 8 20]. WE WILL COMPUTE THE VALUE OF 
TERM #5 THREE TIMES FOR EACH CORRESPONDING VALUE OF I. ONCE FOR 
ALL M/2 OF THE LOWER ACI CHANNELS (have negative summation index) 
BEING PACKED WITH l'S AND THE UPPER ACI CHANNELS ALL BEING PACKED 
WITH O'S. THIS IS THE CASE: (b-=l,b+=0). ONCE FOR ALL THE LOWER 
ACI CHANNELS BEING PACKED WITH ZERO'S AND THE UPPER ACI CHANNELS 
(have a positive summation index) BEING PACKED WITH l’S. THIS 
IS THE CASE: (b-=0,b+=l). FINALLY ONCE FOR BOTH THE UPPER AND 
LOWER ACI CHANNELS BEING PACKED WITH l'S. THIS IS THE CASE: 

(b-=l,b+=l). WE WILL USE BACKGROUND PROCESSING WHICH WILL 
DUMP ALL TWELVE VALUES TO A FILE CALLED term5500.mat 
SINCE THE FILE IS IN BINARY FORMAT WE WILL USE MATLAB 
INTERACTIVE COMMANDS AND THE TEXT EDITOR TO PRODUCE 
READABLE HARDCOPY RESULTS. AFTER THIS IS ALL DONE, WE WILL 
USE THESE VALUES LATER , IN ANOTHER PROGRAM TO COMPUTE THE 
DETECTION STATISTICS AND PROBABILITIES OF BIT ERROR FOR THE 
VARIOUS SIGNAL TO NOISE RATIOS. 


JOHN A. STUDER 
CPT, U.S. ARMY 
550-53-7181 


DATE LAST MODIFIED: 22 OCT 94 




format short e 
rho = 0.99; 
betatau = 500; 

I = [4 5 8 20] ; 

M = [124 98 60 24]; 

g = 0:(betatau-1); 

m = 0:(betatau-1); 

row = rho.~g; 

col = rho.~m; 

rhomatrix = col'*row; 

rowprime = g(ones(betatau,1),:); 

colprime = rowprime'; 

WO = max(rowprime,colprime); 

C = WO/betatau; 

ONEMINUSMAX = 1-C; 


135 










A = colprime-rowprime; ’/, m-g for later complex exponential 

for i = 1:4 

bitmatrix = [-M(i)/2:-l zeros(l,M(i)/2);zeros(l,M(i)/2) l:M(i)/2;... 

-M(i)/2:-1 1:M(i)/2]; 

const = ((2*pi*I(i))/betatau); 
for bitpat = 1:3 

term5sum(i,bitpat) = 0; 
work = bitmatrix(bitpat,:); 
m = work==0; 
work(m) = [] ; 
k = work; 
n = k; 

for kct = 1:length(k) 
for net = 1:length(n) 
if kct == net 

B = j*((2*pi*k(kct)*I(i)*A)/betatau); 

D = exp(B); 

term5sum(i,bitpat) = term5sum(i,bitpat)+... 

sum(sum(rhomatrix.*D.*0NEMINUSMAX)); 


else 

El = exp(j*const*((n(nct)*colprime)-(k(kct)*rowprime))) 

E2 = l-exp(j*const*(k(kct)-n(nct))*WO); 
term5sum(i,bitpat) = term5sum(i,bitpat) +... 
sum(sum((rhomatrix.*E1.*E2)/(j *2*pi*I(i)*(k(kct)-n(nct ))))); 

end 


end 


end 

end 


end 

save term5500 term5sum 
exit 
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< M A T L A B (R) > 

(c) Copyright 1984-93 The MathWorks, Inc. 
All Rights Reserved 
Version 4.1 
Jun 15 1993 


Commands to get started: intro, demo, help help 

Commands for more information: help, whatsnew, info, subscribe 

>> load term5500.mat 
>> who 

Your variables are: 
term5sum 
>> term5sum 
term5sum = 


1.0e+03 * 


1.5423 - O.OOOOi 
0.9419 - O.OOOOi 
0.3307 - O.OOOOi 
0.0431 - O.OOOOi 


1.5423 + O.OOOOi 
0.9419 + O.OOOOi 
0.3307 + O.OOOOi 
0.0431 + O.OOOOi 


1.3825 + O.OOOOi ’/, I = 4 

0.9001 - O.OOOOi % I = 5 

0.3576 + O.OOOOi '/. I = 8 

0.0580 - O.OOOOi 7. I = 20 


» 

'/. CASE: (b-=l,b+=0) CASE: (b-=0,b+=l) 


CASE: (b-=l,b+=l) 


mmmmmmx 

•/. TERM #5 '/. 

'/, betatau = 500 */, 

mXmmmmm 


NOTE: THIS IS THE FILE term5500.out CREATED FROM THE 

BINARY FILE term5500.mat USING MATLAB INTERACTIVE COMMANDS. 
ALL TEXT PRECEDED BY A '/. WAS ADDED LATER USING A 
TEXT EDITOR. 
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X THESIS COMPUTER WORK 

X 

X 

X _ 

X I "l 

X I TERM #6 | 

X I_I 

X 

X COMMENTS: THE IDEA OF THIS PROGRAM IS TO COMPUTE THE FOUR VALUES OF 
X TERM #6 CORRESPONDING TO THE FOUR VALUES OF betatau WHICH ARE: 

X 500,1000,1500,2000. WE WILL USE BACKGROUND PROCESSING TO 

X COMPUTE THE VALUES AND WRITE THESE FOUR VALUES TO A FILE 

X CALLED term6.out. WE WILL USE THESE FOUR VALUES FOR TERM #6 

X TO LATER, IN ANOTHER PROGRAM COMPUTE THE DETECTION STATISTICS 

X AND PROBABILITIES OF BIT ERROR FOR THE VARIOUS SIGNAL TO 

X NOISE RATIOS. 

X 

X JOHN A. STUDER DATE LAST MODIFIED: 16 SEP 94 

X CPT, U.S. ARMY 
X 550-53-7181 


format short e 
rho = 0.99; 

betatau = [500 1000 1500 2000]; 
for i = 1:4 

term6sum(i) = 0; 
for g = l:betatau(i) 

for m = 1:betatau(i)-l 
if (g-1) < m 

term6sum(i) = term6sum(i) + ((rho~((g-l)+m))*(m-(g-l))); 

end 


end 

end 

term6sumfinal(i) = (2*term6sum(i))/betatau(i); 
end 

term6sumfinal 
exit 
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< M A T L A B (R) > 

(c) Copyright 1984-94 The MathWorks, Inc. 

All Rights Reserved 
Version 4.2 
Mar 29 1994 

Commands to get started: intro, demo, help help 

Commands for more information: help, whatsnew, info, subscribe 


term6sumfinal = 

1.8585e+03 

9.9411e+02 

6.6331e+02 

4.9749e+02 

'/, Numbers/arrows/words 
'/, below were added after 
*/, the mat lab bkgd job 

7. 


- 

- 

7, dumped the values to 
y, the file termS.out 

7, 500 

1000 

1500 

2000 

'/, All values/text not 

*/. 

y. - 

-betatau 

values - 


'/, done by Sun Stn is 
'/, preceded by a . 


•/. •/. 

*/. TERM #6 */. 

•/. */. 

mmmmm 


33732508 flops. 
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'/. THESIS COMPUTER WORK 

•/. 

•/, 

*/. _ 

7. I I 

7 . I TERM #7 I 

•/. I_I 

•/. 

•/. COMMENTS: THE IDEA OF THIS PROGRAM IS TO COMPUTE THE FOUR VALUES OF 
7 . TERM #7 CORRESPONDING TO THE FOUR VALUES OF betatau WHICH ARE 

’/„ 500,1000,1500,2000. WE WILL USE BACKGROUND PROCESSING TO COMPUTE 

7 . THE VALUES AND WRITE THESE FOUR VALUES TO A FILE CALLED 

7 . t erm7. out. WE WILL USE THESE FOUR VALUES FOR TERM #7 TO 

*/. LATER, IN ANOTHER PROGRAM COMPUTE THE DETECTION STATISTICS 

7, AND PROBABILITIES OF BIT ERROR FOR THE VARIOUS SIGNAL TO 

'/. NOISE RATIOS. 

7 . 

7 ; JOHN A. STUDER DATE LAST MODIFIED: 16 SEP 94 

•/. CPT, U.S. ARMY 
7 . 550-532-7181 


format short e 
rho = 0.99; 

betatau = [500 1000 1500 2000]; 
for i = 1:4 

term7sum(i) = 0; 
for g = l:betatau(i) 

for m = betatau(i):((2*betatau(i))-l) 
term7sum(i) = term7sura(i) + ... 

((rho~((g-l)+m))*(1-(max(g-l,m-betatau(i))/betatau(i)))); 

end 

end 

term7sumfinal(i) = 2*term7sum(i); 
end 

term7sumfinal 
exit 
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< M A T L A B (R) > 

(c) Copyright 1984-94 The MathWorks, Inc. 

All Rights Reserved 
Version 4.2 
Mar 29 1994 

Commands to get started: intro, demo, help help 

Commands for more information: help, whatsnew, info, subscribe 


term7sumfinal 

9.2657e+01 

7.3500e-01 

5.1105e-03 

3.4503e-05 

7 . Numbers/arrows/words 

7. below were added after 

7, the mat lab bkgd job 

7 . ~ 

- 

— 

- 

7. dumped the values to 

7 . the file term7.out 

7 . 500 

1000 

1500 

2000 

7 . All values/text not 

7 . 

7 . 

7 . - 

- betatau 

values 


7. done by Sun Stn is 

7 . preceded by a " 7 .". 


7.77.7.7.7.7.7.7.7.7.7.7.7.7. 

•/. •/. 


7 TERM #7 7 . 
7 . 7 . 


7 . 7 . 7 . 7 . 7 . 7 . 7 . 77 . 7 . 7 . 7 . 7 . 7 . 7 . 


67510004 flops. 
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’/. THESIS COMPUTER WORK 

% 

% 

y. 

y. _ 

y. i i 

'/. I TERM #8 I 

'/, I betatau = 500 | 

•/. i_i 

7 . 

7 . COMMENTS: THE IDEA OF THIS PROGRAM IS TO COMPUTE TERM #8 
•/. FOR THE VALUE OF betatau = 500. WE WILL COMPUTE 

'/. FOR THE FOUR VALUES OF I (AN INTEGER) . THE VALUES OF I FOR 

# /. betatau = 500 are I = [4 5 8 20] . WE WILL COMPUTE THE VALUE OF 

*/, TERM #8 THREE TIMES FOR EACH CORRESPONDING VALUE OF I. ONCE FOR 

'/. ALL M/2 OF THE LOWER ACI CHANNELS (have negative summation index) 

’/. BEING PACKED WITH l'S AND THE UPPER ACI CHANNELS ALL BEING PACKED 

y. WITH O'S. THIS IS THE CASE: (b-=l,b+=0). ONCE FOR ALL THE LOWER 

•/. ACI CHANNELS BEING PACKED WITH ZERO'S AND THE UPPER ACI CHANNELS 

'/, (have a positive summation index) BEING PACKED WITH l'S. THIS 

*/. IS THE CASE: (b-=0,b+=l). FINALLY ONCE FOR BOTH THE UPPER AND 

% LOWER ACI CHANNELS BEING PACKED WITH l'S. THIS IS THE CASE: 

*/. (b-=l,b+=l) . WE WILL USE BACKGROUND PROCESSING WHICH WILL 

*/. DUMP ALL TWELVE VALUES TO A FILE CALLED term8500. out. WE WILL 

'/. USE THESE VALUES LATER , IN ANOTHER PROGRAM TO COMPUTE THE 

% DETECTION STATISTICS AND PROBABILITIES OF BIT ERROR FOR THE 

y. VARIOUS SIGNAL TO NOISE RATIOS. 

% 

7 . JOHN A. STUDER DATE LAST MODIFIED: 26 SEP 94 

'/. CPT, U.S. ARMY 
’/. 550-53-7181 


format short e 
rho = 0.99; 
betatau = 500; 

I = [4 5 8 20] ; 

M = [124 98 60 24]; 
for i = 1:4 

k = [-M(i)/2:-l l:M(i)/2]; 
term8sum500(i) = 0; 
for ct = 1:M(i) 

for g = 1:betatau 
for m = 1:betatau 

term8sum500(i) = term8sum500(i) + ((rho"((g-l)+(m-l)))/... 

(pi*k(ct)*I(i)))*(sin(2*pi*k(ct)*I(i)*... 
(l-((m-l)/betatau)))-sin(((2*pi*k(ct)*I(i))/... 
bet atau)*(max(g-1,m-1)-(m-1)))); 

end 

end 

if ct == M(i)/2 
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term85001eft(i) = term8sum500(i); 

end 


end 

term85001r(i) = term8sum500(i); 

term8500right(i) = term85001r(i)-term85001eft(i) 


end 

term85001eft 

term8500right 

term85001r 

exit 
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< M A T L A B (R) > 

(c) Copyright 1984-94 The MathWorks, Inc. 

All Rights Reserved 
Version 4.2 
Mar 29 1994 

Commands to get started: intro, demo, help help 

Commands for more information: help, whatsnew, info, subscribe 


term85001eft = 

-3.7098e+02 -2.3852e+02 -9.2820e+01 -1.4109e+01 1 CASE: (b-=l,b+=0) 
term8500right = 

-3.7098e+02 -2.3852e+02 -9.2820e+01 -1.4109e+01 % CASE: (b-=0,b+=l) 


term85001r = 

-7.4196e+02 -4.7703e+02 -1.8564e+02 -2.8217e+01 '/. CASE: (b-=l,b+=l) 

*/. 

'/. 1 = 4 1 = 5 1 = 8 1 = 20 

*/. - 

•/. 

*/. VALUES OF NORMALIZED CHANNEL SPACING INTEGER I 


y.y.y.y.y.y.y.y.y.y.y.y.y.y.y.y.y.y.y.y. 

y. TERM #8 y. 

'/, betatau = 500 '/, 


NOTE: ALL TEXT PRECEDED BY A */. WAS ADDED 

BY A TEXT EDITOR AFTER THE BACKGROUND 
JOB CREATED THE FILE term8500.out. 


2142000318 flops. 
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THESIS COMPUTER WORK 


*/. 

7 . 

7 . 

•/. 

7 . 

*/. 

•/. 

*/. 

•/. 

*/. 

y. 

•/. 

y. 

y. 

y. 

•/. 

•/. 

•/. 

y. 

•/. 

•/. 

y. 

y. 

y. 

y. 

y. 

y. 

y. 

y. 

y. 


TERM #8 

betatau = 1000 


I. 


.1 


COMMENTS: THE IDEA OF THIS PROGRAM IS TO COMPUTE TERM #8 

FOR THE VALUE OF betatau = 1000. WE WILL COMPUTE 
FOR THE FOUR VALUES OF I (AN INTEGER). THE VALUES OF I FOR 
betatau = 1000 are I = [5 6 9 20]. WE WILL COMPUTE THE VALUE OF 
TERM #8 THREE TIMES FOR EACH CORRESPONDING VALUE OF I. ONCE FOR 
ALL M/2 OF THE LOWER ACI CHANNELS (have negative summation index) 
BEING PACKED WITH l'S AND THE UPPER ACI CHANNELS ALL BEING PACKED 
WITH O’S. THIS IS THE CASE: (b-=l,b+=0). ONCE FOR ALL THE LOWER 
ACI CHANNELS BEING PACKED WITH ZERO’S AND THE UPPER ACI CHANNELS 
(have a positive summation index) BEING PACKED WITH l’S. THIS 
IS THE CASE: (b-=0,b+=l). FINALLY ONCE FOR BOTH THE UPPER AND 
LOWER ACI CHANNELS BEING PACKED WITH l’S. THIS IS THE CASE: 

(b-=l,b+=l). WE WILL USE BACKGROUND PROCESSING WHICH WILL 
DUMP ALL TWELVE VALUES TO A FILE CALLED term81k.mat 
SINCE THE FILE IS IN A BINARY FORMAT WE WILL THEN USE 
MATLAB INTERACTIVE COMMANDS AND THE TEXT EDITOR TO PRODUCE 
READABLE HARDCOPY RESULTS. AFTER THIS IS DONE, WE WILL 
USE THESE VALUES LATER , IN ANOTHER PROGRAM TO COMPUTE THE 
DETECTION STATISTICS AND PROBABILITIES OF BIT ERROR FOR THE 
VARIOUS SIGNAL TO NOISE RATIOS. 


JOHN A. STUDER 

CPT, U.S. ARMY 

550-53-7181 
/ */ •/ •/«/ •/«/»/ 0 / •/«/•/ •/ •/ 


DATE LAST MODIFIED: 27 SEP 94 


format short e 
rho = 0.99; 
betatau = 1000; 

I = [5 6 9 20] ; 

M = [198 164 110 48] ; 
for i = 1:4 

k = [-M(i)/2:-l l:M(i)/2]; 
term8sum(i) = 0; 
for ct = l:M(i) 

for g = 1:betatau 
for m = l:betatau 

term8sum(i) = term8sum(i) + ((rho~((g-l)+(m-l)))/... 

(pi*k(ct)*I(i)))*(sin(2*pi*k(ct)*I(i)*... 
(l-((m-l)/betatau)))-sin(((2*pi*k(ct)*I(i))/... 
betatau)*(max(g-1,m-l)-(m-l)))); 












end 


end 

if ct == M(i)/2 

term81kleft(i) = term8sum(i); 

end 

end 

term81klr(i) = term8sum(i); 

term81kright(i) = term81klr(i)-term81kleft(i); 

end 

save term81k term81kleft term81kright term81klr 
exit 
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Inc. 


< M A T L A B (R) > 

(c) Copyright 1984-93 The MathWorks, 

All Rights Reserved 
Version 4.1 
Jun 15 1993 

Commands to get started: intro, demo, help help 
Commands for more information: help, whatsnew, info, subscribe 

>> load term81k 
>> who 

Your variables are: 

term81kleft term81klr term81kright 

>> term81kleft 
term81kleft = 

-464.3344 -327.7899 -148.2912 -29.7545 '/. CASE: (b-=l,b+=0) 

>> term81kright 
term81kright = 

-464.3344 -327.7899 -148.2912 -29.7545 ’/. CASE: (b-=0,b+=l) 

>> term81klr 
term81klr = 

-928.6688 -655.5798 -296.5824 -59.5090 '/. CASE: (b-=l,b+=l) 

» 

7 . 

7 . 1 = 5 1 = 6 1 = 9 1 = 20 

*/. _ 

•/. 

'/. VALUES OF THE NORMALIZED CHANNEL SPACING INTEGER I 


*/. TERM #8 ’/. 

'/, betatau = 1000 */, 

mmxxmmmx 


NOTE: THIS IS THE FILE term81k.out CREATED FROM THE 

BINARY FILE term81k.mat USING MATLAB INTERACTIVE 
COMMANDS. ALL TEXT PRECEDED BY A 7 . WAS ADDED LATER 
USING A TEXT EDITOR. 
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*/. THESIS COMPUTER WORK 

7 . 

7 . 

7 . 

7 . _ 

7 . I I 

’/. I TERM #8 I 

'/, I betatau = 1500 I 

*/. I_I 

7 . 

7 . COMMENTS: THE IDEA OF THIS PROGRAM IS TO COMPUTE TERM #8 

7 . FOR THE VALUE OF betatau = 1500. WE WILL COMPUTE 

7 , FOR THE FOUR VALUES OF I (AN INTEGER) . THE VALUES OF I FOR 

V. betatau = 1500 are I = [7 9 12 20] . WE WILL COMPUTE THE VALUE OF 

*/. TERM #8 THREE TIMES FOR EACH CORRESPONDING VALUE OF I. ONCE FOR 

•/. ALL M/2 OF THE LOWER ACI CHANNELS (have negative summation index) 

*/. BEING PACKED WITH l'S AND THE UPPER ACI CHANNELS ALL BEING PACKED 

*/. WITH O'S. THIS IS THE CASE: (b-=l,b+=0). ONCE FOR ALL THE LOWER 

'/. ACI CHANNELS BEING PACKED WITH ZERO'S AND THE UPPER ACI CHANNELS 

*/, (have a positive summation index) BEING PACKED WITH l'S. THIS 

’/. IS THE CASE: (b-=0,b+=l). FINALLY ONCE FOR BOTH THE UPPER AND 

*/. LOWER ACI CHANNELS BEING PACKED WITH l'S. THIS IS THE CASE: 

*/. (b-=l,b+=l) . WE WILL USE BACKGROUND PROCESSING WHICH WILL 

*/. DUMP ALL TWELVE VALUES TO A FILE CALLED term815k.out WE WILL 

*/. USE THESE VALUES LATER, IN ANOTHER PROGRAM TO COMPUTE THE 

*/. DETECTION STATISTICS AND PROBABILITIES OF BIT ERROR FOR THE 

*/. VARIOUS SIGNAL TO NOISE RATIOS. 

y. 

y. JOHN A. STUDER DATE LAST MODIFIED: 27 SEP 94 

•/. CPT, U.S. ARMY 
'/. 550-53-7181 


format short e 
rho = 0.99; 
betatau = 1500; 

I = [7 9 12 20]; 

M = [212 164 124 74]; 
for i = 1:4 

k = [-M(i)/2:-l l:M(i)/2]; 
term8sum(i) = 0; 
for ct = l:M(i) 

for g = 1:betatau 
for m = 1:betatau 

term8sum(i) = term8sum(i) + ((rho“((g-l)+(m-l)))/... 

(pi*k(ct)*I(i)))*(sin(2*pi*k(ct)*I(i)*... 
(l-((m-l)/betatau)))-sin(((2*pi*k(ct)*I(i))/ . .. 
betatau)*(max(g-1,m-l)-(m-1)))); 

end 

end 

if ct == M(i)/2 
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term815kleft(i) = term8sum(i); 


end 

end 

term815klr(i) = term8sum(i); 

term815kright(i) = term815klr(i)-term815kleft(i); 

end 

term815kleft 

term815kright 

term815klr 

exit 
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< M A T L A B (R) > 

(c) Copyright 1984-94 The MathWorks, Inc. 

All Rights Reserved 
Version 4.2 
Mar 29 1994 

Commands to get started: intro, demo, help help 

Commands for more information: help, whatsnew, info, subscribe 

malloc count matches malloc_debug_count 
malloc count matches malloc_debug_count 
malloc count matches malloc_debug_count 
malloc count matches malloc_debug_count 
malloc count matches malloc_debug_count 
malloc count matches malloc_debug_count 

term815kleft = 

-3.5253e+02 -2.1854e+02 -1.2475e+02 -4.5096e+01 7 . CASE: (b-=l,b+=0) 


term815kright = 

-3.5253e+02 -2.1854e+02 -1.2475e+02 -4.5096e+01 */. CASE: (b-=0,b+=l) 


term815klr = 


-7.0505e+02 -4.3708e+02 -2.4950e+02 -9.0192e+01 ’/. CASE: (b-=l,b+=l) 

X 

*/. 1 = 7 1 = 9 I = 12 I = 20 

•/. _ 

*/. VALUES OF NORMALIZED CHANNEL SPACING INTEGER I 


xxxxxxxxxxxxxxxxxx 

*/. TERM #8 '/. 

7 , betatau = 1500 '/, 

xxxxxxxxxxxxxxxxxx 


NOTE: THIS IS THE FILE term815k.out. ALL TEXT 
PRECEDED BY A X WAS ADDED LATER USING A 
TEXT EDITOR. 


36162000586 flops. 
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THESIS COMPUTER WORK 


TERM #8 

betatau = 2000 


COMMENTS: THE IDEA OF THIS PROGRAM IS TO COMPUTE TERM #8 

FOR THE VALUE OF betatau = 2000. WE WILL COMPUTE 
FOR THE FOUR VALUES OF I (AN INTEGER). THE VALUES OF I FOR 
betatau = 2000 are I = [8 9 12 20]. WE WILL COMPUTE THE VALUE OF 
TERM #8 THREE TIMES FOR EACH CORRESPONDING VALUE OF I. ONCE FOR 
ALL M/2 OF THE LOWER ACI CHANNELS (have negative summation index) 
BEING PACKED WITH l'S AND THE UPPER ACI CHANNELS ALL BEING PACKED 
WITH O'S. THIS IS THE CASE: (b-=l,b+=0). ONCE FOR ALL THE LOWER 
ACI CHANNELS BEING PACKED WITH ZERO'S AND THE UPPER ACI CHANNELS 
(have a positive summation index) BEING PACKED WITH l'S. THIS 
IS THE CASE: (b-=0,b+=l). FINALLY ONCE FOR BOTH THE UPPER AND 
LOWER ACI CHANNELS BEING PACKED WITH l'S. THIS IS THE CASE: 

(b-=l,b+=l). WE WILL USE BACKGROUND PROCESSING WHICH WILL 
DUMP ALL TWELVE VALUES TO A FILE CALLED term82k.mat. WE WILL 
THEN SINCE THE FILE IS IN BINARY FORMAT USE MATLAB INTERACTIVE 
COMMANDS AND THE TEXT EDITOR TO PRODUCE READABLE HARDCOPY 
RESULTS. WE WILL USE THESE VALUES LATER, IN ANOTHER PROGRAM 
TO COMPUTE THE DETECTION STATISTICS AND PROBABILITIES OF 
BIT ERROR FOR THE VARIOUS SIGNAL TO NOISE RATIOS. 


JOHN A. STUDER 
CPT, U.S. ARMY 
550-53-7181 


DATE LAST MODIFIED: 27 SEP 94 




format short e 
rho = 0.99; 
betatau = 2000; 

I = [8 9 12 20]; 

M = [248 220 164 98]; 
for i = 1:4 

k = [-M(i)/2:-l l:M(i)/2]; 
term8sum(i) = 0; 
for ct = l:M(i) 

for g = 1:betatau 
for m = l:betatau 

term8sum(i) = term8sum(i) + ((rho“((g-l)+(m-l)))/... 

(pi*k(ct)*I(i)))*(sin(2*pi*k(ct)*I(i)*... 
(l-((m-l)/betatau)))-sin(((2*pi*k(ct)*I(i))/... 
betatau)*(max(g-1,m-l)-(m-l)))); 

end 
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end 

if ct == M(i)/2 

term82kleft(i) = term8sum(i); 

end 

end 

term82klr(i) = term8sum(i); 

term82kright(i) = term82klr(i)-term82kleft(i); 

end 

save term82k term82kleft term82kright term82klr 
exit 
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< M A T L A B (R) > 

(c) Copyright 1984-93 The MathWorks, Inc. 

All Rights Reserved 
Version 4.1 
Jun 15 1993 

Commands to get started: intro, demo, help help 

Commands for more information: help, whatsnew, info, subscribe 

>> load term82k.mat 
>> who 

Your variables are: 

term82kleft term82klr term82kright 

>> term82kleft 
term82kleft = 

-351.8931 -283.0055 -163.9044 -60.1542 '/. CASE: (b-=l,b+=0) 

>> term82kright 
term82kright = 

-351.8931 -283.0055 -163.9044 -60.1542 ’/. CASE: (b-=0,b+=l) 

>> term82klr 
term82klr = 

-703.7863 -566.0110 -327.8088 -120.3084 '/. CASE: (b-=l,b+=l) 

•/. - 

y. 1 = 8 1 = 9 I = 12 I = 20 

*/. _ 

*/. 

7 , VALUES OF THE NORMALIZED CHANNEL SPACING INTEGER I 


y.y.y.y.y.y.y.y.y.y.y.y.y.y.y.y.y.y. 

y. TERM #8 y, 

'/, betatau = 2000 '/, 

y.y.y.y.y.y.y.y.y.y.y.y.y.y.y.y.y.y. 


NOTE: THIS IS THE FILE term82k.out CREATED FROM THE 

BINARY FILE term82k.mat USING MATLAB INTERACTIVE 
COMMANDS. ALL TEXT PRECEDED BY A */, WAS ADDED LATER 
USING A TEXT EDITOR. 
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7. THESIS COMPUTER WORK 

x 

x 

x 

y. _ 

•/. i i 

X I TERM #9 I 

X I betatau = 500 I 

X I_I 

X 

X COMMENTS: THE IDEA OF THIS PROGRAM IS TO COMPUTE TERM #9 
X FOR THE VALUE OF betatau = 500. WE WILL COMPUTE 

X FOR THE FOUR VALUES OF I (AN INTEGER). THE VALUES OF I FOR 

X betatau = 500 are I = [4 5 8 20]. WE WILL COMPUTE THE VALUE OF 

X TERM #9 THREE TIMES FOR EACH CORRESPONDING VALUE OF I. ONCE FOR 

X ALL M/2 OF THE LOWER ACI CHANNELS (have negative summation index) 

X BEING PACKED WITH l'S AND THE UPPER ACI CHANNELS ALL BEING PACKED 

X WITH O'S. THIS IS THE CASE: (b-=l,b+=0). ONCE FOR ALL THE LOWER 

X ACI CHANNELS BEING PACKED WITH ZERO’S AND THE UPPER ACI CHANNELS 

X (have a positive summation index) BEING PACKED WITH l’S. THIS 

X IS THE CASE: (b-=0,b+=l). FINALLY ONCE FOR BOTH THE UPPER AND 

X LOWER ACI CHANNELS BEING PACKED WITH l'S. THIS IS THE CASE: 

X (b-=l,b+=l). WE WILL USE BACKGROUND PROCESSING WHICH WILL 

X DUMP ALL TWELVE VALUES TO A FILE CALLED term9500.mat. 

X SINCE THE FILE IS IN A BINARY FORMAT WE WILL USE MATLAB 

X INTERACTIVE COMMANDS AND THE TEXT EDITOR TO PRODUCE READABLE 

X HARDCOPY RESULTS. AFTER THIS IS ALL DONE WE WILL 

X USE THESE VALUES LATER, IN ANOTHER PROGRAM TO COMPUTE THE 

X DETECTION STATISTICS AND PROBABILITIES OF BIT ERROR FOR THE 

X VARIOUS SIGNAL TO NOISE RATIOS. 

X 

X JOHN A. STUDER DATE LAST MODIFIED: 28 SEP 94 

X CPT, U.S. ARMY 
X 550-53-7181 


WNNNNNNNNNNNNNNNNNNNNNNNNNNNNNNNNNNNNNNNNNNNNNNNNNNNNNNNNNNNNNNNNNNNNNNNNNNN!^ 


format short e 
rho = 0.99; 
betatau = 500; 

I = [4 5 8 20] ; 

M = [124 98 60 24]; 
for i * 1:4 

k = [-M(i)/2:-1 l:M(i)/2]; 
term9sum(i) = 0; 
for ct = l:M(i) 

for g = l:betatau-l 
for m = 1:betatau 
if m-1 < g 

term9sum(i) = term9sum(i) + ((rho"(g+(m-l)))/(pi*k(ct)... 

*I(i)))*sin(2*pi*k(ct)*I(i)*((g-(m-l))... 
/betatau)); 
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end 


end 

end 

if ct == M(i)/2 

term95001eft(i) = term9sum(i); 

end 

end 

term95001r(i) = term9sum(i); 

term9500right(i) = term95001r(i)-term95001eft(i); 

end 

save term9500 term95001eft term9500right term95001r 
exit 
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<108 sp254201(SunOS) /kepler_u2/studer> matlab 

< M A T L A B (R) > 

(c) Copyright 1984-93 The MathWorks, Inc. 

All Rights Reserved 
Version 4.1 
Jun 15 1993 

Commands to get started: intro, demo, help help 

Commands for more information: help, whatsnew, info, subscribe 

>> load term9500.mat 
» who 

Your variables are: 

term95001eft term9500right 

term95001r 

>> term95001eft 

term95001eft = 

125.1608 80.4704 31.3156 4.7599 */. CASE: (b-=l,b+=0) 

» term9500right 
term9500right = 

125.1608 80.4704 31.3156 4.7599 ’/. CASE: (b-=0,b+=l) 

» term95001r 
term95001r = 

250.3216 160.9408 62.6313 9.5199 */. CASE: (b-=l,b+=l) 

» 

'/. 1=4 1 = 5 1 = 8 1 = 20 

/. _ 

X 

'/. VALUES OF NORMALIZED CHANNEL SPACING INTEGER I 





NOTE: THIS IS THE FILE term9500.out CREATED FROM THE BINARY 

FILE term9500.mat USING MATLAB INTERACTIVE COMMANDS. ALL 
TEXT PRECEDED BY A % WAS ADDED LATER USING A TEXT 
EDITOR. 
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THESIS COMPUTER WORK 


TERM #9 I 

betatau = 1000 | 


COMMENTS: THE IDEA OF THIS PROGRAM IS TO COMPUTE TERM #9 

FOR THE VALUE OF betatau = 1000. WE WILL COMPUTE 
FOR THE FOUR VALUES OF I (AN INTEGER). THE VALUES OF I FOR 
betatau = 1000 are I = [5 6 9 20] . WE WILL COMPUTE THE VALUE OF 
TERM #9 THREE TIMES FOR EACH CORRESPONDING VALUE OF I. ONCE FOR 
ALL M/2 OF THE LOWER ACI CHANNELS (have negative summation index) 
BEING PACKED WITH l'S AND THE UPPER ACI CHANNELS ALL BEING PACKED 
WITH O'S. THIS IS THE CASE: (b-=l,b+=0). ONCE FOR ALL THE LOWER 
ACI CHANNELS BEING PACKED WITH ZERO'S AND THE UPPER ACI CHANNELS 
(have a positive summation index) BEING PACKED WITH l'S. THIS 
IS THE CASE: (b-=0,b+=l). FINALLY ONCE FOR BOTH THE UPPER AND 
LOWER ACI CHANNELS BEING PACKED WITH l'S. THIS IS THE CASE: 

(b-=l,b+=l). WE WILL USE BACKGROUND PROCESSING WHICH WILL 
DUMP ALL TWELVE VALUES TO A FILE CALLED term91k.mat. 

SINCE THE FILE IS IN A BINARY FORMAT WE WILL USE MATLAB 
INTERACTIVE COMMANDS AND THE TEXT EDITOR TO PRODUCE READABLE 
HARDCOPY RESULTS. AFTER THIS IS ALL DONE WE WILL 
USE THESE VALUES LATER, IN ANOTHER PROGRAM TO COMPUTE THE 
DETECTION STATISTICS AND PROBABILITIES OF BIT ERROR FOR THE 
VARIOUS SIGNAL TO NOISE RATIOS. 


JOHN A. STUDER 
CPT, U.S. ARMY 
550-53-7181 


DATE LAST MODIFIED: 28 SEP 94 


format short e 
rho = 0.99; 
betatau = 1000; 

I = [5 6 9 20]; 

M = [198 164 110 48]; 
for i = 1:4 

k = [-M(i)/2:-1 l:M(i)/2]; 
term9sum(i) = 0; 
for ct = l:M(i) 

for g = l:betatau-l 
for m = 1:betatau 
if m-1 < g 

term9sum(i) = ■ 


term9sum(i) + ((rh.o~(g+(m-l)))/(pi*k(ct). 
*I(i)))*sin(2*pi*k(ct)*I(i)*((g-(m-l)) 
/betatau)); 
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end 

end 

end 

if ct == M(i)/2 

term91kleft(i) = term9sum(i); 

end 

end 

term91klr(i) = term9sum(i); 

term91kright(i) = term91klr(i)-term91kleft(i); 

end 

save term91k term91kleft term91kright term91klr 
exit 
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< M A T L A B (R) > 

(c) Copyright 1984-93 The MathWorks, Inc. 

All Rights Reserved 
Version 4.1 
Jun 15 1993 

Commands to get started: intro, demo, help help 

Commands for more information: help, whatsnew, info, subscribe 

» load term91k.mat 
» who 

Your variables are: 

term91kleft term91klr term91kright 

» term91kleft 


term91kleft = 


155.3047 109.6350 

49.5986 

9.9519 

7 . CASE: (b-=l,b+=0) 

» term91kright 

term91kright = 

155.3047 109.6350 

49.5986 

9.9519 

7 . CASE: (b-=0,b+=l) 

» term91klr 

term91klr = 

310.6094 219.2701 

99.1971 

19.9038 

7 . CASE: (b-=l,b+=l) 

» 

7 . 




7 . 1 = 5 1 = 6 

7 . 

1 = 9 

I = 20 



7 . VALUES OF THE NORMALIZED CHANNEL SPACING INTEGER I 


mnmmmmx 

7 . TERM #9 •/. 

'/, betatau = 1000 */, 


NOTE: THIS IS THE FILE term91k.out CREATED FROM THE 

BINARY FILE term91k.mat USING MATLAB INTERACTIVE 
COMMANDS. ALL TEXT PRECEDED BY A 7 . WAS ADDED LATER 
USING A TEXT EDITOR. 
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THESIS COMPUTER WORK 


TERM #9 

betatau = 1500 


COMMENTS: THE IDEA OF THIS PROGRAM IS TO COMPUTE TERM #9 

FOR THE VALUE OF betatau = 1500. WE WILL COMPUTE 
FOR THE FOUR VALUES OF I (AN INTEGER). THE VALUES OF I FOR 
betatau = 1500 are I = [7 9 12 20]. WE WILL COMPUTE THE VALUE OF 
TERM #9 THREE TIMES FOR EACH CORRESPONDING VALUE OF I. ONCE FOR 
ALL M/2 OF THE LOWER ACI CHANNELS (have negative summation index) 
BEING PACKED WITH l'S AND THE UPPER ACI CHANNELS ALL BEING PACKED 
WITH O’S. THIS IS THE CASE: (b-=l,b+=0). ONCE FOR ALL THE LOWER 
ACI CHANNELS BEING PACKED WITH ZERO’S AND THE UPPER ACI CHANNELS 
(have a positive summation index) BEING PACKED WITH l’S. THIS 
IS THE CASE: (b-=0,b+=l). FINALLY ONCE FOR BOTH THE UPPER AND 
LOWER ACI CHANNELS BEING PACKED WITH l’S. THIS IS THE CASE: 
(b-=l,b+=l). WE WILL USE BACKGROUND PROCESSING WHICH WILL 
DUMP ALL TWELVE VALUES TO A FILE CALLED term915k.mat. 

SINCE THE FILE IS IN A BINARY FORMAT WE WILL USE MATLAB 
INTERACTIVE COMMANDS AND THE TEXT EDITOR TO PRODUCE READABLE 
HARDCOPY RESULTS. AFTER THIS IS ALL DONE WE WILL 
USE THESE VALUES LATER, IN ANOTHER PROGRAM TO COMPUTE THE 
DETECTION STATISTICS AND PROBABILITIES OF BIT ERROR FOR THE 
VARIOUS SIGNAL TO NOISE RATIOS. 


JOHN A. STUDER 
CPT, U.S. ARMY 
550-53-7181 


DATE LAST MODIFIED: 28 SEP 94 




format short e 
rho = 0.99; 
betatau = 1500; 

I = [7 9 12 20]; 

M = [212 164 124 74] ; 
for i = 1:4 

k = [-M(i)/2:-l l:M(i)/2]; 
term9sum(i) = 0; 
for ct = l:M(i) 

for g = 1:betatau-1 
for m = 1:betatau 
if m-1 < g 

term9sum(i) = term9sum(i) + ((rho~(g+(m-l)))/(pi*k(ct)... 

*I(i)))*sin(2*pi*k(ct)*I(i)*((g-(m-l))... 
/betatau)); 
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end 


end 

end 

if ct == M(i)/2 

term915kleft(i) = term9sum(i); 

end 

end 

term915klr(i) = term9snm(i); 

term915kright(i) = term915klr(i)-term915kleft(i) 

end 

save term915k term915kleft term915kright term915klr 
exit 
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< M A T L A B (R) > 

(c) Copyright 1984-93 The MathWorks, Inc. 

All Rights Reserved 
Version 4.1 
Jun 15 1993 

Commands to get started: intro, demo, help help 

Commands for more information: help, whatsnew, info, subscribe 

>> load term915k 
>> who 

Your variables are: 

term915kleft term915kright 

term915klr 

>> term915kleft 


term915kleft = 

117.9019 73.0900 41.7221 15.0822 '/. CASE: (b-=l,b+=0) 

>> term915kright 
term915kright = 

117.9019 73.0900 41.7221 15.0822 '/. CASE: (b-=0,b+=l) 

>> term915klr 
term915klr = 

235.8037 146.1800 83.4441 30.1644 ’/. CASE: (b-=l,b+=l) 

•/. 

•/. 1 = 7 1 = 9 I = 12 I = 20 

’/. _ 

7 . 

7. VALUES OF NORMALIZED CHANNEL SPACING INTEGER I 


mmmmmnn 

*/, TERM #9 ’/. 

'/, betatau = 1500 '/, 

mmmmxmm 


NOTE: THIS IS THE FILE term915k.out CREATED FROM THE 

BINARY FILE term915k.mat USING MATLAB INTERACTIVE 
COMMANDS. ALL TEXT PRECEDED BY A '/, WAS ADDED LATER 
USING A TEXT EDITOR. 
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THESIS COMPUTER WORK 


TERM #9 | 

betatau = 2000 j 


COMMENTS: THE IDEA OF THIS PROGRAM IS TO COMPUTE TERM #9 

FOR THE VALUE OF betatau = 2000. WE WILL COMPUTE 
FOR THE FOUR VALUES OF I (AN INTEGER). THE VALUES OF I FOR 
betatau = 2000 are I = [8 9 12 20]. WE WILL COMPUTE THE VALUE OF 
TERM #9 THREE TIMES FOR EACH CORRESPONDING VALUE OF I. ONCE FOR 
ALL M/2 OF THE LOWER ACI CHANNELS (have negative summation index) 
BEING PACKED WITH l'S AND THE UPPER ACI CHANNELS ALL BEING PACKED 
WITH O'S. THIS IS THE CASE: (b-=l,b+=0). ONCE FOR ALL THE LOWER 
ACI CHANNELS BEING PACKED WITH ZERO'S AND THE UPPER ACI CHANNELS 
(have a positive summation index) BEING PACKED WITH l'S. THIS 
IS THE CASE: (b-=0,b+=l). FINALLY ONCE FOR BOTH THE UPPER AND 
LOWER ACI CHANNELS BEING PACKED WITH l'S. THIS IS THE CASE: 
(b-=l,b+=l). WE WILL USE BACKGROUND PROCESSING WHICH WILL 
DUMP ALL TWELVE VALUES TO A FILE CALLED term92k.mat. 

SINCE THE FILE IS IN A BINARY FORMAT WE WILL USE MATLAB 
INTERACTIVE COMMANDS AND THE TEXT EDITOR TO PRODUCE READABLE 
HARDCOPY RESULTS. AFTER THIS IS ALL DONE WE WILL 
USE THESE VALUES LATER, IN ANOTHER PROGRAM TO COMPUTE THE 
DETECTION STATISTICS AND PROBABILITIES OF BIT ERROR FOR THE 
VARIOUS SIGNAL TO NOISE RATIOS. 


JOHN A. STUDER 
CPT, U.S. ARMY 
550-53-7181 


DATE LAST MODIFIED: 28 SEP 94 


format short e 
rho = 0.99; 
betatau = 1500; 

I = [8 9 12 20]; 

M = [248 220 164 98]; 
for i = 1:4 

k = [-M(i)/2:-l l:M(i)/2]; 
term9sum(i) = 0; 
for ct = l:M(i) 

for g = l:betatau-l 
for m = 1:betatau 
if m-1 < g 

term9sum(i) = term9sum(i) + ((rho"(g+(m-l)))/(pi*k(ct)... 

*I(i)))*sin(2*pi*k(ct)*I(i)*((g-(m-l))... 
/betatau)); 
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end 

end 

end 

if ct == M(i)/2 

term92kleft(i) = term9sum(i); 

end 

end 

term92klr(i) = term9sum(i); 

term92kright(i) = term92klr(i)-term92kleft(i); 

end 

save term92k term92kleft term92kright term92klr 
exit 
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< M A T L A B (R) > 

(c) Copyright 1984-93 The MathWorks, Inc. 

All Rights Reserved 
Version 4.1 
Jun 15 1993 

Commands to get started: intro, demo, help help 

Commands for more information: help, whatsnew, info, subscribe 

>> load term92k.mat 
» who 

Your variables are: 

term92kleft term92klr term92kright 

>> term92kleft 


term92kleft = 

91.5258 73.0261 41.6758 15.0544 */. CASE: (b-=l,b+=0) 

» term92kright 
term92kright = 

91.5258 73.0261 41.6758 15.0544 */. CASE: (b-=0,b+=l) 

» term92klr 


term92klr = 

183.0515 146.0523 83.3516 30.1088 */, CASE: (b-=l,b+=l) 


» 

% 

’/. 

7 . 

7 . 

7 . 


1=8 1 = 9 1 = 12 1 = 20 


VALUES OF THE NORMALIZED CHANNEL SPACING INTEGER I 


*/. TERM #9 '/. 

*/, betatau = 2000 '/. 


NOTE: THIS IS THE FILE term92k.out CREATED FROM THE 

BINARY FILE term92k.mat USING MATLAB INTERACTIVE 
COMMANDS. ALL TEXT PRECEDED BY A ’/. WAS ADDED LATER 
USING A TEXT EDITOR. 
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7 . 

7 . 

7 . 

7 . 

*/. 

7 . 

7 . 

7 . 

7 . 

7 . 

7 . 

7 . 

7 . 

7 . 

7 . 

7 . 

7 . 

7 . 

7 . 

7 . 

7 . 

7 . 

7 . 

7 . 

7 . 

7 . 

7 . 

7 . 

7 . 

7 . 


THESIS COMPUTER WORK 


| TERM #10 | 

| betatau = 500 | 


COMMENTS: THE IDEA OF THIS PROGRAM IS TO COMPUTE TERM #10 
FOR THE VALUE OF betatau = 500. WE WILL COMPUTE 
FOR THE FOUR VALUES OF I (AN INTEGER). THE VALUES OF I FOR 
betatau = 500 are I = [4 5 8 20]. WE WILL COMPUTE THE VALUE OF 
TERM #10 THREE TIMES FOR EACH CORRESPONDING VALUE OF I. ONCE FOR 
ALL M/2 OF THE LOWER ACI CHANNELS (have negative summation index) 
BEING PACKED WITH l'S AND THE UPPER ACI CHANNELS ALL BEING PACKED 
WITH O'S. THIS IS THE CASE: (b-=l,b+=0). ONCE FOR ALL THE LOWER 
ACI CHANNELS BEING PACKED WITH ZERO'S AND THE UPPER ACI CHANNELS 
(have a positive summation index) BEING PACKED WITH l'S. THIS 
IS THE CASE: (b-=0,b+=l). FINALLY ONCE FOR BOTH THE UPPER AND 
LOWER ACI CHANNELS BEING PACKED WITH l’S. THIS IS THE CASE: 

(b-=l,b+=l). WE WILL USE BACKGROUND PROCESSING WHICH WILL 
DUMP ALL TWELVE VALUES TO A FILE CALLED terml0500.mat. 

SINCE THE FILE IS IN A BINARY FORMAT, WE WILL USE MATLAB 
INTERACTIVE COMMANDS AND THE TEXT EDITOR TO PRODUCE READABLE 
HARDCOPY RESULTS. AFTER THIS IS ALL DONE, WE WILL 
USE THESE VALUES LATER , IN ANOTHER PROGRAM, TO COMPUTE THE 
DETECTION STATISTICS AND PROBABILITIES OF BIT ERROR FOR THE 
VARIOUS SIGNAL TO NOISE RATIOS. 


7 . 

7. JOHN A. STUDER DATE LAST MODIFIED: 29 SEP 94 

7 . CPT, U.S. ARMY 
7 . 550-53-7181 

•/ •/ v •/ >/ •/•/•/ */ v v •/ •/ •/ •/ •/ •/ •/ •/ •/ •/ •/ •/ •/ v v"/ v v v v v v v v v v v v v v v v v v v v •/ v v */ v v v v v v y v. v. v. v. •/ v. v. v. v. */. 7 y. •/ v. */. y. v. 7 */.'/. 


format short e 
rho = 0.99; 
betatau = 500; 

I = [4 5 8 20] ; 

M = [124 98 60 24]; 
for i = 1:4 

k = [-M(i)/2:-l l:M(i)/2]; 
termlOsum(i) = 0; 
for ct = l:M(i) 

for g = betatau: (2*betatau)-1 
for m = 1:betatau 

termlOsum(i) = termlOsum(i) + ((rho~(g+(m-l)))/... 

(pi*k(ct)*I(i)))*(sin(2*pi*k(ct)*I(i)*... 
(l-((m-l)/betatau)))-sin(((2*pi*k(ct)*I(i))/... 
betatau)*(max(g-betatau,m-l)-(m-l)))); 
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end 

end 

if ct == M(i)/2 

terml05001eft(i) = termlOsum(i); 

end 

end 

terml05001r(i) = termlOsum(i); 

terml0500right(i) = terml05001r(i)-terml05001eft(i); 

end 

save terml0500 terml05001eft terml0500right terml05001r 
exit 
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< M A T L A B (R) > 

(c) Copyright 1984-93 The MathWorks, Inc. 

All Rights Reserved 
Version 4.1 
Jun 15 1993 

Commands to get started: intro, demo, help help 

Commands for more information: help, whatsnew, info, subscribe 

>> load terml0500.mat 
>> who 

Your variables are: 


terml05001eft terml0500right 

terml05001r 


>> terml05001eft 
terml05001eft = 

-2.4375 -1.5672 -0.6099 -0.0927 '/. CASE: (b-=l,b+=0) 

>> terml0500right 
terml0500right = 

-2.4375 -1.5672 -0.6099 -0.0927 */. CASE: (b-=0,b+=l) 

>> terml05001r 
terml05001r = 

-4.8750 -3.1343 -1.2197 -0.1854 */. CASE: (b-=l,b+=l) 


» 

y. 

7 , 1 = 4 1 = 5 1 = 8 1 = 20 

•/. _ 

y. 

y. VALUES OF NORMALIZED CHANNEL SPACING INTEGER I 


•/.y.y.y.y.y.y.y.y.y.y.y.y.y.y.y.y.y.y. 

y. TERM #10 y. 


y, betatau = 500 % 

y.y.y.y.y.y.y.y.y.y.y.y.y.y.y.y.y.y.y. 


NOTE: THIS IS THE FILE terml0500.out CREATED FROM THE BINARY FILE 
terml0500.mat USING MATLAB INTERACTIVE COMMANDS. ALL TEXT 
PRECEDED BY A */. WAS ADDED LATER USING A TEXT EDITOR. 
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•/. 

•/. 

7 . 

7 , 

7 . 

7 . 

7 . 

7 . 

7 . 

7 . 

*/. 

•/. 

*/. 

*/. 

•/. 

•/. 

•/. 

*/. 

•/. 

7 . 

*/. 

7 . 

*/. 

•/. 

*/, 

7 . 

*/. 

•/. 

*/. 

*/. 


THESIS COMPUTER WORK 


TERM #10 | 

betatau = 1000 I 


COMMENTS: THE IDEA OF THIS PROGRAM IS TO COMPUTE TERM #10 
FOR THE VALUE OF betatau = 1000. WE WILL COMPUTE 
FOR THE FOUR VALUES OF I (AN INTEGER). THE VALUES OF I FOR 
betatau = 1000 are I = [5 6 9 20]. WE WILL COMPUTE THE VALUE OF 
TERM #10 THREE TIMES FOR EACH CORRESPONDING VALUE OF I. ONCE FOR 
ALL M/2 OF THE LOWER ACI CHANNELS (have negative summation index) 
BEING PACKED WITH l'S AND THE UPPER ACI CHANNELS ALL BEING PACKED 
WITH O'S. THIS IS THE CASE: (b-=l,b+=0). ONCE FOR ALL THE LOWER 
ACI CHANNELS BEING PACKED WITH ZERO'S AND THE UPPER ACI CHANNELS 
(have a positive summation index) BEING PACKED WITH l'S. THIS 
IS THE CASE: (b-=0,b+=l). FINALLY ONCE FOR BOTH THE UPPER AND 
LOWER ACI CHANNELS BEING PACKED WITH l'S. THIS IS THE CASE: 
(b-=l,b+=l). WE WILL USE BACKGROUND PROCESSING WHICH WILL 
DUMP ALL TWELVE VALUES TO A FILE CALLED termlOlk.mat. 

SINCE THE FILE IS IN A BINARY FORMAT, WE WILL USE MATLAB 
INTERACTIVE COMMANDS AND THE TEXT EDITOR TO PRODUCE READABLE 
HARDCOPY RESULTS. AFTER THIS IS ALL DONE, WE WILL 
USE THESE VALUES LATER , IN ANOTHER PROGRAM, TO COMPUTE THE 
DETECTION STATISTICS AND PROBABILITIES OF BIT ERROR FOR THE 
VARIOUS SIGNAL TO NOISE RATIOS. 

JOHN A. STUDER DATE LAST MODIFIED: 29 SEP 94 

CPT, U.S. ARMY 

550-53-7181 


format short e 
rho = 0.99; 
betatau = 1000; 

I = [5 6 9 20] ; 

M = [198 164 110 48] ; 
for i = 1:4 

k = [—M(i)/2:—1 l:M(i)/2]; 
termlOsum(i) = 0; 
for ct = l:M(i) 

for g = betatau:(2*betatau)-l 
for m = 1:betatau 

termlOsum(i) = termlOsum(i) + ((rho~(g+(m-l)))/... 

(pi*k(ct)*I(i)))*(sin(2*pi*k(ct)*I(i)*... 
(l-((m-l)/betatau)))-sin(((2*pi*k(ct)*I(i))/... 
betatau)*(max(g-betatau,m-l)-(m-l)))); 
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end 


end 

if ct == M(i)/2 

termlOlkleft(i) = termlOsum(i); 

end 

end 

termlOlklr(i) = termlOsum(i); 

termlOlkright(i) = termlOlklr(i)-termlOlkleft(i); 

end 

save termlOlk termlOlkleft termlOlkright termlOlklr 
exit 
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(c) Copyright 1984-93 The MathWorks, Inc. 

All Rights Reserved 
Version 4.1 
Jun 15 1993 

Commands to get started: intro, demo, help help 

Commands for more information: help, whatsnew, info, subscribe 

>> load termlOlk.mat 
>> who 

Your variables are: 

termlOlkleft termlOlkright 

termlOlklr 

>> termlOlkleft 

termlOlkleft = 

-0.0200 -0.0142 -0.0064 -0.0013 % CASE: (b-=l,b+=0) 

>> termlOlkright 
termlOlkright = 

-0.0200 -0.0142 -0.0064 -0.0013 7. CASE: (b-=0,b+=l) 

>> termlOlklr 
termlOlklr = 



NOTE: THIS IS THE FILE termlOlk.out CREATED FROM THE 

BINARY FILE termlOl.mat USING MATLAB INTERACTIVE 
COMMANDS. ALL TEXT PRECEDED BY A 7 . WAS ADDED LATER 
USING A TEXT EDITOR. 


172 





THESIS COMPUTER WORK 


TERM #10 
betatau = 1500 


COMMENTS: THE IDEA OF THIS PROGRAM IS TO COMPUTE TERM #10 
FOR THE VALUE OF betatau = 1500. WE WILL COMPUTE 
FOR THE FOUR VALUES OF I (AN INTEGER). THE VALUES OF I FOR 
betatau = 1500 are I = [7 9 12 20]. WE WILL COMPUTE THE VALUE OF 
TERM #10 THREE TIMES FOR EACH CORRESPONDING VALUE OF I. ONCE FOR 
ALL M/2 OF THE LOWER ACI CHANNELS (have negative summation index) 
BEING PACKED WITH l'S AND THE UPPER ACI CHANNELS ALL BEING PACKED 
WITH O'S. THIS IS THE CASE: (b-=l,b+=0). ONCE FOR ALL THE LOWER 
ACI CHANNELS BEING PACKED WITH ZERO'S AND THE UPPER ACI CHANNELS 
(have a positive summation index) BEING PACKED WITH l’S. THIS 
IS THE CASE: (b-=0,b+=l). FINALLY ONCE FOR BOTH THE UPPER AND 
LOWER ACI CHANNELS BEING PACKED WITH l'S. THIS IS THE CASE: 

(b-=l,b+=l). WE WILL USE BACKGROUND PROCESSING WHICH WILL 
DUMP ALL TWELVE VALUES TO A FILE CALLED terml015k.mat. 

SINCE THE FILE IS IN A BINARY FORMAT, WE WILL USE MATLAB 
INTERACTIVE COMMANDS AND THE TEXT EDITOR TO PRODUCE READABLE 
HARDCOPY RESULTS. AFTER THIS IS ALL DONE, WE WILL 
USE THESE VALUES LATER , IN ANOTHER PROGRAM, TO COMPUTE THE 
DETECTION STATISTICS AND PROBABILITIES OF BIT ERROR FOR THE 
VARIOUS SIGNAL TO NOISE RATIOS. 


JOHN A. STUDER 
CPT, U.S. ARMY 
550-53-7181 


DATE LAST MODIFIED: 29 SEP 94 




format short e 
rho = 0.99; 
betatau = 1500; 

I = [7 9 12 20]; 

M = [212 164 124 74]; 
for i = 1:4 

k = [—M(i)/2:—1 l:M(i)/2]; 
termlOsum(i) = 0; 
for ct = 1:M(i) 

for g = betatau:(2*betatau)-1 
for m = l:betatau 

termlOsum(i) = termlOsum(i) + ((rho“(g+(m-l)))/... 

(pi*k(ct)*I(i)))*(sin(2*pi*k(ct)*I(i)*... 
(l-((m-l)/betatau)))-sin(((2*pi*k(ct)*I(i))/... 
betatau)*(max(g-betatau,m-l)-(m-l)))); 

end 

end 


173 









if ct == M(i)/2 

term!015kleft(i) = termlOsum(i); 

end 

end 

terml015klr(i) = termlOsum(i); 

terml015kright(i) = terml015klr(i)-terml015kleft(i); 

end 

save terml015k terml015kleft terml015kright terml015klr 
exit 


174 



< M A T L A B (R) > 

(c) Copyright 1984-93 The MathWorks, Inc. 

All Rights Reserved 
Version 4.1 
Jun 15 1993 

Commands to get started: intro, demo, help help 

Commands for more information: help, whatsnew, info, subscribe 

>> load terml015k.mat 
>> who 

Your variables are: 

terml015kleft terml015kright 

terml015klr 

>> terml015kleft 

terml015kleft = 

1.0e-04 * 

-1.0000 -0.6199 -0.3539 -0.1279 ’/, CASE: (b-=l,b+=0) 

>> terml015kright 
terml015kright = 

1.0e-04 * 

-1.0000 -0.6199 -0.3539 -0.1279 ’/. CASE: (b-=0,b+=l) 

>> terml015klr 
terml015klr = 

1.0e-03 * 

-0.2000 -0.1240 -0.0708 -0.0256 ’/. CASE: (b-=l,b+=l) 

» 

*/. 

’/, 1 = 7 1 = 9 I = 12 I = 20 

•/. _ 

7 . 

'/. VALUES OF NORMALIZED CHANNEL SPACING INTEGER I 
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NOTE: 


mnmmmmr/. 

% TERM #10 */. 

7 . betatau = 1500 ’/, 


THIS IS THE FILE terml015k.out CREATED FROM THE 
BINARY FILE terml015k.mat USING MATLAB INTERACTIVE 
COMMANDS. ALL TEXT PRECEDED BY A */. WAS ADDED LATER 
USING A TEXT EDITOR. 
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THESIS COMPUTER WORK 


7 . 

7 . 

7 . 

7 . 

7 . 

7 . 

7 . 

7 . 

7 . 

7 . 

7 . 

7 . 

7 . 

7 . 

7 . 

7 . 

7 . 

7 . 

7 . 

7 . 

7 . 

7 . 

7 . 

7 . 

7 . 

7 . 

7 . 

7 . 

7 . 

7 . 

7 . 

7 . 

7 . 


TERM #10 
betatau = 2000 


I. 


COMMENTS: THE IDEA OF THIS PROGRAM IS TO COMPUTE TERM #10 
FOR THE VALUE OF betatau = 2000. WE WILL COMPUTE 
FOR THE FOUR VALUES OF I (AN INTEGER). THE VALUES OF I FOR 
betatau = 2000 are I = [8 9 12 20]. WE WILL COMPUTE THE VALUE OF 
TERM #10 THREE TIMES FOR EACH CORRESPONDING VALUE OF I. ONCE FOR 
ALL M/2 OF THE LOWER ACI CHANNELS (have negative summation index) 
BEING PACKED WITH l’S AND THE UPPER ACI CHANNELS ALL BEING PACKED 
WITH O’S. THIS IS THE CASE: (b-=l,b+=0). ONCE FOR ALL THE LOWER 
ACI CHANNELS BEING PACKED WITH ZERO’S AND THE UPPER ACI CHANNELS 
(have a positive summation index) BEING PACKED WITH l’S. THIS 
IS THE CASE: (b-=0,b+=l). FINALLY ONCE FOR BOTH THE UPPER AND 
LOWER ACI CHANNELS BEING PACKED WITH l’S. THIS IS THE CASE: 

(b-=l,b+=l). WE WILL USE BACKGROUND PROCESSING WHICH WILL 
DUMP ALL TWELVE VALUES TO A FILE CALLED terml02k.mat. 

SINCE THE FILE IS IN A BINARY FORMAT, WE WILL USE MATLAB 
INTERACTIVE COMMANDS AND THE TEXT EDITOR TO PRODUCE READABLE 
HARDCOPY RESULTS. AFTER THIS IS ALL DONE, WE WILL 
USE THESE VALUES LATER , IN ANOTHER PROGRAM, TO COMPUTE THE 
DETECTION STATISTICS AND PROBABILITIES OF BIT ERROR FOR THE 
VARIOUS SIGNAL TO NOISE RATIOS. 


JOHN A. STUDER 
CPT, U.S. ARMY 


DATE LAST MODIFIED: 28 NOV 94 


7 , 550-53-7181 

7 . 7 . 7 . 7 . 7 . 7 . 7 . 7 . 7 . 7 . 7 . 7 . 7 //. 7 . 7 . 7 . 7 . 7 . 7 . 7 . 7 . 7 . 7 . 7 . 7 . 7 . 7 . 7 . 7 . 7 . 7 . 7 . 7 //. 7 . 7 . 7 . 7 . 7 . 7 . 7 //. 7 . 7 . 7 . 7 . 7 . 7 . 7 . 7 . 7 . 7 . 7 . 7 . 7 . 7 . 7 . 7 . 7 . 7 . 7 . 7 . 7 . 7 . 7 . 7 . 7 . 7 . 7 . 7 . 7 . 7 . 7 . 7 . 7 . 7 . 7 . 
format short e 


rho = 0.99; 
betatau = 2000; 

I = [8 9 12 20]; 

M = [248 220 164 98]; 


for i = 1:4 

k = [—M(i)/2:—1 l:M(i)/2]; 
termlOsum(i) = 0; 
for ct = l:M(i) 

for g = betatau:(2*betatau)-1 
for m = 1:betatau 

termlOsum(i) = termlOsum(i) + ((rho“(g+(m-l)))/... 

(pi*k(ct)*I(i)))*(sin(2*pi*k(ct)*I(i)*... 
(l-((m-l)/betatau)))-sin(((2*pi*k(ct)*I(i))/... 
betatau)*(max(g-betatau,m-l)-(m-l)))); 










end 


end 

if ct == M(i)/2 

terml02kleft(i) = termlOsum(i); 

end 

end 

terml02klr(i) = termlOsum(i); 

terml02kright(i) = terml02klr(i)-terml02kleft(i); 

end 

save terml02k terml02kleft terml02kright terml02klr 
exit 


178 



>> terml02kleft8 


terml02kleft8 = 

-6.5584e-07 
>> terml02kright8 
terml02kright8 = 
-6.5584e-07 
» terml02klr8 
terml02klr8 = 
-1.3117e-06 

» 

>> terml02kleft9 
terml02kleft9 = 

-5.2745e-07 
>> terml02kright9 
terml02kright9 = 

-5.2745e-07 
>> terml02klr9 
terml02klr9 = 

-1.0549e-06 

» 

>> terml02kleftl2 
terml02kleftl2 = 

-3.0548e-07 
» terml02krightl2 
terml02krightl2 = 
-3.0548e-07 


'/. CASE: (b—= 1 ,b+=0) 


•/. CASE: (b-=0,b+=l) 


‘/. CASE: (b-=l ,b+=l) 


*/. CASE: (b-=l ,b+=0) 


•/, CASE: (b-=0,b+=l) 


7 . CASE: (b-=l,b+=l) 


*/. CASE: (b—=1 ,b+=0) 


*/. CASE: (b-=0,b+=l) 


mmm 
7 . 1 = 8 •/, 
mmm 


mmm 

•/. I = 9 •/. 

mmm 


mmmy. 

*/. I = 12 ’/. 

7:i:i:i:i:i:i:i:il 




>> terml02klrl2 


terml02klrl2 = 

-6.1096e-07 '/. CASE: (b-=l,b+=l) 

» 

>> terml02kleft20 
terml02kleft20 = 

-1.1211e-07 
» terml02kright20 
terml02kright20 = 

-1.1211e-07 
>> terml02klr20 
terml02klr20 = 

-2.2423e-07 '/. CASE: (b-=l,b+=l) 

» 


’/. CASE: (b—=1 ,b+=0) 


7. CASE: (b-=0,b+=l) 


mxmm 

7 . i = 20 y. 
y.y.y.y.y.y.*/.y.y.y. 


y.y.y.y.*/.y.y.y.y.y.y.y.y.y.y.y. •/.•/. 

y. term #10 y. 

'/, betatau = 2000 '/. 

y.y.y.y.y.y.y.y.y.y.y.y.y.y.y.y.y.y. 


NOTE: THERE WERE SEVERAL MISHAPS IN COMPUTING THIS TERM. POWER OUTAGES, 

MY ACCOUNT BEING SHUT OFF,ETC. IN THE INTERESTS OF BEING ABLE TO 
COMPUTE RESULTS IN A TIMELY MANNER, FOUR COPIES OF THE PROGRAM FOR 
TERM #10, betatau = 2000 WERE MADE. THE I AND M BLOCK WERE MODIFIED 
TO HOLD ONLY ONE VALUE OF I/M, AND THE OUTER LOOP WAS REDUCED TO 
ONE ITERATION. THEN, THE FOUR PROGRAMS WERE RUN ON FOUR SEPARATE 
WORKSTATIONS. IN THE INTERESTS OF BREVITY WE ONLY SHOW THE SINGLE 
PROGRAM FOR COMPUTING THE TERM AS THE FOUR COPIES ARE ESSENTIALLY 
THE EXACT SAME PROGRAM BUT THEY DO ONLY ONE OUTER LOOP EACH. 

THE RESULTS OF THE FOUR PROGRAMS ARE TABULATED HERE USING THE 
A TEXT EDITOR. 
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THESIS COMPUTER WORK 

GENERATION OF GRAPHS 
FOR 

betatau = 500 

DATE: 19 OCTOBER 1994 

COMMENTS: THIS PROGRAM CONTAINS THE RESULTS OF ALL TEN TERMS OF 
APPENDIX B RESULTING IN THE PROGRAM PRODUCING A 
COMPLETE, ACCURATE GRAPH OF THE SYSTEM'S PROBABILITY 
OF BIT ERROR FOR ALL FOUR VALUES OF I. 

I = [4 5 8 20] 


INNNNNNNNNNNNNNNNNNNNNNNNNNNNNNNNNNNNNNNNNNNNNNNNNNNNNNNNNNNNNNNNNNNNNNNNNNNNNNt 


format short e 
rho = 0.99; 
c = (1-rho)“2; 

I = [4 5 8 20] ; 

M = [124 98 60 24]; 

bitmatrix =[000; 001; 010; Oil; 100; 101; 110; 111]; 

term5 = [l.5423e+03 0.9419e+03 0.3307e+03 0.0431e+03 

1.5423e+03 0.9419e+03 0.3307e+03 0.0431e+03 

1.3825e+03 0.9001e+03 0.3576e+03 0.0580e+03]; 

term8 = [-3.7098e+02 -2.3852e+02 -9.2820e+01 -1.4109e+01 

-3.7098e+02 -2.3852e+02 -9.2820e+01 -1.4109e+01 
-7.4196e+02 -4.7703e+02 -1.8564e+02 -2.8217e+0l] ; 

term9 = [125.1608 80.4704 31.3156 4.7599 
125.1608 80.4704 31.3156 4.7599 
250.3216 160.9408 62.6313 9.5199]; 

termlO = [-2.4375 -1.5672 -0.6099 -0.0927 
-2.4375 -1.5672 -0.6099 -0.0927 
-4.8750 -3.1343 -1.2197 -0.1854]; 

for ct = 1:4 
for i = 1:8 

for boo = [1 0] 

XI = boo*c*7.0510e+03; 

X2 = bitmatrix(i,3)*c*959.5662; 

X3 = bitmatrix(i,3)*c*l.2211e+01; 

X4 = bitmatrix(i,3)*c*3.0440e-01; 
if bitmatrix(i,1)==1 k bitmatrix(i,2)==0 
X5 = c*term5(l,ct); 

elseif bitmatrix(i,l)==0 k bitmatrix(i,2)==1 
X5 = c*term5(2,ct); 

elseif bitmatrix(i,1)==1 k bitmatrix(i,2)==1 
X5 = c*term5(3,ct); 


y. 

•/. 

•/. 

y. 

y. 

y. 

y. 

y. 

y. 

y. 

•/. 

y. 

•/. 

y. 







else 

X5 = 0; 

end 

X6 = boo*bitmatrix(i,3)*c*l.8585e+03; 

X7 = boo*bitmatrix(i,3)*c*9.2657e+01; 
if boo == 0 
X8 = 0; 

elseif bitmatrix(i,1)==1 & bitmatrix(i,2)==0 
X8 = c*term8(l,ct); 

elseif bitmatrix(i,1)==0 & bitmatrix(i,2)==1 
X8 = c*term8(2,ct); 

elseif bitmatrix(i, 1) — 1 & bitmatrix(i,2)==1 
X8 = c*term8(3,ct); 
else 

X8 = 0; 

end 

if bitmatrix(i,3) == 0 
X9 = 0; 

elseif bitmatrix(i,1)==1 & bitmatrix(i,2)==0 
X9 = c*term9(l,ct); 

elseif bitmatrix(i,1)==0 & bitmatrix(i,2)==1 
X9 = c*term9(2,ct); 

elseif bitmatrix(i,1)==1 & bitmatrix(i,2)==1 
X9 = c*term9(3,ct); 
else 

X9 = 0; 

end 

if bitmatrix(i,3) ==0 
X10 = 0; 

elseif bitmatrix(i,1)==1 & bitmatrix(i,2)==0 
X10 = c*termlO(l,ct); 

elseif bitmatrix(i,1)==0 & bitmatrix(i,2)==1 
X10 = c*terml0(2,ct); 

elseif bitmatrix(i,1)==1 & bitmatrix(i,2)==1 
X10 = c*terml0(3,ct); 
else 

X10 = 0; 

end 

INTERSUM = X1+X2+X3+X4+X5+X6+X7+X8+X9+X10; 
if boo==l 

XlBETA(ct.i) = INTERSUM; 
else 

XOBETA(ct,i) = INTERSUM; 

end 

end 

end 

XlMIN(ct) = min(XlBETA(ct,:)); 

XOMAX(ct) = max(XOBETA(ct,:)); 

VT(ct) = (XOMAX(ct) + X1MIN(ct))/2; 
end 

RPSQR_TDIVNO_DB = 10:.01:20; 
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RPSQR_TDIVNO = 10.~(RPSQR_TDIVN0_DB*0.10); 

SINGCHAN = 0.5*erfc(RPSQR_TDIVN0/8~0.5); 
for ct = 1:4 

PE(ct,1:1001) = zeros(l,1001); 
for i = 1:8 

PE(ct,:)=PE(ct, : )+0. 25*erfc((RPSQR_TDIVN0/2~0.5)*(XlBETA(ct,i)-VT(ct)))... 
+0.25*erfc((RPSQR_TDIVN0/2~0.5)*(VT(ct)-X0BETA(ct,i))); 

end 

end 

PEFINAL = PE/8; 
figure(1) 

semilogy(RPSQR_TDIVNO_DB,SINGCHAN,'— ’ ,RPSQR.TDIVNO.DB,PEFINAL(1,:),... 

RPSQR_TDIVNO_DB,PEFINAL(2,:),RPSQR_TDIVNO_DB,PEFINAL(3,:), . .. 
RPSQR_TDIVNO_DB,PEFINAL(4,:)) 
xlabelC'Z (dB)'); 
ylabel('Pb’); 
axis([l0 17 10~(—15) 1]) 
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THESIS COMPUTER WORK 


7 , 

I 

7 . 


7 . 
7 . 
7 ,; 

7 . 

7 . 

7 . 

7 . 

7 . 

7 , 

7 . 

7 . 

7 , 

7 . 

7 . 

7 . 

7 . 


GENERATION OF GRAPHS 
FOR 

betatau = 1000 


DATE: 21 OCTOBER 1994 

COMMENTS: THIS GRAPHING PROGRAM DOES NOT CONTAIN THE RESULTS FOR 
TERM #5, APPENDIX B DUE TO AN INABILITY TO ARRIVE AT A 
SOLUTION. THERE WERE MEMORY PROBLEMS WHEN A MATRIX SOLUTION 
WAS TRIED TO SPEED CALCULATION OF THIS TERM, AND THE 
TIME REQUIRED TO CALCULATE USING SIMPLE LOOPS WAS 
PROHIBITIVE. THUS, ON THE RESULTING GRAPH PRODUCED 
ONLY THE I = 20 TRACE IS ACCURATE IN PREDICTING 
THE SYSTEM'S PROBABLILITY OF BIT ERROR. THIS IS BECAUSE 
AT I = 20, THE EFFECTS OF ADJACENT CHANNEL INTERFERENCE (ACI) 
ARE NEGLIGABLE. 




rho = 0.99; 
c = (l-rho)~2; 

I = [5 6 9 20] ; 

M = [198 164 110 48]; 

bitmatrix =[000; 0 0 1; 0 1 0; 0 1 1; 100; 101; 110; 111]; 


term8 = [-464.3344 -327.7899 -148.2912 -29.7545 

-464.3344 -327.7899 -148.2912 -29.7545 

-928.6688 -655.5798 -296.5824 -59.5090]; 


term9 = [155.3047 109.6350 49.5986 9.9519 

155.3047 109.6350 49.5986 9.9519 

310.6094 219.2701 99.1971 19.9038]; 


termlO = [-0.0200 -0.0142 -0.0064 -0.0013 

-0.0200 -0.0142 -0.0064 -0.0013 

-0.0401 -0.0283 -0.0128 -0.0026]; 


for ct = 1:4 
for i = 1:8 

for boo = [1 0] 

XI = boo*c*8.5126e+03; 

X2 = bitmatrix(i,3)*c*492.4271; 

X3 = bitmatrix(i,3)*c*4.2917e-02; 

X4 = bitmatrix(i,3)*c*l.5865e-05; 

*/. NOTE THAT THE RESULTS FOR TERM #5, APPENDIX B ARE 
*/. MISSING. SEE THE COMMENTS IN THE TITLE SECTION 
’/. ABOVE. 

X6 = boo*bitmatrix(i,3)*c*9.9411e+02; 

X7 = boo*bitmatrix(i,3)*c*7.3500e-01; 
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X8 = 0; 

elseif bitmatrix(i,1)==1 Sc bitmatrix(i,2)==0 
X8 = c*term8(l,ct); 

elseif bitmatrix(i,1)==0 Sc bitmatrix(i,2)==1 
X8 = c*term8(2,ct); 

elseif bitmatrix(i,1)==1 & bitmatrix(i,2)==1 
X8 = c*term8(3,ct); 
else 

X8 = 0; 

end 

if bitmatrix(i,3) == 0 
X9 = 0; 

elseif bitmatrix(i,i)==l Sc bitmatrix(i,2)==0 
X9 = c*term9(l,ct); 

elseif bitmatrix(i,1)==0 Sc bitmatrix(i,2)==1 
X9 = c*term9(2,ct) ; 

elseif bitmatrix(i, 1)==1 Sc bitmatrix(i,2)==1 
X9 = c*terra9(3,ct); 
else 

X9 = 0; 

end 

if bitmatrix(i,3) ==0 
X10 = 0; 

elseif bitmatrix(i, 1)==1 Sc bitmatrix(i,2)==0 
X10 = c*termlO(l,ct); 

elseif bitmatrix(i, 1)==0 Sc bitmatrix(i,2)==1 
X10 = c*terml0(2,ct); 

elseif bitmatrix(i, 1)==1 Sc bitmatrix(i,2)==1 
X10 = c*terml0(3,ct) ; 
else 

X10 = 0; 

end 

INTERSUM = X1+X2+X3+X4+X6+X7+X8+X9+X10; 
if boo==l 

XlBETA(ct,i) = INTERSUM; 
else 

XOBETA(ct,i) = INTERSUM; 

end 

end 

end 

XlMIN(ct) = min(XlBETA(ct,:)); 

XOMAX(ct) = max(XOBETA(ct,:)); 

VT(ct) = (XOMAX(ct) + X1MIN(ct))/2; 
end 

RPSQR_TDIVNO_DB = 10:.01:20; 

RPSQR.TDIVNO = 10.*(RPSQR_TDIVN0_DB*0.10); 

SINGCHAN = 0.5*erfc(RPSQR_TDIVN0/8~0.5); 
for ct = 1:4 

PE(ct,1:1001) = zeros(l,1001); 
for i = 1:8 



PE(ct,:)=PE(ct,:)+0.25*erfc((RPSQR_TDIVN0/2~0.5)*(XlBETA(ct,i)-VT(ct))).. 
+0.25*erfc((RPSQR_TDIVN0/2~0.5)*(VT(ct)-X0BETA(ct,i))); 

end 

end 

PEFINAL = PE/8; 
figure(1) 

semilogy(RPSQR_TDIVNO_DB,SINGCHAN,'—',RPSQR_TDIVNO_DB,PEFINAL(1,:),... 

RPSQR_TDIVNO_DB,PEFINAL(2,:),RPSQR_TDIVNO_DB,PEFINAL(3. 
RPSQR_TDIVNO_DB,PEFINAL(4,:)) 
xlabel('Z (dB)O; 
ylabel('Pb'); 
axis([l0 17 10~(-15) 1]) 
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THESIS COMPUTER WORK 


7. 

7. 

7. 

7. 

7. 

7 . 

7 . 

7. 

7. 

7. 

7. 

7. 

7. 

7. 

7. 

7. 

7. 

7. 

7. 


GENERATION OF GRAPHS 
FOR 

betatau = 1500 

DATE: 25 OCTOBER 1994 

COMMENTS: THIS GRAPHING PROGRAM DOES NOT CONTAIN THE RESULTS FOR 
TERM #5, APPENDIX B DUE TO AN INABILITY TO ARRIVE AT A 
SOLUTION. THERE WERE MEMORY PROBLEMS WHEN A MATRIX SOLUTION 
WAS TRIED TO SPEED CALCULATION OF THIS TERM, AND THE 
TIME REQUIRED TO CALCULATE USING SIMPLE LOOPS WAS 
PROHIBITIVE. THUS, ON THE RESULTING GRAPH PRODUCED 
ONLY THE I = 20 TRACE IS ACCURATE IN PREDICTING 
THE SYSTEM'S PROBABLILITY OF BIT ERROR. THIS IS BECAUSE 
AT I = 20, THE EFFECTS OF ADJACENT CHANNEL INTERFERENCE (ACI) 
ARE NEGLIGABLE. 




rho = 0.99; 
c = (l-rho)~2; 

I = [7 9 12 20]; 

M = [212 164 124 74]; 

bitmatrix = [000; 001; 0 10; 0 1 1; 100; 10 1; 1 10; 1 1 1]; 


term8 = [-3.5253e+02 -2.1854e+02 -1.2475e+02 -4.5096e+01 

-3.5253e+02 -2.1854e+02 -1.2475e+02 -4.5096e+01 

-7.0505e+02 -4.3708e+02 -2.4950e+02 -9.0192e+01]; 


term9 = [117.9019 73.0900 41.7221 15.0822 
117.9019 73.0900 41.7221 15.0822 
235.8037 146.1800 83.4441 30.1644]; 


termlO = [-1.0000e-04 -0.6199e-04 -0.3539e-04 -0.1279e-04 
-1.0000e-04 -0.6199e-04 -0.3539e-04 -0.1279e-04 
-0.2000e-03 -0.1240e-03 -0.0708e-03 -0.0256e-03]; 


for ct = 1:4 
for i = 1:8 

for boo = [1 0] 

XI = boo*c*9.0083e+03; 

X2 = bitmatrix(i,3)*c*328.3413; 

X3 = bitmatrix(i,3)*c*l.8815e-04; 

X4 = bitmatrix(i,3)*c*7.2482e-10; 

7. NOTE THAT THE RESULTS FOR TERM #5, APPENDIX B ARE 
7, MISSING. SEE THE COMMENTS IN THE TITLE SECTION 
7. ABOVE. 

X6 = boo*bitmatrix(i,3)*c*6.6331e+02; 

X7 = boo*bitmatrix(i,3)*c*5.1105e-03; 
if boo == 0 







X8 = 0; 

elseif bitmatrixCi,1)==1 k 
X8 = c*term8(l,ct); 
elseif bitmatrixCi,1)==0 k 
X8 = c*term8(2,ct); 
elseif bitmatrix(i,1)==1 k 
X8 = c*term8(3,ct); 
else 

X8 = 0; 

end 

if bitmatrixCi,3) == 0 
X9 = 0; 

elseif bitmatrixCi,1)==1 k 
X9 = c*term9(l,ct); 
elseif bitmatrix(i,1)==0 k 
X9 = c*term9(2,ct); 
elseif bitmatrix(i,1)==1 k 
X9 = c*term9(3,ct); 
else 


bitmatrix(i,2)==0 
bitmatrix(i,2)==1 
bitmatrixCi,2)==1 


bitmatrix(i,2)==0 
bitmatrixCi,2)==1 
bitmatrix(i,2)==1 


X9 = 0; 

end 

if bitmatrix(i,3) ==0 
X10 = 0; 

elseif bitmatrix(i,1)==1 k bitmatrix(i,2)==0 
X10 = c*termlO(l ,ct); 

elseif bitmatrixCi,1)==0 k bitmatrix(i,2)==1 
X10 = c*terml0(2 ,ct); 

elseif bitmatrix(i,1)==1 k bitmatrixCi,2)==1 
X10 = c*termlOC3,ct); 
else 

X10 = 0; 

end 

INTERSUM = X1+X2+X3+X4+X6+X7+X8+X9+X10; 
if boo==l 

XlBETACct.i) = INTERSUM; 
else 


XOBETACct.i) = INTERSUM; 

end 


end 


end 

XlMINCct) = minCXIBETACct,:)); 
XOMAXCct) = maxCXOBETACct,:)); 
VTCct) = CXOMAX(ct) + XlMINCct))/2; 
end 


RPSQR_TDIVNO_DB = 10:.01:20; 

RPSQR.TDIVNO = 10.“(RPSQR_TDIVN0_DB*0.10); 

SINGCHAN = 0.5*erfcCRPSQR_TDIVN0/8"0.5); 
for ct = 1:4 

PECct,1:1001) = zeros Cl ,1001); 
for i = 1:8 

PECct,:)=PECct,:)+0.25*erfcC(RPSQR_TDIVN0/2~0.5)*CXlBETA(ct,i)-VT(ct))).. 
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+0.25*erf c ( (RPSQR_TDIVN0/2''0.5) *(VT(ct) -X0BETA(ct, i) ) ) ; 

end 

end 

PEFINAL = PE/8; 
figure(l) 

semilogy(RPSQR_TDIVNO_DB,SINGCHAN, ’ — ’ ,RPSqR_TDIVNO_DB,PEFINAL(1,:), . .. 

RPSQR_TDIVNO_DB,PEFINAL(2,:),RPSQR_TDIVNO_DB,PEFINAL(3,:) ,. . . 
RPSQR.TDIVNO.DB,PEFINAL(4,:)) 
xlabel('Z (dB)’); 
ylabelC'Pb’); 
axis( [10 17 10“(-15) 1]) 
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THESIS COMPUTER WORK 


*/. 

•/. 

y. 

•/. 
y. 
y. 
y. 
y. 
y. 
y. 
y. 

•/. 
y. 

•/. 
y. 
y. 
y. 
y. 
y. 
y. 
y.: 

rho = 0.99; 
c = (l-rho)~2; 

I = [8 9 12 20] ; 

M = [248 220 164 98]; 
bitmatrix = [000; 0 0 1; 0 1 0; Oil; 100; 101; 110; 111]; 

term8 = [-351.8931 -283.0055 -163.9044 -60.1542 
-351.8931 -283.0055 -163.9044 -60.1542 
-703.7863 -566.0110 -327.8088 -120.3084]; 

term9 = [91.5258 73.0261 41.6758 15.0544 
91.5258 73.0261 41.6758 15.0544 
183.0515 146.0523 83.3516 30.1088]; 

termlO = [-6.5584e-07 -5.2745e-07 -3.0548e-07 -1.1211e-07 
-6.5584e-07 -5.2745e-07 -3.0548e-07 -1.1211e-07 
-1.3117e-06 -1.0549e-06 -6.1096e-07 -2.2423e-07]; 

for ct = 1:4 
for i = 1:8 

for boo = [1 0] 

XI = boo*c*9.2563e+03; 

X2 = bitmatrix(i,3)*c*246.2563; 

X3 = bitmatrix(i,3)*c*9.2720e-07; 

X4 = bitmatrix(i,3)*c*3.2152e-14; 

’/. NOTE THAT THE RESULTS FOR TERM #5, APPENDIX B ARE 
*/. MISSING. SEE THE COMMENTS IN THE TITLE SECTION 
'/. ABOVE. 

X6 = boo*bitmatrix(i,3)*c*4.9749e+02; 

X7 = boo*bitmatrix(i,3)*c*3.4503e-05; 


GENERATION OF GRAPHS 
FOR 

betatau = 2000 

DATE: 15 NOVEMBER 1994 

COMMENTS: THIS GRAPHING PROGRAM DOES NOT CONTAIN THE RESULTS FOR 
TERM #5, APPENDIX B DUE TO AN INABILITY TO ARRIVE AT A 
SOLUTION. THERE WERE MEMORY PROBLEMS WHEN A MATRIX SOLUTION 
WAS TRIED TO SPEED CALCULATION OF THIS TERM, AND THE 
TIME REQUIRED TO CALCULATE USING SIMPLE LOOPS WAS 
PROHIBITIVE. THUS, ON THE RESULTING GRAPH PRODUCED 
ONLY THE I = 20 TRACE IS ACCURATE IN PREDICTING 
THE SYSTEM'S PROBABLILITY OF BIT ERROR. THIS IS BECAUSE 
AT I = 20, THE EFFECTS OF ADJACENT CHANNEL INTERFERENCE (ACI) 
ARE NEGLIGABLE. 
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X8 = 0; 

elseif bitmatrix(i,l)==l k bitmatrix(i,2)==0 
X8 = c*term8(l,ct); 

elseif bitmatrix(i,l)==0 k bitmatrix(i,2)==1 
X8 = c*term8(2,ct); 

elseif bitmatrix(i,l)==l k bitmatrix(i,2)==1 
X8 = c*term8(3,ct); 
else 

X8 = 0; 

end 

if bitmatrix(i,3) == 0 
X9 = 0; 

elseif bitmatrix(i, 1)“1 k bitmatrix(i,2)==0 
X9 = c*term9(l,ct); 

elseif bitmatrix(i,1)==0 k bitmatrix(i,2)==1 
X9 = c*term9(2,ct); 

elseif bitmatrix(i,l)==l k bitmatrix(i,2)==l 
X9 = c*term9(3,ct); 
else 

X9 = 0; 

end 

if bitmatrix(i,3) ==0 
X10 = 0; 

elseif bitmatrix(i,l)==l k bitmatrix(i,2)==0 
X10 = c*termlO(l,ct); 

elseif bitmatrix(i,l)==0 & bitmatrix(i,2)==1 
X10 = c*terml0(2,ct); 

elseif bitmatrix(i,l)==l k bitmatrix(i,2)==1 
X10 = c*terml0(3,ct); 
else 

X10 = 0; 

end 

INTERSUM = X1+X2+X3+X4+X6+X7+X8+X9+X10; 
if boo==l 

XlBETA(ct,i) = INTERSUM; 

©1 s © 

X0BETA(ct,i) = INTERSUM; 

end 

end 

end 

XlMIN(ct) = min(XlBETA(ct,:)); 

XOMAX(ct) = max(XOBETA(ct,:)); 

VT(ct) = (XOMAX(ct) + XlMIN(ct))/2; 
end 

RPSQR_TDIVNO_DB = 10:.01:20; 

RPSQR.TDIVNO = 10.“(RPSQR_TDIVN0_DB*O.10); 

SINGCHAN = 0.5*erfc(RPSQR_TDIVN0/8"0.5); 
for ct = 1:4 

PE(ct,1:1001) = zeros(l,1001); 
for i = 1:8 




PE(ct,:)=PE(ct,:)+0.25*erfc((RPSQR_TDIVN0/2~0.5)*(X1BETA(ct,i)-VT(ct))).. 
+0.25*erfc((RPSQR_TDIVN0/2~0.5)*(VT(ct)-X0BETA(ct,i))); 

end 

end 

PEFINAL = PE/8; 
figure(1) 

semilogy(RPSQR_TDIVNO_DB,SINGCHAN,’—',RPSQR_TDIVNQ_DB,PEFINAL(1. 

RPSQR_TDIVNO_DB,PEFINAL(2,:),RPSQR_TDIVNO_DB,PEFINAL(3,:),... 
RPSQR_TDIVNO_DB,PEFINAL(4,:)) 
xlabel('Z (dB)'); 
ylabelC'Pb'); 
axis([!0 17 lO''(—15) 1]) 
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